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FOREWORD 


This  is  a  final  report  covering  three  years  of  research 
on  the  thermal- physical  parameters  of  materials.  The  purpose 
of  the  research  was  o  establish,  on  the  basis  of  theoretical 
studies,  a  simplified,  mathematical  model  of  the  fabric-skin 
system  using  heat  flow  equations  and  to  solve  these  equations 
by  means  of  programming  on  a  Control  Data  G-15D  digital  com¬ 
puter.  The  results  of  the  calculations,  v.hen  compared  with 
literature  values  and  with  measurements  made  in  the  laboratory 
on  fabric  systems  as  to  surface  temperature  and  heat  transfer, 
will  provide  better  insight  into  the  chemical  degradation 
processes  taking  place  in  materials  exposed  io  high  intensity 
thermal  energy  pulses.  This  will  lead  eventually  to  new  mate¬ 
rials  affording  higher  thermal  protection  than  that  provided 
by  present  conventional  materials. 

The  work  was  begun  by  the  Aerospace  Research  Department 
of  General  Mills'  Electronics  Division  under  Contract  No. 
DA-19-129-QM-2998.  On  11  September  1963,  the  Aerospace 
Engineering  and  Research  Departments  of  General  Mills'  Elec¬ 
tronics  Division  were  purchased  by  Litton  Industries  and  are 
now  known  as  the  Applied  Science  Division  (ASD),  Litton 
Systems,  Inc.  Since  the  sale  transferred  all  facilities, 
personnel  and  capabilities  involved  in  this  program,  the 
research  continued  without  interruption  under  the  same  con¬ 
tract.  The  first  annual  report  was  issued  by  General  Mills; 
the  second  annual  report  and  this  final  report  by  the  Applied 
Science  Division,  Litton  Systems,  Inc.  The  work  was  directed 
during  the  first  year  by  Hr.  Sheldon  Steinberg,  and  in  the 
second  and  third  years  by  Mr.  R.  E,  Larson  with  Mr.  A.  R.  Xydd 
assisting. 

The  contract  was  initiated  under  Project  No.  1K0~2/;401-A113 
and  was  administered  unde1'  the  Materials  Research  Branch  of  the 
Clothing  and  Organic  Materials  Division  of  the  U.  S.  Army  Natick 
Laboratories,  with  Mr.  David  Feldman  acting  as  Project  Officer 
and  Dr.  William  E.  C.  Yelland  as  Alternate  Project  Officer. 
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ABSTRACT 


This  is  the  final  report  of  a  three-year  program  of  research  on 
the  thermal -physical  parameters  of  materials  necessary  for  protection 
against  intense  thermal  radiation.  It  briefly  summarizes  the  theoretical 
and  experimental  studies  of  the  first  two  years,  and  includes  a  more 
comprehensive  discussion  of  the  work  performed  during  the  final  year 
of  the  contract. 

In  the  theoretical  studies,  a  simplified  model  of  the  fabric -skin 
system  was  chosen  and  the  heat  flow  equations  were  solved  using  an 
explicit  finite  difference  calculation  procedure  programmed  for  solu¬ 
tion  on  a  Control  Data  G-15D  digital  computer.  Several  sample  cal¬ 
culations  were  compared  with  results  existing  in  the  literature  and  a 
limited  parameter  study  was  performed  with  variation  of  the  fabric 
surface  heat  transfer  coefficients. 

Fabric  test  samples  were  irradiated  in  an  arc  image  furnace  test 
facility  and  the  heating  and  degradation  processes  were  qualitatively 
and  quantitatively  studied  with  a  Mach-Zehnder  interferometer. 

This  instrument  provided  fabric  surface  temperature  data,  heat 
transfer  coefficient  information,  and  a  better  insight  into  the  chemical 
processes  taking  place. 
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1.  INTRODUCTION 

During  recent  years,  considerable  research  effort  has  been  ex¬ 
pended  on  studying  the  effects  of  intense  thermal  radiation,  such  as  that 
from  released  nuclear  energy,  on  military  clothing  fabrics.  A  variety 
of  experimental  techniques  and  procedures  have  been  used;  screening 
tests  on  all  types  of  fabric  s  hav  ?  been  carried  out  for  various  thermal 
radiation  exposures.  Several  theoretical  models  ha’  e  been  postulated 
and  studied.  The  influence  of  such  effects  as  fabric  diathermancy, 
moisture  content,  surface  heat  losses,  and  time  dependent  heat  inputs 
have  been  investigated.  But  at  the  present  time,  a  basic  insight  into 
the  detailed  processes  involved  in  protection  against  intense  thermal 
radiation  has  not  been  obtained. 

The  general  scope  of  the  study  repo rted  herein  is  to  determine  the 
physical  properties  required  in  a  new  material  that  could  be  used  as  a 
clothing  fabric  for  protection  against  intense  thermal  radiation.  Op¬ 
timum  values  are  to  be  determined  for  these  properties  on  the  basis 
of  a  theoretical  analysis  and  parameter  study.  To  complement  the 
theoretical  study,  an  experimental  investigation  is  required  to  perfect 
procedures  and  techniques  to  measure  thermal,  physical,  and  optical 
properties  of  existing  fabric  materials. 

To  accomplish  these  objectives  a  five-year  program  of  research 
was  originally  outlined  -  three  years  of  which  have  been  completed. 
Reductions  in  DASA  funds  curtailed  the  completion  of  the  entire  five- 
year  research  program  This  tinal  report  summarizes  our  accom¬ 
plishments  of  the  above  specific  objectives  to  the  extent  possible 

During  the  first  year  (as  reported  in  Reference  1),  a  simplified 
mathematical  model  was  chosen,  and  the  heat-flow  equation  was  solved 
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using  an  explicit  finite  difference  technique.  Considered  in  the  analysis 
were  a  time -dependent  heat  input,  fabric  diathermancy,  and  nonlinear 
radiation  and  convection  boundary  conditions. 

Not  included  in  this  theoretical  model  are  various  physicochemical 
effects  which  occur  in  a  real  fabric  as  the  temperature  increases.  These 
can  include  thermal  property  variations  and  changes  in  the  fabric  optical 
properties.  In  addition,  heat  liberation  or  absorption  processes  which 
accompany  some  chemical  reactions  effectively  create  heat  sources  or 
sinks  in  the  fabric  which  must  be  accounted  for. 

Several  sample  calculations  were  performed  for  comparison  with 
existing  theories.  A  concurrent  experimental  investigation  was  initiated 
to  provide  much  needed  information  pertaining  to  the  fabric  heat-flow 
problem.  A  major  purpose  of  this  investigation  was  to  provide  infor¬ 
mation  pertaining  to  the  heat-flow  processes  at  the  front  surface  of  the 
fabric  and  in  the  air-gap  region.  In  addition,  further  information  per¬ 
taining  to  the  thermal,  optical,  and  physical  properties  of  various 
fabrics  were  also  needed. 

Theoretical  considerations  were  given  to  the  flash  discharge  method 
for  determining  thermal  diffusnuty  of  diathermanous  materials.  An 
equation  describing  the  temperature  time  variation  of  a  specimen  rear 
surface  for  energy  absorption  in  depth,  as  opposed  to  absorption  at  the 
front  surface,  was  derived.  An  extension  of  this  developed  equation 
was  derived  to  consider  front  surface  heat  losses. 

Laboratory  test  experiments  with  the  flash  tube  apparatus  were  con¬ 
ducted  and  thermal  property  data  were  obtained  for  several  opaque  con¬ 
ductors  (metals)  and  for  opaque  and  diathermanous  non-conductors 
(plastics).  The  temperature  time  histories  compared  favorably  with 
theoretical  predicted  values. 

During  the  second  year  of  this  program  theoretical  procedures  were 
extended  and  a  computational  study  was  performed  to  better  understand  the 
relative  significance  of  various  parameters  isolated  auring  the  first  year's 


efforts.  Explicit  finite  difference  calculations  were  performed  on  a 
Control  Data  G-l  5D  digital  computer.  (Refer  to  Reference  2). 

Emphasis  on  experimental  efforts  was  increased  during  the  second 
year.  An  arc  image  furnace,  thermal  test  chamber,  and  optical 
visualization  equipment  were  set  up  to  study  heating  effects  of  fabrics 
exposed  to  square  wave  radiant  energy  heat  pulses.  Techniques 
utilizing  a  Mach-Zehnder  Interferometer  were  developed  to  measure 
fabric  surface  temperature.  The  interference  events  were  recorded 
with  motion  picture  cameras  allowing  accurate  descriptions  of  tempera¬ 
ture-time  variations  for  the  transient  heating  phenomena.  A  6 -ounce 
cotton  poplin,  an  8.  8-ounce  cotton  sateen  and  an  8-ounce  wool  (OG107) 
were  selected  for  considerable  experimental  testing. 

Complete  description  of  equipment  and  experimental  procedures, 
with  major  emphasis  given  the  interferometer,  and  results  of  calibra¬ 
tion  tests  and  preliminary  temperature  time  variation  recordings  are 
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included  in  the  Second  Annual  Report. 

During  the  third  and  final  year  major  emphasis  was  given  experi¬ 
mental  investigations  using  interferometer  techniques.  Experiments  were 
conducted  on  fabric -ai r space -skin  models  and  temperatur  e -time  variations 
compared  with  the  theoretical  ncn-dimensional  heat  transfer  response 
characteristics  for  the  model  developed  during  the  first  year.  Results 
from  interference  temperature  recordings  are  presented  in  terms  of 
various  heat  transfer  parameters. 

Further  interferometer  techniques  were  developed  to  enable  de¬ 
termination  of  degradation  gases  and  concentrations  as  a  function  of 
temperature  and  time.  The  necessity  for  considering  chemical  proper¬ 
ties  evolved  from  requirements  f^r  obtaining  accurate  uterferometer 
temperature -heat  transfer  data.  Refractive  index  variations  due  to  the 
presence  of  unknown  degradation  gases  were  studied.  To  analyze  these 
variations  a  multi -wavelength  technique  and  a  multi -environmental  test 
technique  were  considered  m  detail.  Expressions  for  interference 
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fringe  shift  displacements  in  terms  of  temperature,  gas  molecular  refrac- 
tivity,  and  gas  concentration  were  developed  for  both  techniques  and  pro¬ 
grammed  for  solution  on  a  Control  Data  1604  Digital  Computer.  Com¬ 
prehensive  tests  were  made  on  fabric  systems  in  separate  environments 

of  air,  helium,  nitrogen,  argon,  and  C07.  For  water  saturated  fabrics 

L  2 

under  heat  loads  of  approximately  6  cal/cm  -sec  concentrations  of  vapor 
at  the  fabric  surface  were  determined  as  a  function  of  temperature  and 
time. 

Fabric  front  and  rear  surface  heat  transfer  coefficients  were 
measured  both  for  steady  state  and  transient  heat  flow  conditions.  These 
tests  showed  that  extreme  caution  must  be  shown  in  properly  evaluating 
heat  transfer  coefficients  for  conditions  of  rapidly  changing  surface 
temperatures. 


2.  THEORETICAL  STUDIES 


2.  I  Previous  Theoretical  Studies 

Solutions  to  problems  involving  the  conduction  of  heat  in  solids  can 
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be  found  in  such  textbooks  as  Carslaw  and  Jaeger.  They  present 
several  solutions  for  semi -infinite  slabs  in  which  the  use  of  the  Laplace 
transform  makes  the  solution  relatively  simple.  Boundary  conditions 
involving  radiation  (convection  and  linearized  thermal  radiation)  into 
a  medium  with  temperature  varying  with  time  are  reduced  to  the  case 
of  constant  environment  conditions  by  Duhamel  s  theorem.  Exact 
solutions  for  nonlinear  boundary  conditions  or  when  properties  are  a 
function  of  temperature  are  not  presented. 

The  fabric -skin  heat  transfer  problem  was  studied  by  Chen  and 
4 

Jensen,  who  emphasized  a  numerical  graphical  technique  patterned 
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after  Schn  dt  s  original  method  (with  modifications  of  initial  condi¬ 
tions  and  method  of  boundary  subdivision),  linearized  radiation  and 
convection  losses  at  the  fabric  surfaces  and  combined  them  into 
single  heat-loss  terms.  Their  method  was  developed  for  square - 
wave  radiation  heat  input,  although  they  showed  how  it  could  be  ap¬ 
plied  to  a  time -viependeut  heat  pulse.  Fabric  diathermancy  was 
considered,  and  optical  and  physical  properties  were  assumed  con¬ 
stant.  The  results  were  presented  in  a  general  nondimensional  form 
for  the  case  of  fabric  opacity,  and  a  "thr ee -cu rve "  method  was  pre¬ 
sented  which  provided  fairly  accurate  results  for  various  fabric 
diathermancies. 

In  another  report,  Hottel,  et  al.,^  extended  the  work  of  Reference  4 
to  include  the  case  of  a  moist  cloth -skin  system  and  discussed  the 
computat icna!  methods  in  great  d<  tail.  A  method  was  developed  for 
square-wave  radiation  impinging  on  opaque  fabric,  but  the  possibilities 
of  extending  the  calculations  to  include  diathermanous  fabric  and  a 
time -variant  heat  input  were  mentioned.  Experiments  were  carried 
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out  to  substantiate  both  the  dry-  and  moist-cloth  models,  and  results 
agreed  fairly  well  with  the  calculations.  The  authors  reported  that 
the  greatest  uncertainties  existed  in  the  evaluation  of  heat  and  mass 
transfer  coefficients  aud  water  vapor  diffusivity. 

7 

In  a  research  program  at  Armour,  an  analytical  study  of  single 
and  composite  slab  systems  irradiated  by  square -wave  and  variable  in¬ 
tensity  radiation  was  reported,,  Working  charts  to  allow  prediction  of 
the  degree  of  burn  were  presented  for  a  composite  slab  exposed  to 
variable  intensity  radiation,  neglecting  the  effects  of  surface  re- 
radiation. 

Calculations  that  included  the  re-radiation  term  were  made  for  a 
protecting  material  of  finite  thickness  in  perfect  thermal  contact  with 
an  infinite  slab  (the  skin).  This  nonlinear  boundary  condition  was 
handled  by  formulating  a  function  which  described  the  temperature - 
time  variation  at  the  outer  surface.  Two  parameters  in  the  function 
were  evaluated  by  trial  and  error. 

8 

The  A.  D.  Little  Company  has  aiso  completed  a  theoretical  and 
experimental  investigation  of  heat  transfer  tc  skin  through  a  cloth  layer. 
One  of  the  models  they  considered  was  multilayered  with  the  inner 
layer  having  a  low  thermal  <  onductivity  (air),  overlayed  first  by  a  low 
emissivity  layer,  and  finally  an  outer  layer  which  evaporates,  and  in 
so  doing,  absorbs  heat.  A  constant  property  solution  was  obtained  and 
square  wave  and  triangular  approximations  of  the  thermal  pulse  were 
assumed.  Their  first  model  assumed  no  temperature  gradients  normal 
to  the  outer  two  layers  and  a  constant  temperature  at  the  bottom  of  the 
air  space.  They  also  linearized  the  radiation  boundary  condition  term 
for  one  calculation  and  retained  .t  for  another.  An  improved  model 
made  it  possible  to  vary  the  tempi  rature  both  in  the  air  space  and  in 
the  skin  at  the  bottom  of  the  dir  space.  Their  final  model  allowed  the 
temperature  to  vary  with  depth  in  the  outer  <1  -thing  layer.  They  re¬ 
tained  the  nonlinear  radiation  boundary  condition  and  used  a  finite  dif¬ 
ference  method  of  solution.  They  did  not  consider  any  dathe  rrna.n  y 
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effects  through  the  cloth.  Their  final  model  considered  the  effective¬ 
ness  of  smoke  generation  in  reducing  heat  transfer  to  the  skin. 

The  heat  conduction  equation  has  been  solved  for  many  configura¬ 
tions  using  both  analytical  and  numerical  methods.  However,  the 
boundary  conditions  used  were  not  general  enough  to  apply  to  our 
problem. 

Additional  theoretical  studies  pertaining  to  the  present  problem 
are  included  in  Reference  1.  In  addition,  new  results  have  since 
been  reported  in  the  literature.  The  new  studies  will  not  be  dis¬ 
cussed,  however,  as  they  do  not  appreciably  advance  the  present 
state  of  the  knowledge.  The  greatest  improvment  in  the  fabric -skin 
theory  now  available  will  come  when  chemical  effects  can  be  included 
in  the  mathematical  analysis. 

2.  2  Development  of  Theoretical  Model 

The  studies  discussed  in  2.  1,  above,  and  in  the  First  Annual  Report* 
have  contributed  in  varying  degrees  towara  the  attainment  of  a  basic 
understanding  of  the  fabric -skin  heat  transfer  problem.  Various  assump¬ 
tions  have  been  made  regarding  the  proper  heat  transfer  mechanisms 
to  insert  into  the  surface  boundary  conditions.  Factors  such  as  fabric 
diathermancy,  moisture  content,  and  time -dt  pendent  heat  inputs  have 
been  considered.  But  further  work,  both  analytical  and  experimental, 
is  needed  to  better  define  the  heat  tlow  processes  within  the  fabric  and 
at  the  surfaces.  Very  little  information  is  available  regarding  actual 
temperature  distributions  through  the  fabric  during  a  heat  pulse  because 
of  the  experimental  difficulties  involved  Usually  the  effectiveness  of 
a  protective  fabric  layer  is  judged  by  the  heating  effects  experienced  by 
a  skm  simulant  v.  r  calorimeter  placed  behind  the  fabric. 

Under  conditions  of  high  radiant  loading  and  associated  fabric 
temperature*,,  chemical  r»  actions  influence  the  heat  flow  process  to  a 
considerable  extent.  The  mtmn.te  interrelationship  between  heat  transfer 


and  chemical  reactions  (such  as  ablation,  sublimation,  oxidation  and 
reduction,  and  degradat’.on)  have  long  been  recognized  and  studied  in 
great  detail  in  nose  cone  re-entry  problems.  To  our  knowledge,  the 
effects  of  chemical  reactions  in  the  fabric-skin  heat  transfer  problem 
have  not  been  included  in  any  analysis.  These  reactions  complicate 
any  mathematical  or  experimental  studies  because  they  can  give  rise 
to  heat  absorption  or  generation,  physical  and  optical  property  variation 
and  surface  degradation  that  greatly  influences  both  absorbed  and  re¬ 
radiated  energy. 

Fabric  moisture  content  is  important  (as  shown  in  Reference  6), 
and  additional  analytical  and  experimental  efforts  are  needed  to  under¬ 
stand  the  transfer  processes  involved.  Heat  transfer  processes  be¬ 
tween  the  fabric  rear  surface  and  skin,  and  convection  and  reradiation 
at  the  front  surface  should  be  investigated  further.  Additional  information 
is  needed  regarding  the  temperature  dependence  of  a  fabric's  physical 
and  optical  properties. 

The  model  initially  chosen  for  the  present  investigation  is  fairly 
simple;  it  does  not  allow  inclusion  of  all  the  basic  mechanisms  and 
processes  mentioned  above.  It  will,  however,  be  improved  as  fur¬ 
ther  insight  is  obtained  from  our  experiments  and  parameter  studies. 

More  of  the  fabric’s  physical  and  optical  properties  will  appear  during 
the  formulation  of  the  heat  flow  equations.  A  meaningful  insight  into 
the  problem  can  only  be  obtained  1*  the  relative  importance  of  parameters 
formulated  from  these  proper  tuns  is  <ht<  rmmed.  Initial  emphasis  will 
be  placed  on  the  fabric  itself,  rather  than  on  the  heat  transfer  and 
degree  of  burn  experienced  by  the  skm. 

For  our  problem,  two  one -dimensional  models  are  postulated  (see 
Figure  I).  Both  involve  infinite  flat  surfaces  of  finite  thickness.  In  the 
first,  the  flat  plate  is  separated  from  another  structure,  the  skin,  bv  an 
air  space;  in  the  seem  d  model,  the  fabric  is  in  intimate  contact  with  the 
skin.  In  the  actual  case,  combinations  of  these  two  limiting  conditions 
.vuuic!  probably  be  obtained. 
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Figure  1.  Heat  Transfer  Models 
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The  first  model  actually  involves  solution  of  a  three -layer  problem, 
the  fabric  and  air  space  of  finite  thickness  and  the  skin  layer  extending 
to  infinity.  We  have  assumed  that  a  heat  input  varying  with  time  is  im¬ 
pinging  on  the  fabric  outer  surface.  A  certain  portion  will  be  reflected 
and  the  rest  absorbed  and  diffused  into  the  fabric.  Chemical  reactions 
will  not  be  considered ,so  no  exothermic  or  endothermic  heat  processes 
will  take  place  and  the  surface  properties  will  remain  constant. 

Of  the  energy  absorbed,  a  portion  will  be  reradiated  and  convected 
from  the  outer  surface.  The  rest  will  diffuse  into  the  fabric  and  make 
its  presence  felt  by  raising  the  back  surface  temperature.  Depending 
on  the  configuration,  heat  transfer  will  then  take  place  to  the  skin  by 
processes  of  conduction,  connection  and  radiation.  The  degree  of  burn 
experienced  by  the  cutaneous  skin  layer  will  be  a  function  of  this  heat 
flux. 


2.  3  Heat  Flow  Equations  and  Boundary  Conditions 

The  problem  of  a  time -dependent  radiant -heat  input  into  a  fabric - 
skin  system  is  extremely  complex  and  practically  unsolvable  if  rigid 
adherence  to  the  true  physical  situation  is  shown.  To  simplify  the 
problem,  the  following  assumptu  ns  have  been  made  for  the  first 
model ; 

1)  The  fabri  c  -a ’  r space -ski  n  system  is  assumed  to  be 
one  dimensional. 

2)  Physical  and  optical  properties  of  the  system  re¬ 
main  constant  during  the  radiant  exposure. 

3)  No  chemical  reactions  arc  taking  place, 

4)  The  fabric  is  only  s'lght'y  duithe rmanous  to 
thermal  radiation,  and  the  skin  is  completely 
opaque. 

5}  The  fabric  has  .no  moisture  eoutent. 

6)  The  fabric  is  treated  as  a  homogenous  solid. 
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Heat  losses  from  the  front  surface  of  the  fabric  consist  of  free 
convection  and  reradiation.  Both  terms  appear  in  the  front  surface 
boundary  condition,  and  no  attempt  is  made  to  linearize  them  or  to 
group  them  into  a  single  linearized  heat-loss  term.  Similarly,  at  the 
back  surface  of  the  fabric,  radiation  interchange  is  allowed  to  take 
place  between  the  fabric  and  the  skin,  and  heat  transfer  by  conduction 
or  free  convection  through  the  air  space  is  also  assumed.  This  allows 
explicit  inclusion  of  fabric  emicsivity  and  absorptivity  in  the  equations, 
and  the  heat-transfer  processes  and  coefficients  can  be  studied 
separately.  If  a  diathermanous  material  *s  assumed,  impinging 
short-wave  radiation  is  absorbed  in  depth  across  the  fabric,  and 
the  one -dimensional,  heat-flow  equation  must  include  a  heat  source 
term.  Assuming  absorption  according  to  the  Beer -Lambert  lav/,  the 
equation  is 
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where  1  =  i  (1  -  In  the  skin,  which  is  assumed  to  be  completely 

opaque  to  thermal  radiation,  the  equation  without  the  heat  source  term 
is 
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with  the  assumption  of  constant  properties  in  both  equations. 
The  boundary  conditions  are  as  fellows, 
l)  At  the  outer  surface  of  the  fabric 
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where  the  first  three  terms  are  the  conduction,  convection,  and  re- 
radiation  at  the  surface,  and  the  last  is  the  radiation  input  term.  From  a 
strictly  mathematical  viewpoint,  the  diathermancy  factor  D  is  either 


0  or  1.  Inclusion  of  the  diathermancy  term  in  Equation  (1}  with  a 
defined  value  of  the  extinction  coefficient  requires  selection  of  the 
former.  If  the  material  is  completely  opaque,  the  diathermancy  term 
vanishes  as  in  Equation  (2).  and  D  becomes  1  as  is  usually  assumed 
for  electrically  conductive  materials.  When  the  equations  are  cast 
into  a  finite -difference  form,  it  will  be  shown  that  selection  of  a  D 
other  than  0  or  1  is  physically  realistic,  even  though  doing  so  is 
mathematically  improper  when  a  differential  representation  is  used. 

2)  At  the  inner  surface  of  the  fabric 
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where  the  second  term  represents  heat  transfer  by  conduction  or  free 
convection  through  the  air  gap,  and  the  last  term  is  the  radiation  inter¬ 
change  between  fabric  and  skin.  Further  discussion  of  the  effective  air- 

gap,  heat-transfer  coefficient  h  is  included  in  Section  3.  4. 
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3)  And  at  the  skin  layer 
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with  the  last  term  representing  the  radiation  directly  transmitted  through 
the  fabric  and  absorbed  by  the  skin.  ihe  factor  0^  is  the  absorptivity 
of  the  skin  for  short -wavelength  radiation.  Values  of  a  SQ  less  than  unity 
give  rise  to  reflected  short-wavelength  radiation  which  is  absorbed  in 
depth  in  the  fabric.  This  effect  is  neglected  m  the  present  development. 

A  detailed  discussion  of  the  above  equations,  sample  calculations, 
and  a  parameter  study  are  presented  in  Appendix  A. 


3.  EXPERIMENTAL  STUDIES 


3.1  (  j  e  n  e  r  a  1 

An  important  finding  of  the  theoretical  and  experimental  studies 
reported  in  References  1  and  2  is  the  lack  of  information  pertaining  to 
heat  transfer  processes  at  the  fabric  front  surface  and  in  the  air  space 
between  the  fabr  ic  and  skin.  Because  of  the  difficulties  associated  with 
the  measurement  of  fabric  temperatures  under  transient  conditions, 
experimental  data  ar  e  not  available  to  compare  with  theoretically  pre¬ 
dicted  venations.  Thus  new  techniques  such  as  utilization  of  the 
Mach-Zehnder  1  nterferometer  become  essential. 

The  interferometer  prov  d<  s  a:i  ms 'got  into  the  conduction  and  con¬ 
vection  heat  transfer  processes  occurring  at  both  surfaces  of  the  fabric 
and  m  the  an  spaces  between  layer  systems.  This  instrument  also 
permits  study  of  the  Vo  ious  degradation  processes  and  eventual  com¬ 
bustion  of  fabric  s. 

The  interferometer  as  well  as  the  schiieren  and  shadowgraph, 
operates  on  the  charac  ter .  Stic  phenomenon  that  density  and  pressure 
changes  in  a  gas  produce  c  flanges  in  the  refractive  index  of  that  gas. 
Whereas  the  schiieren  apparatus  s  sensitive  to  density  gradients  and 
the  shadowgraph  to  ch-tnges  in  the  density  gradients,  the  interferometer 
indicates  a  direct  change  in  dens.t.  itself.  In  a  constant  pressure 
environment  a  direct  measure  of  temperature  changes  may  be  obtained. 

The  Mach  Zehnder  nte rf e romet <•  r  has  been  used  to  great  extent  in 

aerodynamic  studies.  It  h-«  *  1..  .*n  used  less  extensive  ly  in  heat  trans- 
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fer  work,  Ec  kei  ;  u  d  S'  m  hngeo  iiitt  heuuard  were  among  the  first 

to  use  t  lie  Mii>  h  Zehnder  ’  nte  •  I  e  romett  r  to  nvestigate  heat  transfer 
phenomena..  ('tier  n  r  -  ••  e  re  e  1  stud  e  s  •  i  n  !  u  d  e  that  o  I  ( i  a  r  1  s  o  n .  ^  ( i  o !  d  - 
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A  bi  iff  review  of  some  "t  these  works  is  me  luded  in  Reference  2  and 

in  Sec  t ion  3.4  of  tins  report. 
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In  our  review  of  the  literature  we  found  no  studies  where  the 
interferometer  was  used  to  record  phenomena  surrounding  a  body 
heated  by  an  intense  radiant  source  such  as  the  arc  image  furnace. 
Furthermore  we  have  found  no  investigations  that  concern  interfero¬ 
metric  heat  transfer  studies  involving  fabrics.  In  general,  little 
attention  has  been  given  to  transient  problems  in  various  heat  flow 
processes  and  the  use  of  the  interferometer  for  recording  heat  trans¬ 
fer  or  mass  diffusion  about  a  fabric-airspace- skin  system  under  high 
thermal  loadings  represents  a  new  area  of  investigation. 

As  described  herein  the  combustion  processes  of  materials  under 
exposure  to  high  intensity  radiation  can  be  examined  with  the  inter¬ 
ference  technique.  The  possibility  of  using  the  interferometer  to 
arrive  at  gas  composition  and  concentrations  during  the  heating  stages 
and  eventual  combustion  and  buring  of  an  exposed  fabric  is  an  impor¬ 
tant  consideration  in  this  program. 

3.2  Experimental  Equipment  and  Techniques. 

3.2.  1  Thermal  Radiation  Test  Facility 

The  basic  component  of  this  facility  is  the  carbon  arc-image  fur¬ 
nace.  The  furnace  is  used  as  a  high  temperature  radiation  source  to 

simulate  solar  e.  ergy  and  thermal  radiation  from  nuclear  energy 

2 

release;  it  has  a  heat  llux  capability  of  75  cal/sec  (11  cal/cm- sec  i.  A 
U.  S.  Navy  24-inch  carbon  arc  searchlight  was  acquired  and  modified 
for  use  as  the  furnace.  The  parabulv  mirror  was  replaced  by  an  el¬ 
liptical  mirror  to  permit  focal  imaging  of  the  arc.  Power  for  the 
furnace  is  supplied  by  a  15  kw  motor  geiieratoi  set.  Figure  2  shows 
an  overall  view  of  fhe  far  ihty. 

Flux  intensities  are  measured  with  various  standard  and  miniature 
slug  calorimeters  developed  by  ASD.  Heat  pulse  initiation  ai  d  dura¬ 
tion  are  automatically  controlled  and  measuted  by  -i  glow-tube  counter 
(adjusted  to  0.01  sec  t  and  by  timing  systems  (0.001  sc  <  timing  incre¬ 
ments)  coupled  to  the  furnace  signal  key  and  motor  gear  assembly. 

<  2 
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Figure  3.  Fabric  Holder 
and  Flux  Redistributor 
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Flux  intensity  and  uniformityare  controlled  with  various  flux  re¬ 
distributor  systems  that  transform  the  circular  focal  image  of  the  arc 
into  a  rectangular  pattern  on  the  test  specimen.  This  permits  closer 
approximation  of  one -dimensional  heating  processes  and  easier  corre¬ 
lation  with  theory.  (Refer  to  Figure  3.  ) 

An  integral  part  of  this  facility  is  a  special  test  chamber  which 
can  enclose  in  a  controlled  environment  samples,  sample  holders, 
slug  calorimeters,  flux  distributors,  and  other  associated  instrumen¬ 
tation.  A  quartz  window  at  the  front  of  the  chamber  transmits  the  arc 
beam  onto  the  test  specimen,  and  two  optical  quality  viewing  windows 
allow  observation  of  the  heating  effects.  The  chamber  is  shown  in 
Figure  4. 

3.2.2  Mach-Zehnder  Interferometer  and  Associated  Equipment 

The  Mach-Zehnder  interferometer  and  associated  optical  equip¬ 
ment  are  used  to  obtain  qualitative  and  quantitative  measurements  of 
the  heating  effects  or.  fabrics.  Included  are  an  optical  bench,  Fastax 
and  still  cameras,  various  light  sources  and  other  optical  and  elec¬ 
tronic  instruments. 

The  interferometer  used  in  the  present  investigations  was  originally 
developed  for  wind  tunnel  flow  visualization  studies.  Without  modifica¬ 
tion,  this  instrument  was  used  to  study  the  he-it  transfer  processes  at 
fabric  surfaces  and  between  fabric  and  skin  surfaces  for  this  present 
study.  The  opt:*.  1 1  elements  of  the  interft  ron.eter  which  are  mounted 
on  a  2  -  1  ^2-c  entimetor  Urn  k  aluminum  plate,  ar>  3  cm  m  diameter  and 
provide  an  elliptical  viewing  area  with  dimension. s  of  approximately 
3  cm  x  5  cm. 

A  schematic  diagram  of  the  interferometer  is  shown  in  Figure  5 
and  Figure  2  also  shows  the  interferometer  and  some  of  the  photo¬ 
graphic  equipment. 
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3.3  Interferometric  Photographic  Studies 

3.3.1  Single  Environment  Tests 

3.3.  1.  1  Still  Photographic  Measurements 

A  series  of  heating  tests  on  various  fabrics  exposed  to  approxi- 

2 

mately  5  cal/cm  -sec  radiant  heating  loads  of  various  pulse  duration 
were  performed  using  still  photographic  techniques.  The  arc  shutter 
timing  mechanism  was  set  for  a  given  pulse  duration  and  immediately 
following  the  close  of  the  arc  shutters  nterferograms  were  taken  at 
1 /400th  of  a  second  using  plus  X  fiim  (ASA  160).  Figure  6  shows 
interference  fringe  patterns  for  undisturbed  heated  and  burning  condi¬ 
tions  of  a  white  cotton  twill  weave  suede  cloth.  When  burning  occurs 
the  combustion  enhances  the  free  convection  boundary  layer  thickness 
considerably  as  shown  by  the  fringe  bulging  in  Figure  6(c).  The  flame 
front  extends  well  into  the  exit  opening  of  the  flux  redistributor  and 
distortion  of  the  burning  patterns  are  clearly  shown. 

From  Figure  6(b);  the  temperature  field  distributions  on  both  sides 
of  the  fabric  have  been  determined  and  plotted  in  Figure  7.  The  con¬ 
stant  temperature  bands  are  conveniently  measured  fringe  shift  divi¬ 
sions  and  account  for  the  particular  temperature  readings  shown.  The 
boundary  thickness  temperature  distribution  for  the  center  of  heated 
area  is  shown  in  the  scaled  insert  to  the  right,  evidencing  the  higher 
temperature  of  the  front  surface,  ir.  sunsequent  motion  picture  se¬ 
quences  from  which  interference  temperature  data  were  extracted,  the 
differences  from  the  fron'  to  rear  surfaces  were  found  to  be  considera¬ 
bly  greater  for  8.8  oz  cotton  sateen  than  that  indicated  in  Figure  7. 

Interferometric  recordings  were  made  at  different  instances  m  the 
heating  phase  in  the  manner  desi  nhed  abo\e  for  a  variety  of  fabric 
materials.  For  cottons,  although  the  combustion  temperature  depends 
on  moisture  content  rate  of  heating,  fabric  purity,  etc.,  the  ignition 
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Figure  6,  Interference  Fringe  Displacement  Effects  for  Vertically 
Heated  White  Cotton  Twill -W  ;ave  Suede  Cloth 


temperature  as  indicated  from  these  interferograms  corresponds  quite 

1  8 

well  with  values  presented  by  MacHattie  in  his  measurements  of  fabric 
rear  surface  temperature  using  a  radiometric;  method.  In  Figure  8  a 
two-layer  cotton  system  (spacing  0.  6  cm)  is  shown  before  heating  (a), 
after  heating  for  1  second  (b),  and  after  burning  occurs  (c).  It  is 
interesting  to  observe  the  effects  heating  or  burning  of  a  front  layer 
has  on  a  second  layer.  When  the  first  fabric  layer  ignites,  the  flame 
front,  shown  qualitatively  in  the  interference  fringe  pattern,  appears 
squeezed  between  the  fabrics  prior  to  ignition  of  the  second  layer.  As 
the  spacing  decreases  the  flame  front  appears  further  contained  and 
does  not  develop  fully  as  at  a  free  surface.  Finally,  for  close  spacing, 
to  the  point  of  contact,  burning  envelops  both  layers  with  combustion  of 
these  layers  occurring  after  a  greater  heat  input  than  required  for 
combustion  of  a  single  test  layer  or  for  two  test  layers  of  greater 
spacing. 

It  must  be  realized  that  the  laboratory  results  of  a  multi-layer 
diathermanous  fabric  system  cannot  be  directly  related  to  field  condi¬ 
tions.  Due  to  scattering,  and  since  the  incident  energy  impinges  on  a 
relatively  small  area,  the  direct  energy  received  by  the  second  layer 
is  spread  over  a  larger  area  with  a  resultant  decrease  in  flux  density. 
Under  field  conditions  the  incident  flux  is  distributed  over  a  large  area 
and  localized  effects  do  not  o'- cur. 

3.  3,  1.2  Motion  Picture  Measurements 

Interferometric  temperatur*  time  data  were  taken  from  motion 
picture  recordings  of  the  arc  heated  fabrics.  Sequences  were  filmed 
at  speeds  of  64  frames  per  second  to  9000  frames  per  second.  For 
sequences  of  several  seconds  duration  recordings  at  64  frames  per 
second  were  adequate  for  data  analysis  purposes.  During  the  final 
heating  stages  and  combustion  instances  recordings  were  made  at  1000 
to  9000  frames  per  second  using  the  400- foot,  lb- mm  film  capacity  of  the 
Fastax  camera.  Timing  marks  at  O.OUl -second  inter  'als  permitted 
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Figure  8.  Fringe  Displacement  Effects  for  Two  Layer,  Light  Green 
Cotton  System  (Spacing  0.  6  cm) 


accurate  time  measurements  through  the  event.  Selected  frames  for 
the  heating  and  combustion  sequences  for  a  heated  cotton  specimen  "are 
included  as  Figure  9.  Because  of  the  rapid  sequence  of  events  during 
the  intermediate  phases  of  combustion,  photographs  at  0.001 -second 
increments  are  shown  starting  at  1.  705  to  1.710  seconds.  The  remain¬ 
der  of  the  photographs  are  presented  at  0.  005  -  second  increments. 

The  indicated  fabric  surface  temperature  for  a  6-oz  OG-107  cotton 
fabric,  as  a  function  of  time  after  heat  pulse  initiation,  is  shown  in 
Figure  10.  The  dashed  curve  is  a  correction  for  errors  in  the  mea¬ 
surements  attributed  to  end  effects.  These  errors,  associated  with 
bourn  ary  layer  gj  owth  at  the  edges  of  the  fabric  (increasing  r  efractive 
path  length)  was  determined  from  the  calibration  curve  (Figure  B-14 
of  Appendix  B)  for  a  nrctal  disc  0.  3  cm  thick.  Since  the  fabric  thickness 
varies  it  was  determined  that  the  correction  does  not  truly  apply  (see 
Reference  11). 

These  tests  in  general  were  carried  out  until  severe  fabric  damage 
is  visually  discerned.  Test  conditions  were  taken  at  ambient  air  and 
room  humidity  (approximately  28  percent  relative  humidity).  In  other 
>'\periments,  oven-dried  fabric  heating  response  characteristics  were 
compared  with  room  humidity  tests. 

Indicated  temperature  values  were  found  true  during  the  early  and 
intermediate  heating  stages;  thereafter  errors  due  to  the  presence  of 
degradation  gases  and  combustion  produces  were  too  great  and  analysis 
v -a  s  not  made. 

In  Figures  11  through  14,  front  and  rear  surface  indicated  tempera¬ 
ture  values  are  plotted  for  representative  military  style  fabrics.  The 
heavier  cotton  (8.  8-oz  sateen)  demonstrates  a  greater  temperature 
cliff*  'nee  from  front  to  rear  for  a  given  heat  load  over  that  shown  for 
the  6-oz  cotton  poplin  variety  as  would  be  expected.  The  temperature 
difference  indications  for  these  fabrics  are  considerably  greater  than 
that  noted  m  earlier  still-photo  recordings  where  mte rferugrarns  were 
taken  upon  close  of  the  arc-furnace  shutter.  No  corrections  for  end 
effect  were  made. 
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Figure  10.  Temperature-Time  Variation  of  Cotton  Poplin  Fabric  Surface 
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Figure  12.  Temperature-Time  Variations  of  8.8  oz  Cotton  Fabric  Surfaces 
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Figure  13.  Temperature-Time  Variation  of  Dynel  Fabric  Surfaces 
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3.3.2  Multi-Environment  Tests 


3.  3.2.  1  General  Test  Conditions 

Interferometric  recordings  of  the  fabric  heating  degradation  pro¬ 
cesses  were  made  with  the  motion  picture  camera  at  64  frames  per 
second.  High-speed  motion  picture  sequences  at  1000  frames  per 
second  or  more  were  not  needed  because  the  heating  rate  in  current 
tests  only  required  data  analysis  at  approximately  0.  10- second  inter¬ 
vals  for  sufficient  resolution. 

3.  3.  2.  2  Fabric  Tests  -  Ambient  Humidity 

The  test  procedure  involved  evacuation  of  the  thermal  test  cham¬ 
ber  and  introduction  of  a  particular  environment.  To  speed  the  test 
processes  and  the  number  of  tests  per  given  time,  a  vacuum  of  5  mm  Hg 
was  obtained  and  the  new  gas  flushed  into  the  chamber  while  the  vacuum 
pump  was  still  pumping.  It  was  anticipated  that  purging  was  not  com¬ 
plete  and  the  euvir'  ament  conceivably  contaminated  with  up  to  1/2  per¬ 
cent  air.  This  would  result  in  a  small,  but  undetermined  error  in  the 
assumed  molecular  refractivity  of  the  environmental  gas. 

The  first  series  of  tests  in  the  different  environments  were  con¬ 
ducted  on  cottons  and  wool  at  room  humidity  conditions.  Figure  15  is 
a  plot  of  the  mterferenc  e  fringe  shift  variations  with  time  (raw  data 
form)  for  a  heated  6-oz  cotton  poplin  specimen  in  a  helium  environ¬ 
ment.  Included  in  this  plot  are  rear  boundary  layer  fringe  values 
during  degradation.  The  depression  in  the  ft  mge  shif*  values  for 
helium  prior  to  complete  degradation  (from  0.2  to  0.7  second)  suggests 
that  some  small  degradation  reactions  are  occurring.  It  was  deter¬ 
mined  that  mainly  water  vapor  (rcsidu d  moisture  content)  and  air 
(entrapped  m  the  fabric  during  evacuation  procedures)  were  evidenced 
during  the  heating  cycle. 
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ayer  in  a  Helium  Environment 


The  presence  of  both  water  vapor  and  air  were  noted  also  in  the 
helium  environment  tests  on  6-oz  cotton  poplin,  8.  8-oz  cotton  sateen, 
and  on  8-oz  wool  specimens.  See  Figures  16  through  18.  The  fringe 
shift  variations  in  other  environments  do  not  readily  indicate  the  pre¬ 
sence  of  these  gases  (representative  curves  presented  in  Figures  19 
and  20).  Argon  and  nitrogen  have  molecular  refractivities  close  to 
that  of  air  and  thus  display  quite  similar  fringe  shift  curves  for  some 
tests  even  during  the  degradation  reactions  where  it  would  be  expected 
that  combustion  (in  air)  should  differ  substantially  from  the  processes 

encountered  in  an  inert  environment.  From  tests  in  CO.,  and  with 

c. 

identical  tests  in  helium  and  in  nitrogen,  an  indication  of  the  composite 
foreign  gas  molecular  refractivity  is  obtained. 

One  would  imagine  a  greater  moisture  content  or  more  air  en¬ 
trapped  in  the  heavier  8.  8-oz  cotton  sateen  over  that  of  the  light  6-oz 
cotton  and  possibly  an  even  greater  amount  of  both  constituents  in  wool 
over  that  of  the  cottons.  To  some  extent  the  fringe  values  observed  at 
the  rear  facing  surface  for  the  three  materials  substantiate  this,  wool 
having  a  negative  fringe  value  of  approximately  1.9  before  rising  posi¬ 
tively,  the  8.  8-oz  cotton  sateen  fringe  value  of  approximately  -1.0  and 
the  6-oz  cotton  poplin  a  value  of  -0.  6.  The  rear  surface  variations  in 
helium  are  compared  to  a  nearly  identical  test  in  air  (Figure  21). 

In  the  air  environment  during  the  initial  and  intermediate  stages  of 
heating  the  small  amount  of  water  present  having  a  molecular  refrac¬ 
tivity  rot  greatly  different  than  air  did  not  adversely  affect  the  inter¬ 
ference  fringe  number  or  the  determined  temperature  values  for  the 
6-oz  cotton  specimen.  As  calculated  in  Appendix  H  to;  a  1  00  percent 
saturation  of  water  vapor  in  air  at  100'C  (water  vapor  effectively  con¬ 
stituting  the  entire  gaseous  phase)  the  indicated  temperature  would  be 
too  high  by  70°C  For  this  test  with  the  6-oz  cotton  (room  humidity), 
water  vapor  did  not  exceed  a  fractional  concentration  of  it)  percent  as 
compared  with  the  complet'd/  sat.. ra;  d  cest.  With  an  hen  terminate 
percentage  concentration  of  entrapped  air  which  wouh!  not  ailed  the 
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Fringe  Shift  Variations  for  Dry  o-oz  Cotton  Poplin  in  Helium  Environment 


Figure  21.  Fringe  Shift  Variations  at  Rear  Surface  of  Dry,  6-oz 
Cotton  Poplin  in  Beth  Air  and  Helium  Environments 
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temperature  indication,  the  resultant  temperature  error  could  not 
exceed  1°C,  with  the  indicated  temperatures  more  likely  being  higher 
than  actual  temperature  by  less  than  5°C. 

Temperature  indications  deviate  more  for  heavier  cotton  and  for 
the  wool  specimens.  At  approximately  30  percent  relative  humidity 
(25  °C)  the  fiber  moisture  concentration  or  adsorption  regain  (Refer¬ 
ence  19)  for  wool  is  7.84  grams  per  100  grams  dry  fiber,  whereas  cot¬ 
ton  has  a  value  of  3.  8.  As  relative  humidity  increases  the  adsorption 
regain  values  increase  (see  Table  1).  The  ratio  of  wool  to  cotton  inter¬ 
ference  fringe  shift  values  in  helium  (lowest  negative  fringe  values 
during  intermediate  heat  stages)  suggests  an  approximation  to  the  ratio 
of  moisture  regain  values  for  the  two  fiber  materials. 


Table  i.  Adsorption  Regains  (Reference  19) 

(grn  adsorbed  water  / 1 00  gm  dry  fiber) 

-  Cotton  (0..3 19  cal/gm-  C) 
c  -  Wool  (0.325  cal/gm-  °C) 

_ Adsorption  Regains 

Relative  Wool  Cotton  Nylon 

Humidity _ (2  5*C)  (2  5  °C ) _  ( 2 4  ° C ) 
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From  these  tests  two  important  considerations  can  be  noted. 

First,  the  interferometer  temperature  indications  in  air  at  relatively 
low  humidity  levels  are  not  influenced  greatly  by  fiber  moisture  con¬ 
tent  and  may  be  considered  reasonably  accurate.  Secondly,  although 
the  measurement  of  temperature  itself  is  representative  the  effects  of 
fiber  moisture  content  on  the  temperature-time  history,  especially  for 
wool,  would  be  quite  pronounced.  Considering  the  specific  heats  of  the 
fiber  material  and  of  the  water,  and  also  the  heat  of  vaporization  for 
the  amount  of  water  in  the  fiber  material  at  a  given  relative  humidity, 
the  resultant  heat  storage  and  adsorption  of  heat  during  vaporization 
points  out  the  influence  of  water.  For  example,  the  heat  storage 
associated  with  100  grams  of  wool  {c  =  0.325  cal/gm-°C)  at  50  per¬ 
cent  relative  humidity  when  increasing  the  temperature  from  20 °C  to 
100°C  is  24  calories  for  the  fiber,  9.2  calories  for  raising  the  water 
content  to  100°C  and  61  calories  for  the  subsequent  vaporization. 

3.  3.  2.  3  Wet  Fabric  Tests 

In  the  water- saturated  fabric  heating  tests  the  fringe  shift  values 
reach  a  characteristic  plateau  at  about  lOO^C  and  remain  at  this  fringe 
value  until  the  water  has  full,  evaporated.  As  noted  earlier,  the  indi¬ 
cated  temperature  level  is  higher  by  approximately  70  degrees  for  the 
condition  of  a  saturated  (100%)  vapor  concentration  at  the  surface. 

The  exact  amount  of  water  in  the  wetted  fabric  was  not  determined  for 
each  heating  test.  In  environments  other  than  air  the  evacuation  pro¬ 
cedure  was  carried  out  with  the  saturated  specimen  housed  in  the 
thermal  test  chamber.  In  the  light  6-oz  cotton  specimens  much  of  the 
water  contained  within  the  fibers  was  outgased  during  the  evacuation 
processes,  subsequently  the  heating  tests  on  this  fabric  display  some 
fringe  shift  irregularities.  Figures  22,  21  and  2d  are  representative 
fringe  shift  variations  for  a  heated  fabric  in  the  different  environments. 
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Netted  8.  8-oz  Cotton  Sateen  in  Nitrogen  Environment 


M  §  ‘*J!4S  aSuiaj 


Figure  23.  Fringe  Shift  Variations  for  Wetted  8.  8-oz  Cotton  Sateen  in  Argon  Environment 
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Figure  24.  Fringe  Shift  Variations  for  Wetted  8.  8-oz  Cotton  Sateen  in  Helium  Environment 
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The  fringe  shift  variations  with  time  are  plotted  to  show  primarily 
the  plateau  regions  for  the  various  environments.  Since  the  moisture 
content  was  not  controlled  the  precipitation  of  more  complete  degrada¬ 
tion  after  the  water  had  been  fully  evaporated  occurred  at  different 
times  in  the  heating  phase. 


i  Since  the  molecular  refractivity  of  water  vapor  is  known,  concen¬ 

trations  of  vapor  at  the  fabric  surfaces,  and  the  surface  temperature 
can  be  calculated  from  tests  in  two  environments.  Comparisons  can 
be  made  between  combinations  of  environmental  tests  such  as  air- 
helium,  argon-helium,  and  nitrcgen-helium. 

One  such  pairing  (air-helium)  is  shown  in  Figure  25.  Water  vapor 
concentrations  cannot  accurately  be  predicted  beyond  the  vaporization 
plateau  since  the  amount  of  moisture  is  different  in  each  test  and  the 
time  function  of  fringe  shift  variations  are  not  equivalent.  The  initial 
rise  curve  may  also  vary  slightly  due  to  the  heat  storage  and  dissipa¬ 
tion  associated  with  variable  water  content. 


To  resolve  variations  introduced  by  heat  and  mass  transfer  dis¬ 
similarities,  which  are  functions  of  the  environmental  gases,  compari¬ 
sons  for  different  gas  pairs  were  made.  Representative  determinations 
of  temperature  and  concentration  at  the  rea  *  surface  for  the  8-oz  cotton 
and  8-oz  wool  are  shown  in  Figures  26  and  27.  Values  obtained  for 
helium  and  argon  compare  favorably  with  the  other  two  pairings 
(helium-air,  helium-nitrogen).  Little  difference  was  noted  here 
despite  the  fact  that  argon  (high  molecular  weight)  would  conceivably 
have  different  temperature  and  velocity  boundary  profiles  from  that  of 
the  other  two  environmental  pairings. 

One  important  finding  can  be  recognized.  At  this  level  of  heat 

2 

loading  (5  cal/cm**-s ec )  the  rear  surface  temperature  of  the  h-oz 
cotton  during  the  vaporization  plateau  period  is  effective  100°C.  This 
also  is  the  temperature  level  for  tti«'  wool  (see  Figure  27a);  (note  that 
the  concentration  of  water  vapor  at  the  rear  surface  for  wool,  Figure  27b, 
is  the  same  as  for  the  cotton)  For  lower  heat  loadings  the  rear 
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Figure  25.  Fringe  Shift  Variations  of  Water -Saturated  6-oz  Cotton  Poplin 
Rear  Surface  in  Both  Air  and  Helium  Environments 


3-37 


r  Surface 


Rear  Surface 


.JVR/tr 


u 

4, 

(a 

x 

OB 

£ 

*4 

vO 

c 

0 


i) 

« 

H 

3 

ft. 

i 

*-> 

<4 

V 

X 

v* 


v. 

js 

£ 

3 

2 

v 

£ 

v> 

£*4 


T) 

« 

c 

£ 


V 
■»-> 
4> 

Q 

OB 

c 

o 

•r* 

(« 

»-) 

v. 

<4 

> 

0) 

£ 

•1-4 

H 

i 

c 

o 

•••* 

4-* 

<4 

Vi 

c 

4) 

U 

C 

O 

u 

u 

o 

Ok 

(4 

> 

Vi 

V 
*■> 
<4 

* 

4) 

V 

(4 

V* 

v> 

3 

ifl 

iu 

(4 

4» 

06 


* J 

V 

V 

3 

ftt 


iotj«.tjuiouo3  4<>d«^  i*}»M 


3-41 


from  Identical  Tests  in  Helium  and  Argon  Environments 


surface  temperature  during  this  vaporization  equilibrium  period  is 
expected  to  be  approximately  100  °C  also.  However,  as  the  heat  flux  is 
increased  a  critical  input  energy  will  be  reached  where  higher  equili¬ 
brium  temperatures  will  be  evident  thereafter.  In  two  later  tc  ts,  one 
in  helium  and  the  other  in  COy,  at  incident  radiant  heat  loads  of 
approximately  6.  5  cal/cm  -sec,  the  temperature  values  obtained  for 
this  environmental  pa>r  were  higher  (Figure  28)  then  that  deter¬ 
mined  from  the  earlier  tests.  As  shown  in  Figure  29a,  the  fringe 
values  in  helium,  as  expected,  are  quite  different  depending  on  the 
degree  of  heat  loading.  The  concentration  values  also  are  higher 
(Figure  29(b))  for  the  higher  intensity  radiant  heat  flux. 

The  front  surface  temperature  and  concentration  values  are  ex¬ 
pected  to  be  higher  th'n  rear  surface  values.  This  was  shown  to  be 
true  as  indicated  in  Figure  30. 

A  desired  series  of  experiments  would  be  to  vary  the  heat  flux 
under  similar  type  experimental  conditions  and  characterize  vapor  con¬ 
centration  values  and  temperature  levels  for  a  series  of  representative 
fabric  materials.  For  the  wet  fabric  skin  complex  a  further  experi¬ 
mental  study  could  involve  the  important  mass  diffusion  processes  tind 
transport  mechanisms. 

3.4  Measurement  of  Fabric  Front  and  Rear  Surface  Heat  Transfer 
Coeff  ic  ients 

3.4.1  General 

Proper  evaluation  of  the  fabric  front  and  rear  heat  transfer  co¬ 
efficients  is  necessary  when  comparing  experimental  data  with 
theoretical  predictions  of  the  response  of  a  fabric -skin  system  to  a 
radiant  heat  input. 

The  heat  - 1  ransler  model  under  consideration  is  composed  of  a 
skin  layer  separated  from  the  clothing  ensemble  by  an  airspace.  The 
flow  of  heat  from  the  outer  layer  of  the  fabric  to  the  skin  is  impeded 
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by  the  insulating  effects  of  the  air  layer  and  several  layers  of  clothing. 
Heat  transfer  from  the  inner  layer  of  the  clothing  to  the  skin  is  accom¬ 
plished  by  one  or  a  combination  of  several  modes  of  heat  transfer-- 
namely,  convection,  conduction,  and  radiation.  The  radiation  environ¬ 
ment  presented  to  the  skin  is  formed  by  the  inner  layer  of  protective 
clothing.  The  amount  of  heat  transferred  by  the  mechanism  is 
governed  by  the  temperatures  of  the  surfaces,  by  the  radiation  inter¬ 
change  coefficients,  and  by  view  factors. 


The  convective  and  conductive  modes  of  heat  transfer  combine  to 
form  a  fairly  complicated  mechanism  of  heat  transfer.  If  the  air 
between  the  two  surfaces  were  still,  the  heat  would  be  transferred 
purely  by  a  conductive  process  across  the  air  gap.  The  amount  of  heat 
transferred  would  then  be  a  function  of  the  temperature  gradient  and 
thermal  conductivity  of  the  air.  This  situation  can  be  approximated  if 
the  warm  surface  or  fabric  is  above  the  skin,  and  if  no  circulation  of 
air  occurs.  If  the  fabric  is  below  the  skin,  buoyancy  forces  can  result, 
and  free  convection  will  occur.  Forced  convection  has  to  be  considered 
if  air  movement  occurs  laterally  through  the  air  gap. 

The  external  environment  plays  an  important  role  in  fixing  the 
amount  of  heat  lost  from  the  outer  surfaces  of  the  fabric  by  the  com¬ 
bined  transfer  processes.  Heat  is  lost  from  the  outer  layer  of  the 
clothing  by  a  combined  forced  and  free  convection  process  and  by  radia¬ 
tion.  The  convection  process  is  strongly  influenced  by  air  temperature 
and  surface  w.  'ds. 


3.4.2  Free  Convection  Processes 

Convection  is  a  heat-transfer  process  that  occurs  because  of  the 
mixing  of  one  parcel  of  fluid  with  another.  Forced  convection  occurs 
when  the  fluid  motion  is  produced  by  mechanical  means,  whereas  free 
convection  is  the  result  of  body  forces  produced  on  fluid  particles 
because  of  differences  in  density.  Free  convection  occurs  when  heat 
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is  transferred  to  or  from  a  surface  to  the  .adjoining  fluid  layer.  Tem¬ 
perature  changes  in  the  fluid  layer  cause  coi  responding  density  changes, 
resulting  in  fluid  movement  because  of  the  body  forces  created  by  the 
acceleration  of  gravity.  Heat  is  then  conducted  through  the  layers  and 
carried  away  by  bulk  motion  or  convection.  By  the  use  of  integral 
methods,  problems  of  laminar  and  turbulent,  free  convection  on  a 
vertical,  flat  plate  have  been  solved.  For  the  laminar  case,  the  velo¬ 
city  and  temperature-boundary  layers  were  approximated  by  second- 
order  parabolic  equations.  Using  these  methods,  the  following  equa¬ 
tion  was  developed  for  laminar,  free  convection  on  a  vertical,  flat 
plate:  20 

Nu  =  0.  508  Pr1/2  (0.  952  +  Pr)'1/4  (Gr  )1/4  ,  (6) 

y  Y 

where: 

Nu  =  Nusselt  number 

y 

Pr  =  Prandtl  number 

Gr  =  Grashof  number. 

y 

Considering  air  as  an  ideal  gas  and  assuming  a  Prandtl  number  of 
0.714,  the  following  simplified  relation  results: 

Nu  =  0.  378  (Gr  |1/4  .  (7) 

F'or  turbulent,  free  convection  on  a  vertical  plate,  a  more  com¬ 
plicated  relationship  must  be  used  for  the  velocity  and  temperature 
profiles  in  the  boundary  layer.  Using  methods  similar  to  those  used 
for  forced  convection  boundary  layers,  shear  and  heat  flow  conditions 
at  the  wall  are  modified  in  the  light  of  experimental  data.  This  method 
then  provides  the  shear  stress  at  the  wall  that  can  be  converted  to  heat 

flow  by  use  of  the  Reynolds  analogy,  and  results  in  the  following 

20 

equation: 
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Nu  =  0.0295 


(Gry) 


2/5 


(Pr) 


7/5 


1  +  0.  494  (Pr) 


2/3' 


-2/5 


(8) 


Using  simplified  relationships  for  air,  the  amount  of  heat  loss  by 
free  convection  from  a  vertical,  flat  plate  can  be  represented  by  the 
following  equation: 
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cal 
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cm  sec 


(9) 


where  the  laminar  and  turbulent  heat  transfer  coefficients  are  respectively 


and 


h 

c 


0.  27 


h 

c 


0.22  (  t  ~  t  )1/3  . 
w  e 


(10) 


Laminar  flow  exists  when  the  product  of  Grashof  and  Prandtl  numbers 

5  7  7 

is  10  to  2  x  10  ,  while  turbulent  flow  exists  when  it  is  equal  to  2x10 

to  3  x  1  0* (Reference  5). 

For  a  gas,  the  Grashof  number  is  determined  from 


G 

r 


*  1  *w  -  *e  >  y' 
2- 
V  t 


(H) 


where: 

g  =  gravitational  constant 
y  =  characteristic  length 


t  -  wall  temperature 

t  =  environment  temperature, 
e 

and 

the  Prandtl  number  from  Pr  ~  — 

k 

where: 

k  =  thermal  conductivity 
Cp  *  specific  heat  ac  constant  pressure 
pi  =  viscosity. 

A  more  complicated  situation  exists  when  die  free-convection 

heat-transfer  process  occurs  in  a  fluid  enclosed  between  two  plane 

walls.  For  horizontal  walls  no  free  convection  occurs  when  the  upper 

wall  is  heated  and  the  lower  cooled.  Here  the  buoyancy  forces  tend  to 

stabilize  the  situation.  Heat  will  be  transferred  by  conduction  only, 

and  the  temperature  profile  will  be  linear.  When  the  lower  plate  is 

heated,  however,  an  extremely  complicated  flow  situation  occurs. 

Heat  now  flows  from  the  lower  hot  plate  to  the  upper  cooler  plate,  and 

a  situation  exists  in  which  colder  fluid  particles  are  on  top  of  warmer 

ones.  An  unstable  situation  thus  exists,  and  convection  currents  are 

established  when  the  product  of  Grashof  and  Prandtl  number  reaches 

ZQ 

1700.  For  suv  li  a  situation,  Siedcntopf  made  the  flow  visible  by  tiny 
aluminum  particles.  He  found  that  a  peculiar  cellular  structure  occurs 
with  air  moving  upward  in  the  interior  of  the  cells  and  downward  along 
the  outer  portions.  This  flow  situation  oc«  uvs  until  the  Grashof- Prandtl 
number  product  reaches  around  47.000.  Above  this  value,  irregular 
turbulence  occurs. 
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It  can  be  demonstrated  by  dimensional  analysis  that  for  an  air 
layer  enclosed  between  two  vertical  walls,  the  following  relationship 
holds: 


Nu  =  f(Gr  ,  Pr,  H/L) 

where  H/L  is  the  ratio  of  air  layer  height  to  plate  spacing. 

The  above  relationship  has  been  verified  by  early  experiments. 

As  reported  by  Carlson,  *  observations  were  made  by  different  investi¬ 
gators  of  the  heat  transfer  across  plane  vertical  air  layers,  and  as 
early  as  1909  Nusselt  concluded  that  the  transfer  depends  on  the  dis¬ 
tance  between  enclosed  surfaces  and  that  for  wider  spacings,  convec¬ 
tion  effects  were  dominant  over  that  of  conduction.  It  was  shcwr.  by 
others  later  that  such  effects  as  temperature  difference  between 
bounding  surfaces,  width  to  spacing  ratio,  and  height  to  spacing  ratio 
must  be  considered.  Subsequent,  empirical  expressions  for  Nusselt 
number  as  a  function  of  these  variables  were  developed.  These  early 
studies  were  concerned  with  the  measurement  of  the  overall  heat 
transfer  occurring  across  air  enclosures  where  a  steady-state  equili¬ 
brium  temperature  existed  or  where  enclosed  surface  temperatures 
were  maintained  at  different  but  constant  values,  Carlson's  work  was 
an  extension  of  the  early  investigations  including,  in  part,  an  examina¬ 
tion  of  local  heat  transfer  phenomena. 

21 

Inapaperby  Eckert  and  Carlson,  which  was  an  extension  o* 
Carlson's  thesis,*  local  heat  transfer  coefficients  for  a  vertical-plate 
flow  situation  were  derived  from  the  temperature  gradients  in  the  air 
normal  to  the  plate  surfaces.  For  values  of  Crashed’  number  below  a 
certain  value,  a  linear  temperature  gradient  existed  oc tween  the  two 
walls  and  the  heat  was  transferred  only  by  conduction.  A s  the  Grashof 
number  was  increased,  distinct  boundary  layers  were  seen  to  form  on 
the  two  walls.  The  air  was  seen  to  flow  in  the  upper  direction  on  tne 
heated  wall  and  in  a  downward  direction  or,  the  cooler  wail.  In  the 
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central  core  between  the  two  wall  coundary  layers,  the  temperature 
was  found  to  be  essentially  constant  in  the  direction  normal  to  the  walls, 
although  it  gradually  increased  in  the  vertical  direction. 

Various  methods  have  been  used  to  calculate  the  net  flow  of  heat 

between  a  cold  and  hot  surface  in  such  a  vertical  arrangement. 

Carlson^  showed  that  the  boundary-layer  thicknesses  and  heat-transfer 

coefficients  for  this  situation  were  within  20  percent  of  those  calculated 

for  individual  heated  and  cooled  walls  assuming  as  a  free- stream  tem- 

20 

perature  the  temperature  m  the  c  ore.  Kraussold  defined  an  equiva¬ 
lent  heat  conductivity  k^,  which  when  combined  with  the  distance  between 
the  plate  anc  the  temperature  difference  allowed  a  fairly  accurate  calcu¬ 
lation  of  the  heat  flow.  A  log- log  plot  of  equivalent  thermal  conductivity 
to  molecular  conductivity  ratio  as  a  function  of  the  product  of  Grashof 
and  Prandtl  numbers  showed  that  the  conductivity  ratio  remained  at  0 
up  to  an  abscissa  value  of  3  and  then  steadily  increased  as  the  product 
was  increased.  This  same  relationship  can  be  used  for  two  concentric 
cylinders  as  well.  This  increase  in  effective  conductivity  indicates  the 
contribution  of  free  convection  to  the  heat  flow  process. 

Flow  situations  m  which  both  free  and  forced  convection  exist  can 
be  extremely  difficult  to  calculate.  In  addition  to  the  Grashof  number 
and  other  parameters  that  desi  nbe  tee  boundaries  and  orientation  of 
the  flow  system,  a  Reynolds  number  must  also  be  considered.  For 
large  Reynolds  numbers  and  smal'  Gtashoi  numbers  tin-  forced  convec¬ 
tion  predominates,  whereas  with  targe  Grashol  numbers  and  small 
Reynolds  numbers  the  free  v  «nvt“c  t  *>n  becomes  dominant.  Certain 
flow  situations  exis* ,  however  m  who  h  nether  free  convection  nor 
forced  convection  predominates.  K*.  kei  t  ^  showed  that  this  mixed -flow 
region  (d  fined  as  a  region  *n  who  h  the  heat  transfer  differs  by  more 
than  10  percent  from  either  tree  >»r  lor*  rti-  tonv  -ct: on  relations)  in 
quite  small.  A  rule  proposed  by  McAdams  to  calculate  the  heat  flow 
by  both  methods  and  use  the  one  that  provides  the  larger  heat  (low.  ^ 
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Bodoia  and  Osterle  studied  free  convection  between  two  vertical 
plates  for  conditions  when  the  flow  was  not  fully  developed.  Such  a 
situation  exists  for  short  plates  and  near  the  leading  edges  of  longer 
plates.  The  equations  of  motion  were  converted  to  a  finite  difference 
form  and  programmed  for  solution  on  a  digital  computer.  As  com¬ 
mented  on  by  E.  M.  Sparrow,  the  flow  near  the  entrance  to  the  plate  is 
essentially  one  of  forced  convection  on  which  is  superimposed  a  free- 
convection  process.  This  situation  was  studied  in  more  detail  by  Pau- 

Cheng  Lu,  who  investigated  the  case  of  a  cylinder  with  its  axis  parallel 

23  24, 

to  the  direction  of  body  forces  of  the  fluid.  In  two  papers,  Ostrach 
also  obtained  detailed  solutions  for  parallel-  plate  channels. 

3.4.3  Non-Steady  Heat  Transfer 

When  condition:',  are  changing  with  time,  some  concern  has  been 

expressed  about  whether  the  steady -state  equations  are  applicable. 

Sparrow  and  Gregg  investigated  the  case  of  free  convection  where  sur- 

2.6 

face  temperature  varied  with  time.  Under  these  conditions,  heat- 
transfer  calculations  have  often  been  made  under  the  assumption  that 
quasi-steady  conditions  prevailed.  That  is,  it  is  assumed  that  at  each 
moment  in  time  an  instantaneous  steady- state  condition  exists  and  as  a 
result  a  steady-state  heat  transfer  relationship  can  be  used.  For  a 
desired  accuracy  of  5  percent,  they  proposed  as  a  criterion  a  parameter 
that  must  be  less  than  a  certain  desired  value.  This  parameter  is  a 
function  of  temperature  change  with  time,  temperature  potential,  dis¬ 
tance,  and  free- stream  temperature.  For  many  heat-flow  situations 
this  criterion  is  not  of  concern;  but  for  such  cases  as  a  fabric -air space- 
skin  system  irradiated  by  a  time- dependent  heat  pulse,  such  as  from  a 
nuclear-energy  release,  caution  must  be  shown  m  utilizing  steady- 
state  f ree-c unvection  relationships  under  these  transient  conditions. 
Changing  ventilating  flow  rates  due  to  relative  motions  of  fabric  and 
skin  can  also  give  rise  to  non  •  stationary  effects. 
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Non- steady  phenomena  were  studied  experimentally  by  Goldstein 

and  Eckert  on  a  uniformly  heated,  vertical  flat  plate  using  a  Mach- 

2  7 

Zehnder  interferometer.  These  tests  showed  that  the  temperature 
field  in  the  fluid  developed  initially  in  the  same  way  as  heat  conduction 
into  a  semi-infinite  solid.  A  short  transition  period  then  followed,  and 
finally  a  steady- state  condition  was  reached.  Siegel  carried  out  a 
theoretical  study  of  the  same  situation.  He  used  a  method  of  charac¬ 
teristics  to  obtain  solutions  to  the  time-dependent,  free-convection 
equations  placed  in  an  integral  form  after  the  Karman-  Pohlhausen 
method.  This  analysis  predicted  that  the  boundary  layer  grows  initially 
in  thickness  with  time  and  decreases  before  it  finally  settles  on  its 
steady- state  value.  The  heat  transfer  coefficient  also  follows  this 
change;  that  is,  it  initially  decreases  with  time,  passes  to  a  minimum 
and  then  rises  again  to  its  steady-state  value.  The  results  obtained  by 
Goldstein  and  Eckert  substantiated  this  prediction.  These  studies  again 
indicate  that  caution  must  be  exercised  in  applying  steady- state  rela¬ 
tions  to  transient  free-convection  situations. 

3.4.4  Experimental  Results 

3.4.4.  1  Fabric  Surface,  Air  Gap,  and  Skin  Temperatures 

Our  experiments  involve  non- steady  or  quasi- steady  state  conditions 
for  two  vertical  facing  surfaces  using  interferometry  techniques.  To  our 
knowledge  this  configuration  has  not  been  studied  previously  with  the 
interferometer.  Idealized  correlations  of  the  heat  transfer  functions 
for  our  fabric -skin  model  cannot  easily  be  made  with  those  of  the  above 
earlier  studies  because  of  experimental  dissimilarities.  Early  investi¬ 
gators  chose  large  plate  areas  and  electrical  heat  input.  Our  specimen 
size  was  necessarily  small  to  permit  radiant  exposures  of  sufficient 
intensity.  The  height  and  width  dimensions  were  lie  Id  constant;  only 
the  space  between  the  fabric,  and  skin  was  varied.  With  a  height  dimen¬ 
sion  of  1.27  cm  the  height  to  spacing  ratios  were  varied  from  3.  7  to 
13.  5  For  a  fabric  exposure  width  of  3.18  cm  the  width  to  spacing  ratios 
thus  varied  from  9.  2  to  33.  7. 
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The  fabric-skin  configuration  is  shown  in  Figure  31.  To  simulate 
the  skin,  a  material  to  accept  heat  at  a  rate  identical  to  that  of  the  skin 
itself  was  chosen,  which  requires  matching  ol  the  respective  kc  P  and 
k/ c  P  values.  Although  skin  is  diathermanous  to  thermal  radiation  the 
simulant  need  rot  be,  provided  the  fabric  is  opaque.  For  these  experi¬ 
ments  the  8.  8-oz  cotton  sateen,  which  has  a  reflectance  of  18.  8  percent, 
transmittance  of  1.5  percent,  and  absorptance  of  79.7  percent  (for 
carbon-arc  source)  was  used.  The  amount  of  energy  transmitted  was 
considered  small  to  enable  the  choice  of  an  opaque  skin  simulant  if  de¬ 
sired. 

-4  2  4 

For  human  skin  kco  =8.48x10  cal  /cm  -sec  and  k/cp  = 

-3  2 

1.  325  x  10  cm  /sec.  Pyrex  and  fused  silica  have  thermal  properties 
similar  to  that  of  skin. 

Two  skin  simulant  models,  each  having  the  same  exposure  area  as 
that  of  the  fabric,  were  made  from  quartz.  The  thickness  of  the  quartz 
skin  simulants  were  0.  34  cm  and  0.  70  cm.  Experiments  were  made 
with  both  simulant  models  to  determine  if  a  semi-infinite  skin  assump¬ 
tion  was  valid  for  the  heating  periods  studied.  Heat  flux  was  measured 
2 

at  6  ±0,  5  cal/cm  -sec.  Figure  32  shows  interferograms  for  the  fabric- 
skin  system  at  various  time  intervals  and  spacings.  Figure  33  shows 
temperature -time  variations  for  the  fabric  front  and  back  surfaces, 
center  of  airgap,  and  the  surface  of  the  skin  simulant.  Measurements 
were  made  at  different  vertical  positions  in  the  fabric-skin  system  but 
for  the  present  discussion  we  will  present  now  temperature  data  at  the 
center  position  of  the  exposed  area  (y  =  0.  5H), 

For  these  tests  with  z  blackened  radiation  absorbing  skin,  the 
effects  of  fabric  diathermancy  can  be  seen,  the  initial  portions  of  the 
temperature  rise  curve  at  the  skin  surface  have  higher  values  than  that 
of  the  fabric  rear  suriace.  The  temperature  of  the  skin  does  not  in¬ 
crease  for  the  short  interval  after  heat  pulse  initiation  to  0.  1  second 
(conduction  into  airspace  not  felt  at  skin)  but  gradually  increases  there¬ 
after  (irregularly  for  some  tests)  to  approximately  20  to  40°C  above 
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a.  Photograph  Showing  Flux  Redistributor,  Sample  Holder, 
and  Fabric-Skin  Model  Assembly 


b.  Schematic  Drawing 


Figure  31.  Fabric -Skin  Configuration 
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t  =  0.  6  sec.  x  -0.8  sec. 

Figure  32a.  Sequence  of  Interferograms  Showing  Fringe  Shift  Variations 
for  Different  Fabric-Skin  Spacings  (i  =6  cal/cm^-sec) 

L  =  0.  094  cm  ° 
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x  =  0.  2  sec. 


X  =  0.  4  sec. 


X  -  0.  6  sec.  x  =0.8  sec. 


Figure  32d.  Sequence  of  Interferograms  Showing  Fringe  Shift  Variations 
for  Different  Fabric -Skin  Spacings  (i  =  6  cal/ cra^ sec) 

L  =  0.  245  cm  ' 
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ambient,  at  0.8  second--this  temperature  increase  depending  upon 
spacing.  The  unblackened  skin  surface  does  not  demonstrate  anin- 
itially  higher  temperature  level  as  does  the  camphor  black  coated 
skin  but  essentially  responds  to  conduction  and  convection  heat  trans¬ 
fer  in  like  manner  (refer  to  Figure  34). 

The  temperature  histories  of  the  fabric  surfaces  and  air  gap  are 
smooth  curves.  The  skin  surface  did  not  always  display  a  continuous 
temperature  increase;  frequently  a  varying  rate  was  observed  as  shown 
in  Figure  33. 

The  influence  of  fabric-skin  spacing  on  temperature  level  is  shown 
in  Figure  35  with  time  as  a  parameter.  The  effect  of  a  decrease  in 
spacing  is  to  increase  the  skin  temperature.  The  fabric  rear  surface 
temperature  at  any  instant  generally  was  higher  for  increased  spacing. 
The  lower  temperature  values  at  closer  spacings  show  further  the  in¬ 
fluence  of  an  adjacent  skin. 

The  effects  of  fabric  moisture  content  on  temperature  response 
were  studied.  In  Figures  >6a  and  36b  the  temperature  response  of  fabric 
and  skin  are  displayed  for  a  dry  fabric  specimen  and  for  one  at  room 
humidity  conditions  (relative  humidity  28  percent).  The  dry  specimen 
had  a  finite  but  undetermined  moisture  content  since  final  preparations 
for  testing  were  conducted  at  ambient  conditions.  Although  the  moisture 
content  of  cotton  at  28  percent  relative  humidity  (3.  8  percent  by  weight) 
is  not  appreciable,  the  test  results  show  tha’  the  temperature  response 
for  a  heated  room  humidity  specimen  lags  thd»  for  a  heated  dry  fabric. 
Temperature  measurements  on  skin  simulators  reported  in  Reference  6 
indicated  a  strong  influence  of  fabric  moisture  content  on  skin  response. 
To  investigate  these  effects  further,  a  more  rigorously  controlled  series 
of  tests  were  conducted.  Specimens  w(  re  housed  for  lengthy  periods  in 
a  dry  air  environment  (dew  poi  at  -60JC).  Upon  heating,  these  dry  air 
specimens  did  not  always  display  higher  temperature  response  behavior 
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34.  Temperature -Time  Variations  of  Heated  Fabric-Skin  System,  Unblackened 
Skin  Surface 


over  that  of  specimens  at  room  ambient  conditions  (sec  Figure  37). 

In  fact,  some  of  the  room-humidity  controlled  specimens  indicated 
temperature  values  somewhat  higher  than  that  of  the  dry  specimens. 

It  must  be  emphasized  that  the  temperature  values  for  the  room- 
humidity  control  tests  are  not  corrected  for  the  effects  of  water  vapor. 
Our  earlier  multi-environmental  tr-st  determinations  demonstrate 
that  due  to  the  presence  of  water  vapor  (at  approximately  28  percent 
relative  humidity),  the  temperature  indications  at  a  considered  in¬ 
stance  (at  0.  3-0.  4  sec)  are  high  by  as  much  as  7®C.  Application  of 
such  a  correction  for  the  entire  profile  would  probably  show  that  little 
statistical  difference  would  be  evident  be  tween  the  room  humidity  and 
the  dry  specimen  tests. 

Final  conclusions  as  to  the  effect  of  fabric  moisture  content  on 
the  temperature  response  characteristics  cannot  be  made  at  this  time. 
However,  a  tentative  conclusion  is  that  the  gross  behavior  of  both  dry 
and  room  humidity  prepared  sp<  cimens  appear  to  be  much  the  same, 
and  even  when  compared  with  the  wetted  fabric  tests,  the  differences 
in  temperature  response  during  the  initial  phases  of  heating  are  not 
as  large  as  one  might  expect. 

The  plot  of  data  points  in  Figure  37,  although  not  demonstrating 
any  measurable  differences  between  dry  and  room-humidity  specimens 
during  the  period  of  temperature  increase,  does  show  that  combustion 
or  more  complete  degradation  occurs  earlier  for  the  dry  samples  than 
for  the  room -humidity  test  samples  by  approximately  0.2  sec.  This  is 
a  significant  point,  although  eas.ly  overlooked  in  the  log-log  profiles, 
since  the  total  energy  delivered  prior  to  canbustion  may  be  different, 
in  the  different  tests,  by  20  percent  or  more.  The  effects  of  water 
content  thus  may  be  only  evideni  i u  this  maimer. 

From  these  conside rations  it  becomes  necessary  to  perform  ex¬ 
perimental  tests  over  a  wide  range  of  room-humidity  conditions  to 
better  assess  the  influence  of  moisture  content.  In  addition, 
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Figure  3?.  Temperature -Time  Variations  of  Heated  8.  8-oz  Cotton  Specimens  at 
Dry  and  at  Room  Humidity  Conditions 


the  multi -wavelength  or  multi -environmental  interferometric  techniques 
must  be  utilized  to  determine  water  vapor  concentrations  at  the  fabric 
surfaces  and  thus  accurately  correct  indicated  temperature  profiles. 

For  the  multi-environmental  tests,  consideration  must  be  given  the 
water  vapor  diffusion  processes  in  the  different  environments.  A  vapor 
concentration  level  at  the  fabric  surface  may  be  quite  different  in  the 
separate  environments. 

The  further  important  processes  involving  transfer  of  vapor  or 
gaseous  products  frcrri  the  heated  fabric  to  the  skin  need  to  be  understood. 
Although  the  influence  of  fabric  moisture  content  on  the  temperature 
response  of  an  underlying  or  adjacent  skin  (considering  only  heat  transfer) 
are  important,  the  effects  of  vapor  cr  mass  transfer  may  be  equally  im¬ 
portant,  and  for  prescribed  conditions  even  overshadow  the  effects  of 
pure  conduction  heat  transfer.  More  comprehensive  tests  thus  are 
needed  to  isolate  these  influencing  factors. 

In  Figure  38  the  temperature  differences  from  fabric  front  to  rear 
surfaces  for  two  fabric-airspace -skin  conditions  are  plotted  to  demon¬ 
strate  the  constancy  of  this  difference  beyond  0.  2  second.  The  im¬ 
portance  of  this  value  and  of  the  actual  temperature  recordings  would 
be  recognized  if  additional  information  on  fabric  diathermancy  (ab¬ 
sorption  in  depth)  were  available.  Since  front  and  rear  surface  tem¬ 
perature  values  can  be  obtained  from  interferometric  observations  it 
would  appear  that  determinations  of  fabric  diathermancy  would  also  be 
possible.  Using  a  modified  flash  discharge  technique  and  the  interferometer 
this  could  be  accomplished.  The  apparent  thermal  diffusivity,  thermal 
conductivity  and  heat  capacity  (important  properties  of  the  material  in 
the  absorption  and  dissipation  of  heat)  can  also  be  determined  for  non- 
homogeneous  materials.  These  property  measurements  were  not  made 
during  this  program  period  although  art  suggested  for  future  investiga¬ 
tions. 
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The  above  considerations  have  been  applied  to  the  center  position 
of  the  vertical  fabric -airspace -skin  complex.  At  different  vertical 
positions,  the  surface  temperature  histories  and  thermal  boundary 
profiles  would  also  be  different.  As  shown  in  Figure  39,  these  dif¬ 
ferences  are  markedly  evident;  the  fabric  front  surface  temperature 
varied,  for  example,  by  as  much  as  100°C  for  an  y/H  =  .  14  to 
y/H  =.  60  at  0.  6  second.  In  Goldstein's  nonsteady  state  tests  with  a 
brass  plate,  considerable  variation  in  temperature  along  the  vertical 
plane  (approximately  10  cm  height)  was  also  noted.  For  large  verti¬ 
cal  surface  areas,  the  local  heat  losses  would  decrease  with  height 
up  to  a  certain  position,  characterized  by  the  flow  conditions  and  sur¬ 
face  temperature;  thereafter  increase  slightly  with  height  for  some 
distance  and  finally  become  constant.  It  has  been  postulated  that  two 

heat  transfer  regimes  exist:  pure  conduction  and  laminar  and  turbulent 

1112 

convection  flow.  Both  Goldstein  and  Carlson,  ’  among  others, 
have  attempted  to  characterize  these  heat  flow  regimes  in  their 
studies. 

The  degree  of  temperature  difference  with  vertical  dimension 
suggests  that  considerable  importance  must  be  attached  to  heat  losses 
by  conduction  and  convection  even  for  short  periods  of  exposure  as 
considered  in  our  tests.  The  combustion  or  more  complete  degrada¬ 
tion  occurrences  were  undoubtedly  initiated  near  the  top  of  the  fabric 
specimen;  since  here  the  heat  loss  was  minimum,  resulting  in  a  locally 
higher  temperature. 


Figure  39  (continued) 


3.4.4.  2  Calculation  of  Heat  Transfer  Coefficients  from  Measured 


3. 4.  4.  2,  1  General.  -  The  local  convective  heat  transfer  from  the 

v  r  -  -  r 

fabric  may  be  determined  from  the  boundary  layer  temperature  dis¬ 
tribution  at  any  instant.  The  most  important  value  required  in  deter¬ 
mining  the  heat  loss  is  the  film  heat  transfer  coefficient,  h,  defined  in 
the  terms  of  the  heat  flow  q  per  unit  area  and  the  temperature  differ¬ 
ence  At. 


q  =  h  At 
nc  c 


Since  the  air  layer  at  the  fabric  surface  is  at  rest,  the  heat  flow  may 
also  be  calculated  from  the  conduction  equation 


'  k 


where  is  the  thermal  conductivity  of  the  air,  generally  considered 
at  the  wall  temperature,  and  the  dtMxJ  the  temperature  gradient  at 
the  surface.  By  equating  these  expressions  and  setting  the  terms  in 
form  of  the  non-dimensional  heat  transfer  coefficient,  Nusseit  number. 


we  obtain 


y  d  t 
At  dx 


where  y  is  a  characteristic  dimension. 


3.4.4.  2.2  Comparison  with  Existing,  Data.  -  Figure  40  shows  a 
plot  of  h  and  Nusseit  number  variation  with  time  for  the  front  surface  of 
a  heated  8.  8-oz  cotton  fabric.  The  characteristic  lengthy  in  the 
Nusseit  number  is  approximately  0.  5  H  at  which  horizontal  plane  the 
measurements  were  made.  As  shown,  the  coefficients  decrease  with 
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Figure  40,  Variation  of  Surface  Heat  Transfer  Coefficients  with 
Time  after  Heat-Pulse  Initiation 
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time  during  the  initial  pure  conduction  phas  .  Thereafter  a  quasi 
steady-state  condition  is  reached  where  h  is  expected  to  increase  only 
as  At  increases,  then  finally  becoming  constant  where  a  true  steady- 
state  exists.  A  true  steady- state  was  not  obtained  in  our  tests  since 
the  degradation  processes  occurred  prior  to  equilibrium  heat  in-heat 
loss  conditions.  At  the  minimum  value  of  h,  or  somewhat  thereafter, 
the  thermal  boundary  layer  was  found  to  have  a  maximum  thickness, 
decreasing  slightly  from  thereon  to  the  ignition  temperature. 

The  local  non-dimensional  representation  of  heat  flow  (Nu)  occur¬ 
ring  at  the  fabric  front  surface  in  terms  of  the  local  non-dimensional 
Gr^  number  is  plotted  in  Figure  41.  Calculations  were  made  at  five 
vertical  positions  for  times  0.  1  to  0.  8  second.  Theoretical  and 
empirical  relationships  for  quasi- steady  convection  heat  flow  predict 
for  a  fluid  of  given  Prandtl  number  the  function 

Nu  --  K(Gr)1/4  (15) 


where  K  is  a  constant. 

This  relationship  is  not  entirely  evident  from  the  plotted  data 

points.  If  we  chose  to  ignore  certain  points  (for  example,  those  at  low 

values  of  Gr  where  theoretical  and  experimental  agreements  have  been 

y 

suggested  as  being  poor)  we  may  obtain  exponent  values  from  1/8  to 

1/3.  These  data  may  be  plotted  as  shown  in  Figure  42  to  examine  the 

1  /  4 

variation  of  Nu  / (Gr  )  from  a  constant  value. 

y  y 

In  providing  an  important  non-dimensional  relationship,  such  as 

this,  consideration  must  be  given  the  most  suitable  property  values  to 

3 

use  in  the  calculations.  In  calculating  Gr  -  (gpAt  y  )/v,  the  expansion 

coefficient  of  a  gas  P  was  determined  directly  from  1/t  ;  the  kinematic 

0 

viscosity  v  was  evaluated  at  t  ,  whereas  the  wall  temperature  value 

1 2 

was  used  for  the  evaluation  of  k.  Goldstein  suggests  that  for  a  non- 
isothermal  wall,  as  in  these  studies,  l/te  for  the  expansion  coeffi¬ 
cient  is  appropriate;  however,  an  evaluation  of  y  at  wall  temperature 
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Figure  41.  Variation  of  Front  Surface  Nusselt  Number  with  Grashof  Number 


Run  No.  Ap-31 


Figure  42.  Variation  of  Front  Surface  Local  Nondimensional 

Heat  Transfer  Parameter  with  Time  (8.  8-oz  Cotton) 


calculated  for  an  assumed 


values  would  be  more  correct.  Sparrow 
ideal  gas  that  in  evaluating  Nu  number  as  a  function  ol'  Gr  the  refer- 

°  °  y  y 

ence  temperature  t  =  t  -  0.  38  (t  -  tg)  should  be  used  in  evaluating 
v,  and  k.  If  the  above  representation  required  further  accuracy, 
comprehensive  data  could  be  compiled  and  a  statistical  treatment 
followed. 


At  the  rear  surface,  for  small  spacing s  between  fabric  and  skin, 
the  presence  of  the  skin  surface  is  felt  very  rapidly.  Conduction 
occurs  across  the  air  gap  and  fully  developed  convective  flow  does  not 
occur.  The  skin  surface  temperature  is  higher  for  smaller  spacing. 

A  minimum  heat  loss  from  the  fabric  occurs  at  some  intermediate 
spacing  where  fully  developed  convective  flow  cannot  occur  and  where 
conduction  is  not  entirely  predominant.  For  our  fabric  skin  spacing 
0.  09  to  0.35  cm  the  transfer  of  heat  from  the  fabric  rear  surface  was 
by  conduction,  at  least  along  the  centerline  measurements.  At  the 
close  spacing s,  top  and  bottom  end  effects  were  also  negligible;  con¬ 
duction  heat  transfer  occurring  entirely. 


For  a  fabric  spacing  of  0.  181  cm  the  local  fabric  rear  surface 
non-dimensional  heat  transfer  coefficient  Nu  has  been  calculated  for 

y 

five  different  vertical  positions  and  different  time  instances  as  a  func- 

1/3 

tion  of  Gr^.  The  predicted  relationship  Nu^/(Gr^)  approaches  a 

constant  at  0.  6  to  0.  8  second  (steady-state)  as  shown  in  Figure  43.  At 

0.2  second  the  exponent  m  =  1/3  already  holds  precisely  since  Nu  (for 

,i/3  y 

pure  conduction)  is  equal  to  (Gr  /Gr  )  .  If  the  characteristic  length 

y 

in  the  non-dimensional  coefficient  were  chosen  to  be  the  fabric-skin 
spacing  dimension  L,  then  Nu^  =  1  for  all  values  of  Gr^  within  the 
pure  conduction  regime.  If  the  spacing  had  been  greater  (values  of 
GrT  greater  than  10^)  then  Nu  =K(GrT)*^. 


Grashof  Numbe 


3.  4.  4.  2.  3  Measurement  of  Non-Dimensional  Heat  Transfer 

Coefficients  jor  Utilization  in  'theoretical  Development 

As  presented  in  Appendix  A,  the  fabric  front  and  rear  surface  heat 
transfer  coefficients  were  non-dimensionalized  as  follows: 


and 


H. 


haeff  Lf 


10 


(16) 


(17) 


It  is  to  be  noted  that  in  the  above  relationships  the  thermal  conductivity 
is  that  of  the  fabric  rather  than  air.  In  addition,  the  characteristic 
length  is  the  fabric  thickness,  L^..  In  this  form  (Biot  number),  the 
process  of  heat  transfer  from  ait  fo  fabric  is  more  intimately  related 
to  conduction  heat  transfer  in  the  fabric  because  of  the  conductivity. 

The  variations  of  these  coefficients  with  time  (and  for  the  rear  surface, 
the  variation  with  air  gap  spacing  also)  may  be  determined.  Figure  44 
is  a  plot  of  the  non-dimensional  front  surface  heat  transfer  coefficient  I 
(Biot  number)  for  three  separate  heating  tests  on  the  8.  8-oz  cotton 
sateen.  The  film  heat  transfer  coefficient  h^  was  determined  from 
the  expression 


lfo 


Vo  “  V  i 


dt 

Sx 


fo 


where  k  was  chosen  at  t  • 

e 

4 

The  effective  fabric  value  for  k^/L^  was  taken  as  a  constant. 

(No  variation  with  fabric  temperature  was  assumed.  )  The  plot  follows 
the  predicted  variation  for  the  different  heat  transfer  regimes. 


fo 
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Dimensionless  Front  Surface  Heat  Transfer  Coefficient 


For  values  of  k,  in  the  above  calculations,  taken  at  t£Q  instead  of 
at  ,  the  more  representative  but  similar  variation  of  h^,  Figure  45, 
would  result. 

In  Figure  46  the  rear  surface  values  of  H.q  are  plotted  to  show 
variations  with  time  and  with  spacing.  Since  a  conduction  flow  regime 
exists  for  the  small  fabric  skin  spacings,  the  non-dimensional  co¬ 
efficient  is  calculated  as 


H. 

10 


k. 


\  f 


(18) 


It  should  be  noted  that  for  high  fabric  rear  surface  temperatures 
the  heat  conductivity  change  is  sufficient  to  prevent  the  formation  of  a 
linear  temperature  gradient  for  a  truly  one-dimensional  heat  flow  situa¬ 
tion.  As  the  deviations  from  a  linear  distribution  are  not  as  large  as 
would  be  expected  for  the  higher  temperatures,  it  is  possible  that 
localized  areas  of  free  convection  may  be  occurring. 


3 .  5  Comparison  of  Measured  Temperature-Time  Variation  with 
Theory 


In  Appendix  A  four  sample  calculations  are  presented  for  the 
following  conditions: 


1)  Case  I.  Opaque  Fabric- Airspace-Skin  System 

1)  Case  II.  Diathermaneous  Fabric -Airspace-Skin 
System 

3)  Case  III.  Flat  Plate  Radiating  to  a  Deep  Space 

Environment 

4)  Case  IV,  Flat  Plate  Irradiated  by  a  Time-Dependent 

Heat  Pulse. 


These  calculations  were  performed  to  allow  a  comparison  with  results 
of  Reference  4  (Case  I  and  Case  II),  Reference  29  (Case  III),  and 


i-H  ( 


Dimensionless  Air  Gap  Heat  Transfer  Coefficient 


References  30  and  31  (Case  IV).  A  limited  parameter  study  (Case  V) 
is  also  presented  to  show  the  effects  of  variation  of  the  various  heat 
transfer  coefficients. 

To  our  knowledge,  no  experimental  verification  has  existed  pre¬ 
viously  pertaining  to  the  fabric  response.  For  instance,  in  the  study 
reported  in  Reference  4,  complete  attention  was  shown  to  the  skin 
simulant  response  due  to  lack  of  an  accurate  method  of  measuring 
fabric  surface  temperatures.  The  interferometric  techniques  described 
herein  have  allowed  measurement  of  the  fabric  response.  For  this 
reason,  major  emphasis  is  shown  to  the  fabric  behavior  and  a  skin 
simulant  is  only  needed  to  properly  fix  the  rear  surface  heat  losses. 

Presented  in  Figure  47  are  experimental  results  for  fabric -skin 
simulant  responses  at  ambient  and  dry  moisture  conditions.  Particu¬ 
lar  care  was  exercised  in  conditioning  the  fabrics  for  these  tests. 

Shown  also  for  comparison  are  theoretical  predictions  for  our  particu¬ 
lar  test  conditions  (refer  to  Appendix  A  for  details  of  the  theoretical 
calculations). 

Results  for  both  test  conditions  fail  below  the  theoretical  predic¬ 
tions  for  the  fabric  response  with  reasonable  agreement  between 
theory  and  ex’-^riment.  A  similar  agreement  exists  for  the  skin 
simulant  response,  with  the  experimental  data,  being  above  the  the¬ 
oretical  curve.  The  experimental  data  tend  to  level  off  as  the  igni¬ 
tion  point  is  approached,  while  the  theoretical  predictions,  which  are 
based  on  inert  behavior,  keep  on  increasing. 

These  results  indicate  that  additional  modifications  are  needed 
in  the  present  mathematical  model.  Further,  other  effects  which 
cannot  be  identified  at  this  time,  tend  to  mask  deviations  in  fabric 
response  caused  by  moisture  conten*.  Additional  experiments  at  var¬ 
ious  heating  rates,  fabric  mpisttire  content,  and  fabric-skin  simu¬ 
lant  spacings,  are  required  to  properly  describe  the  observed  be¬ 
havior. 
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Temperature 


4.  SUMMARY  AND  CONCLUSIONS 


The  general  objective  of  the  study  reported  herein  was  to  deter¬ 
mine  the  physical  properties  required  in  a  new  material  that  could  be 
used  as  a  clothing  fabric  for  protection  against  intense  thermal  radia¬ 
tion.  To  accomplish  this  objective,  a  five-year  research  program  was 
planned.  A  curtailment  of  funds,  however,  terminated  the  program 
after  three  years  of  research  efforts.  This  Final  Report  summarizes 
findings  reported  in  the  first  two  Annual  Reports,  and  includes  a  de¬ 
tailed  description  of  the  third  year's  efforts. 

A  combined  theoretical  and  experimental  approach  was  followed  to 
allow  attainment  of  the  program  objectives.  The  studies  that  have  been 
carried  out  have  isolated  many  of  the  salient  features  of  the  behavior 
of  a  fabric-airspace- skin  system  when  exposed  to  intense  thermal 
radiation.  During  this  program  the  following  tasks  were  accomplished: 


4.  1  Theoretical  Studies 

1)  An  extensive  and  critical  review  of  previous  theoreti¬ 
cal  and  experimental  studies  reported  in  the  literature 
was  undertaken.  This  study  was  continued  during  the 
second  and  third  years  of  the  program. 

2)  To  obtain  a  basic  insight  into  the  heat  flow  processes 
involved,  a  simplified  mathematical  model  was  chosen 
and  the  heat  flow  equations  were  solved  using  an  ex¬ 
plicit  finite  difference  technique.  This  procedure  was 
programmed  for  solution  on  a  Control  Data  G-15D 
digital  computer.  Considered  in  the  analysis  are  a 
time-dependent  heat  input,  fabric  diathermancy,  and 
non-linear  radiation  and  convection  boundary  conditions. 

3)  Portions  of  the  computer  program  were  examined  by 
comparing  our  results  with  those  existing  in  the  litera¬ 
ture.  Included  in  the  calculations  are  the  temperature¬ 
time  responses  of  opaque  and  diathermanous  fabric - 
skin  systems  irradiated  by  a  constant  heat  pulse,  a 
flat  plate  radiating  to  a  heat  sink  environment,  and  a 
flat  plate  irradiated  by  a  time -dependent  heat  pulse. 


Iii  addition,  calculations  were  performed  for  an  opaque 
fabric -skin  system  exposed  to  a  constant  heat  input 
assuming  various  values  of  the  heat  transfer  coeffi¬ 
cients.  Envelopes  showing  the  limits  of  these  re¬ 
sponses  for  order  of  magnitude  changes  of  the  heat 
transfer  coefficients  were  obtained. 

4)  A  critical  examination  was  made  of  the  flash  technique 
as  a  means  of  measuring  thermal  properties  of 
materials.  We  extended  this  method,  as  originally 
developed  at  NRDL,  for  application  to  diathermanous 
materials  such  as  fabrics.  In  addition,  the  theory  was 
improved  to  allow  consideration  of  radiative  and/or 
convective  front  surface  heat  losses. 


4.  2  Experimental  Studies 

1)  Interferometer  techniques  were  perfected  for  study  of 
fabric  materials  exposed  to  high  intensity  thermal 
radiation.  Combined  with  high  speed  photographic 
methods,  these  techniques  allow  both  qualitative  and 
quantitative  measurements  during  the  initial  phases  of 
heating,  as  degradation  processes  were  initiated,  and 
when  combustion  occurred. 

2)  Interferometric  temperature  measurements  at  the 
fabric  surfaces,  in  the  air  gap,  and  at  the  skin  simu¬ 
lant  surfaces,  were  obtained  for  various  fabric -skin 
spacings,  fabric  water  content,  and  fabric  types. 

Radiant  energy  pulses  from  an  arc-image  furnace  pro¬ 
vided  constant  intensity  heat  loadings  up  to  7  cal/cm^ 
sec . 

3)  The  temperature  data  were  utilized  to  calculate  fabric 
front  surface  and  air  gap  convective  heat  transfer  co¬ 
efficients.  Instantaneous  values  of  these  coefficients 
during  the  heating  phases  were  calculated,  allowing 
determination  of  the  effects  of  transient  flow  conditions. 
The  experimentally  obtained  temperature-time  data 

for  the  fabric-skin  systems  were  further  utilized  in 
preliminary  comparisons  of  the  theoretically  predicted 
values. 

4)  Multi-environment  and  multi-wave  length  interferometer 
techniques  were  studied  and  utilized  as  methods  of 
detecting  the  release  and  determining  the  concentra¬ 
tions  of  fabric  degradation  pi-oduM.,.  Comprehensive 
heating  tests  on  dry  and  wetted  cotton  and  wool 
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fabrics  in  separate  inert  gas  environments  of  helium, 
nitrogen,  argon  and  CO2  allowed  separation  of  tem¬ 
perature  and  water  vapor  concentration  at  fabric  front 
and  rear  surfaces.  To  conduct  these  tests  a  special 
fabric  test  chamber  was  designed  and  fabricated  which 
provided  a  controlled  environment  and  made  possible 
the  optical  observation  of  fabric  samples  while  exposed 
to  radiation  from  the  arc -image  furnace. 

5)  The  flash  discharge  technique,  originally  developed  at 
NRDL,  was  perfected  in  our  laboratories.  Measure¬ 
ments  of  thermal  dilfusivity  were  made  initially  on 
opaque  electrically-conducting  samples.  The  studies 
were  then  extended  to  include  opaque  non-conductors 
and  finally,  diathermanous  materials. 


As  a  result  of  the  theoretical  and  experimental  investigations,  the 
following  conclusions  have  been  reached: 

1)  Information  obtained  during  our  studies  indicate  that 
separation  of  the  inert  behavior  of  an  "ideal"  fabric 
from  the  true  behavior  of  an  actual  fabric  is  not 
realistic.  A  coupling  of  the  semi-ideal  response  of  a 
dry  fabric  during  the  early  phases  of  heating  with  the 
actual  performance  of  a  real  fabric  during  degradation 
and  eventual  combustion  must  be  accomplished.  Only 
then  will  a  true  understanding  of  the  behavior  of  a 
clothing  fabric  under  field  conditions  be  possible. 

2)  Studies  of  the  heat  flow  processes  in  a  fabric-airspace- 
skin  system  irradiated  by  an  intense  source,  indicated 
the  usefulness  of  the  Mach-Zehnder  interferometer  in 
providing  experimental  data  previously  unobtainable. 

We  are  able  to  gain  information  pertaining  to  t»  e  fabric 
surface  temperatures,  free  convection  flow  processes 
at  the  fabric  surfaces  and  in  the  air  gap,  degradation 
phenomena,  and  a  description  of  eventual  combustion. 

From  these  measurements  maximum  temperature  and 
critical  energy  levels  of  fabric  systems  can  be  deter¬ 
mined.  Thermal  diffusion  processes  through  multi¬ 
layer  fabrics  and  the  behavior  and  interaction  of  fabric 
types  on  heat  flow  cen  be  assessed. 

3)  The  multi-wavelength  and  multi- environment  inter¬ 
ferometric  techniques  show  promise  in  identifying 
degradation  products  and  determining  concentrations. 
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Eventual  perfection  of  techniques  such  as  these  will 
allow  measurement  of  the  mass  diffusion  and  transfer 
processes  and  provide  chemical  effect  evaluations. 

4)  Other  techniques  such  as  differential  thermal  analysis 
are  needed  to  complement  interferometric  methods  in 
providing  additional  information  pertaining  to  exother¬ 
mic  and  endothermic  reactions  which  occur  in  fabrics 
during  heating.  Important  also  would  be  a  study  of 
reaction  rates,  both  for  steady-state  conditions,  and 
as  a  function  of  temperature  changes  with  time. 

Further  experimental  efforts  must  be  expended  to  pro¬ 
vide  a  full  description  of  the  heat  flow  processes. 
Variations  in  heat  loading,  fabric  type  and  configura¬ 
tion,  conditions  of  moisture  content  and  chemical  addi¬ 
tion,  and  ventilation  effects  and  layer  -  spacing  systems 
are  a  few  of  the  desired  experimental  considerations. 

5)  The  mathematical  model  should  be  improved  to  include 
consideration  of  exothermic  and  endothermic  chemical 
effects  and  other  real  fabric  behavior.  Until  this  in¬ 
formation  is  available,  other  improvements  in  the 
model  are  not  feasible,  as  the  chemical  effects  over¬ 
power  the  "ideal"  behavior  assumed  in  present  analy¬ 
ses.  The  initiation  of  effective  heat  "sources"  and 
"sinks"  caused  by  chemical  effects  can  be  incorporated 
into  the  finite  difference  calculation  algorithm. 


Additional  theoretical  and  experimental  studies  are  required  before 
a  proper  understanding  of  the  complicated  fabric  heating  processes  is 
obtained.  An  integrated  program  combining  the  research  efforts  of 
specialists  in  physical  chemistry,  heat  transfer,  mathematical  analy¬ 
sis,  and  instrumentation  is  required  The  necessary  experimental  and 
theoretical  techniques  have  now  been  sufficiently  developed  to  allow 
reasonable  expectation  of  success  in  the  near  future. 
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I.  SELECTION  OF  COMPUTATION  METHOD 

The  heat  flow  equation  and  boundary  conditions  describing  the 
fabric -skin  heat  transfer  processes  as  presented  in  Section  2. 3  are 
as  follows: 

A.  Heat  Flow  Equations  in  Fabric  and  Skin 
1.  Fabric 


(A-l) 


(A-2) 


B.  Boundary  Conditions 


1.  At  the  Outer  Surface  of  the  Fabric 
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2.  At  the  Inner  Surface  of  the  Fabric 
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3.  At  the  Skin  Layer 
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The  system  of  equations  and  boundary  conditions  under  present 
consideration  is  not  subject  to  exact  solution  by  known  analytic  tech¬ 
niques.  If  the  problem  is  attacked  in  its  most  general  form,  the  combi 
nation  of  nonlinear  boundary  conditions  and  temperature -dependent^ 
physical  properties  is  sufficient  to  force  the  use  of  numerical  methods. 
When  an  analytical  sclution  of  the  heat-flow  equation  is  possible, 
temperature  can  be  expressed  as  a  continuous  function  of  the  time 
and  space  variables.  In  a  numerical  solution,  the  differential  repre¬ 
sentation  is  replaced  by  a  discrete  or  finite  representation.  A  method 
of  spatial  subdivision  is  chosen,  and  temperatures  are  calculated  only 
for  certain  pre-selected  points  that  are  assumed  to  be  representative 
of  the  region,  including  the  point.  In  a  transient  heat -flow  problem 
proper  subdivisions  must  also  be  made  in  the  time  variable. 

Over  the  years,  various  computation  systems  have  been  developed 
for  numerical  analysis  of  transient  Leat-flow  problems.  In  three  of 
these  methods,  the  partial  differential  equations,  and  the  boundary  and 
initial  conditions  are  replaced  totally  or  partially  by  difference  repre¬ 
sentation  of  space  and  time  variables.  In  one  method,  the  partial 
differential  equations  are  replaced  by  a  system  of  equations  with  a 
difference  representation  of  the  space  variable,  while  the  time  variable 
remains  in  differential  form.  *  The  other  two  methods  utilize  differ¬ 
ence  representations  of  both  the  space  and  time  variables.  The  first 
of  these,  the  explicit  method,  allows  calculation  of  the  temperature  at 
a  future  time  on  the  basis  of  temperatures  and  heat  inputs  or  heat 
generation  existing  during  the  present  time  increment. 

'  The  finite  difference  form  of  Equation  (A -2)  can  be  obtained  by 
performing  a  heat  balance  on  an  internal  element,  as  shown  in 
Figure  A-l.  The  net  heat  transferred  to  an  element  can  be  equated  to 
the  change  of  elemental  heat  content.  In  a  solid,  the  heat  added  to  a 
body  is  equal  to  its  increase  in  enthalpy,  which  for  a  constant  specific 
heat  can  be  expressed  as  a  product  of  specific  heat  and  temperature 
change.  For  internal  points  in  the  system,  the  primary  mechanism 


A  -3 


of  heat  transfer  to  and  from  the  element  is  by  conduction,  except  when 
the  material  is  diather manors,  and  the  temperatures  are  high  enough 
to  cause  ihterelement  radiation  interchange.  If  heat  .,  is  being  gener¬ 
ated  in  the  element  due  to  chemical  reaction  or  diathermancy,  an 
additional  term  must  be  included  in  the  calculation.  For  the  case  of 
conduction  heat  transfer  only  and  no  heat  generation, 
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where  Ax  and  AT  are  increments  of  the  space  and  time  variables  and 
the  set  of  points  in  the  x,  t  plane  are  given  by  x  *  jAx  and  t  *  nAx  . 
Rearranging  slightly, 
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Equation  (A-7)  is  an  explicit  representation  of  Equation  (A-2)  using 
forward  time  differences;  i.  e. ,  the  temperatures  at  the  beginning  of  the 
time  increment  A  Tare  used  to  calculate  the  temperature  at  the  e  id  of 
the  time  increment.  The  second,  or  implicit  method,  includes  con¬ 
sideration  of  conditions  at  the  end  of  the  time  increment.  Expressing 
Equation  (A-2)  in  an  implicit  form  utilizing  backward  time  differences 
we  obtain 
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which  is  similar  in  form  to  Equation  (A-7),  except  that  the  second 
derivative  in  space  is  taken  at  the  end  of  the  time  increment.  Applica¬ 
tion  of  this  method  requires  solution  of  a  set  of  simultaneous  equations 
at  each  step  in  time. 
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Considering  further  the  explicit  finite  difference  method  and 

solving  for  t.  in  Equation  (A- 7),  one  obtains 
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where 
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a  dimensionless  modulus.  If  a  value  of  2  is  chosen  for  the  modulus  M, 
the  temperature  of  the  element  at  the  end  of  the  time  increment  at  is 
simply  the  average  of  the  previous  temperatures  of  the  adjoining  ele¬ 
ments.  The  calculation  can  then  be  carried  out  graphically  using  the 
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familiar  Binder -Schmidt  method.  Nessi  and  Nissole  discussed  the 

4 

use  of  other  moduli,  and  Dusinberre  presented  a  graphical  con- 
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struction  for  M  *  4.  Paul  showed  how  the  method  could  be  generalized 
even  further  for  other  moduli. 

It  is  known  that  implicit  analogs  are  stabl  j  for  all  division  of 
space  and  time  while  the  semi -discrete  analog  is  stable  for  all  divi¬ 
sions  of  the  space  increment  and  allows  the  use  of  continuous  or 
analog  time  computers  for  certain  heat  flow  problems.  *  Explicit 
calculations,  however,  are  subject  to  several  restrictions  to  insure 
stability.  The  stability  condition  is  defined  as 


M 


If  stability  restrictions  are  not  observed,  large  undamped  fluctuations 
can  sometimes  appear  in  the  solution  that  make  the  calculation  worth¬ 
less.  Calculations  are  usually  made  with  a  finite  net  size;  if  the 
calculation  system  is  unstable,  the  symptoms  appear  early.  The 
stability  condition  equation  indicates  that  making  the  space  divisions 
small  without  also  decreasing  the  ume  increment  could  '■-vse  in¬ 
stability  problems. 
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At  the  boundary,  an  additional  restriction  can  be  place!  on  M, 
depending  on  the  method  of  subdivision.  If  convective  heat  transfer 
is  occurring  at  the  surface  and  a  division  like  that  in  Figure  A- lb 
is  used,  the  heat  balance  is 


(t 

k~ 


A  x 


‘  h^0  ‘  V  *  Cp 


A  x 

T37 


(A-10 


To  simplify  the  notation,  the  temperature  at  the  end  of  the  time  in¬ 
crement  is  designated  by  a  prime  and  the  unprimed  values  represent 
temperatures  at  the  beginning  of  the  increment.  Solving  for  t^',  one 
obtains 


t  ' 
o 


2Bi  t  +  2t.  +  t  (M  -  2Bi  -  21 
e  1  o 

M 


(A- 12) 


where  the  Biot  number  Bi,  is  defined  as  —  .  Equation  (A- 12)  intro¬ 
duces  the  requirement  that  M*  (2Bi  +  2).  If  this  requirement  is  not  met 

4 

(according  to  Dusinberre  ),  an  absurd  physical  situation  develops, 
namely,  the  warmer  region  0  is,  the  colder  it  will  be  after  the  time 
intervalAt . 

The  above  limitiation  can  be  circumvented  if  a  subdivision  sug- 
4 

gested  by  Eh. ainberre  and  shown  in  Figure  A-lc  is  used.  As  no  heat 
capacity  is  associated  with  point  0,  the  surface  temperature  is  cal¬ 
culated  using  a  steady  state  relationship,  or, 

(tl  •  ‘o' 

-H(to  -  te)  .  0  (A-13) 

and 

t  ♦  2t .  /  Bi 

*0  ■  -TTE7Br 

This  t  is  then  used  to  calculate  the  transient  t,  after  the  time  incrementAr  . 
o  * 
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Thought  must  be  given  not  only  to  stability  considerations  but  also 
to  other  errors  which  can  arise.  Truncation  error  arises  because  the 
difference  representation  differs  from  the  original  problem,  a  fact 
which  places  the  mathematical  burden  on  the  user,  who  must  prove 
that  the  solution  of  the  difference  equation  converges  to  that  of  the 
differential  equation  as  the  space  and  time  increments  are  refined. 
Round  off  error  appears  because  the  difference  equation  is  not  solved 
exactly,  the  computation  being  accomplished  with  numbers  of  finite 
length.  Rounding  errors  are  not  usually  of  much  consequence.  If 
stable  equations  are  used,  the  rounding  errors  are  not  amplified  with 

time  and  accumulate  roughly  as  the  square  root  of  the  number  of  cal- 

,  .  6 

culation  steps. 

Analytical  formulation  of  stability  criteria  and  truncation  error 
estimates  for  a  finite  difference  approximation  of  the  heat  flow  equation 
depend  upon  the  particular  computing  algorithm  used.  If  an  analytical 
solution  is  available  or  if  experimental  data  are  taken  in  such  a  way 
as  to  satisfy  all  requirements  of  the  original  equations,  a  good  insight 
into  these  errors  can  be  obtained.  If  sufficient  computer  time  is 
available,  the  rate  of  convergence  of  the  finite  difference  approxima¬ 
tion  to  the  "true"  solution  can  be  studied  by  refining  the  calculation 
net.  System  instability  can  easily  be  seen  by  examination  of  the 
temperature -time  curves.  The  major  danger  m  consideration  of 
the  appropriateness  of  the  finite  difference  approximation  is  to 
assume  that  refinement  of  the  calculation  mesh  is  a  sufficient  condi¬ 
tion  to  insure  accuracy  of  solution.  Even  though  the  exact  solution 
closely  approaches  the  difference  results  as  the  space  and  time  in¬ 
crements  are  refined,  it  is  possible  that  the  exact  solution  of  the 
difference  equations  can  depart  considerably  from  the  true  solution 
of  the  initial  value  problem,^ 

The  choice  of  whether  to  use  a’1  explicit  or  implicit  method  of 
calculation  depends  mamly  on  the  problem  to  be  solved  and  the  avail¬ 
ability  of  high-speed  computers.  Very  often  stability  requirements 
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imposed  on  explicit  methods  dictate  that  an  excessive  number  of  time 
and  space  increments  must  be  used  to  provide  a  meaningful  solution. 
Then  hand  calculations  become  impossible  and  even  machine  computa¬ 
tion  time  must  be  considered.  It  is  then  possible  that  an  implicit 
scheme  could  be  utilized  with  more  efficiency.  In  many  cases  the 
advantage  of  fewer  increments  in  time  and  space  because  of  less  rigid 
stability  requirements  overrides  the  disadvantage  of  more  complex 
calculations  at  each  time  increment. 

In  the  results  reported  in  Reference  7,  a  graphical  explicit  finite 
difference  technique  was  emphasized,  and  the  results  so  obtained  were 
in  good  agreement  with  the  experiments.  Our  initial  requirements 
were  also  satisfied  by  an  explicit  technique  designed  for  computer 
solution,  but  it  is  felt  that  an  implicit  technique  would  allow  easier 
handling  of  temperature -dependent  properties  by  relaxing  the  stability 
restrictions.  The  details  of  transforming  the  heat  flow  equations  into 
an  explicit  finite  difference  form  are  presented  as  follows. 


II.  EXPLICIT  FINITE  DIFFERENCE  EQUATIONS 


A.  Development  of  Equations 

To  uncover  the  important  parameters  describing  the  fabric -skin 
heat  transfer  situation,  and  to  obtain  an  insight  into  the  salient  features 
of  finte -difference  techniques,  a  simplified  explicit  finite -difference 
procedure  was  chosen  for  our  initial  studies.  The  stability  re¬ 
strictions  imposed  on  this  method  tend  to  counterbalance  its  main 
virtue  of  simplicity  by  requiring  an  excessively  large  number  of  calculation 
steps  with  small  time  increments.  The  stability  of  such  a  technique*  for 
conditions  of  variable  properties,  is  also  open  to  question.  For  the  above 
reasons,  it  is  felt  that  an  implicit  procedure  would  be  more  suitable  for  future 
studies.  However,  the  insight  obtained  into  the  relative  importance 
of  the  various  parameters  and  heat  transfer  boundary  conditions  more 
than  compensate  for  the  effort  expended* 

For  the  p  resent  problem  the  boundary  subdivision  shown  in  Figure 
A-2  is  used.  In  addition  to  convection  at  the  surface  0,  radiant  input  and 
reradiation  as  a  function  of  surface  temperature  are  considered.  In 
addition,  a  moderate  amount  of  fabric  diathermancy  is  assumed. 

The  radiation  absorbed  by  element  1  is  assumed  to  be  released  as  heat 
at  the  surface  0  and  is  included  in  the  heat  balance.  For  the  other  ele¬ 
ments  the  absorbed  radiation  is  assumed  to  be  released  at  their  centers. 


Figure  Fi'Ute  Difference  Subdivisions 
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In  the  revised  Schmidt  graphical  method,  a  half  slice  having  no 
heat  capacity  is  added  outside  the  surface.  The  temperature  at  the 
outer  surface  of  this  element  is  obtained  graphically  as  a  function  of 
surface  and  environmental  temperatures.  This  procedure  is  discussed 
in  detail  in  References  2  and  7  and  will  not  be  presented  here.  This 
concept  allows  a  simple  graphic  determination  of  the  transient  temper¬ 
atures  of  element  1,  as  a  function  of  the  fictional  element  surface 
temperature  and  element  2  temperature,  using  the  familiar  Schmidt 
representation.  The  accuracy  of  this  method  is  compromised  some¬ 
what  because  in  the  calculation  of  t^'  the  fictional  element  is  assumed 
to  possess  the  thermal  properties  of  the  real  system.  This  assump¬ 
tion  enters  implicitly  when  Equation  (A-9)  is  used  to  calculate  tj'.  In 
addition,  the  initial  surface  temperature  is  assumed  equal  to  the  en¬ 
vironmental  temperature,  which  is  not  completely  compatible  with  the 
other  assumptions  made.  A  steady  state  heat  balance  is  made  at  the 
surface,  which  has  no  heat  capacity,  and  which  should  immediately 
sense  any  heat  transfer  and  react  instantaneously. 

The  method  used  in  Reference  7  is  essentially  as  discussed  above, 
except  that  the  initial  surface  temperature  is  allowed  to  seek  a  value 
different  from  the  environmental  value,  which  tends  to  counterbalance 
the  error  arising  in  the  calculation  of  t ^ using  the  Schmidt  equation. 

With  a  modulus  of  2,  it  can  be  shown  that  the  equations  for  the 
surface  temperature  and  transient  temperatures  for  element  1  derived 
herein  are  equivalent  to  those  presented  in  Reference  1,  although  some 
of  the  assumptions  and  internal  details  of  the  developments  differ 
slightly.  The  method  presented  herein  is  not  intended  for  graphical 
solution  and  allows  the  use  of  moduli  other  then  2. 

The  heat  balance  at  the  outer  surface  becomes  (see  Figure  A-l) 
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which  is  the  finite  difference  form  of  Equation  (A-3).  The  finite  differ¬ 
ence  form  of  Equation  (A-l)  for  element  1  is 


(A- 16) 


For  element  2  and  all  others  in  the  fabric  the  diathermancy  term  is 
included.  For  all  elements  except  the  rear  element  the  general  or  jth 
term  is 


kf  +  kf  (  +  1  exp  [  -**£>] 
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The  last  or  kth  element  equation  is  of  the  form  of  Equation  (A -17)  ex¬ 
cept  that  the  first  term  is  multiplied  by  2,  as  in  Equation  (A- 16),  because 
of  the  half  increment  distance  used  in  the  subdivision: 
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The  heat  balance  at  the  inner  surface  of  the  fabric,  Equation  (A-18), 
in  finite  difference  form  is 
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where  the  radiation  interchange  term  is  calculated  on  the  basis  of  parallel 
infinite  flat  plates  radiating  as  grey  bodies  (refer  to  Table  A.  1). 
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Table  A.  1  Dimensionless  Groups 


Radiation  input  into  system  I  =  - - ,  where 


:V3«»«;- 


At  the  skin,  surface  Equation  (A-19)  becomes 
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Because  the  skin  is  assumed  to  be  opaque,  all  of  the  transmitted  radia¬ 
tion  is  absorbed  or  reflected  at  the  surface.  For  these  calculations, 
the  reflected  part  is  neglected.  The  equation  for  the  first  skin  element  is 
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The  equations  for  all  succeeding  elements  are  of  the  form 
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If  the  properties  are  allowed  to  vary  with  temperature,  the  finite 
difference  calculations  must  be  made  with  the  above  equations.  If  the 
properties  are  assumed  constant,  the  solution  can  be  made  more  general 
by  nondime nsionalizing  the  equations  with  a  group  of  dimensionless 
parameters.  Unfortunately,  the  diathermancy  term  remains  ex¬ 
plicitly  in  the  equations.  As  shown  in  Section  III,  separate  calculations 
also  have  to  be  made  for  various  bomb  yields  when  a  time -dependent 
heat  input  from  a  nuclear  energy  release  is  considered.  Introducing 
the  dimensionless  groups  presented  in  Table  A.  1,  the  equations  for 
the  fabric  and  skin  become 


1)  Fabric  outer  surface 
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2)  First  element  of  fabric 


(A-21a) 


8)  General  skin  element 


(T.  ,  -  T.)  +  (T.  ,  -  T.)  +  (T,1  -  T.)  (A-22a) 
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When  the  fabric  ie  in  good  thermal  contact  with  the  skin,  which  is  r> 

assumed  to  be  opaque,  any  radiation  transmitted  through  the  fabric  will 
be  released  as  heat  at  the  interface.  The  steady-state  heat  balance  at  the 
interface  i  is, 
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It  might  be  noted  here  chat  transmitted  short-wavelength  radiation  re¬ 
flected  from  the  skin  is  again  assumed  negligible  and  is  not  considered 
in  the  fabric  heat  balances.  The  above  equation  in  dimensionless  form  is 

2(T,  -  T.).  2(T.  -  T  ) 

*so  e*p(-fLf)  .  0  (A-23a) 

X  s 

For  this  configuration  the  above  equation  replaces  Equations  (A- 19a)  and 

(A -20a), 


B.  Solution  of  Equations 

Equations  (A-l6a),  (A-17a),  (A-18a),  (A-21a),  and  (A -22a)  are  in 
an  explicit  finite  difference  form.  That  is,  tho  temperatures  at  the  end 
of  the  time  increment  A  9  are  based  on  conditions  existing  at  the  beginning 
of  the  time  increment.  The  boundary  conditions,  as  represented  by 
Equations  (A- 15a).  (A- 19a),  and  (A-20a)  are  of  a  steady-state  form. 

This  comes  about  because  of  the  particular  method  of  boundary  sub¬ 
division  utilized.  As  shown  in  Figure  A-i,  no  heat  capacity  is 
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associated  with  surface  0.  Thus,  the  heat  balance  for  this  surface 
does  not  include  consideration  of  a  heat  storage  term,  and  the  surface 
temperature  is  a  function  of  the  heat  transfer  occurring  at  the  same 
instant  in  time. 

Equations  (A-15a),  (A-19a),  and  (A-20a)  being  nonlinear  were 

g 

solved  using  the  Newton-Raphson  method  of  determining  roots.  This 
is  an  iterative  scheme  which  converges  quite  rapidly  for  our  particular 
problem  because  of  the  smoothness  of  the  temperature -time  function. 
Equation  (A -15a)  involves  only  one  variable,  while  Equations  (A -19a) 
and  (A-20a)  must  be  solved  simultaneously  as  the  fabric  and  skin 
temperatures  both  appear  in  each  equation.  The  Newton-Raphson 
method  utilizes  a  recurrence  formula  of  the  type 


k+1  k  f(x^) 
f'(xk) 

to  obtain  a  real  solution  of  an  equation  of  the  form 


(A -24) 


f(x)  =  0 


(A-25) 


Successive  applications  of  this  iteration  result  in  a  value  of  x  which 
satisfies  Equation  (A-25).  This  method  w<ts  applied  to  obtain  a  solution 
for  Equation  (A-I5a). 

The  above  method  was  also  uo-  i  to  obtain  solutions  tor  the  two 
simultaneous  equations  represented  by  Equations  (A -19a)  and  (A-20a). 
These  equations  are  of  the  form 
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These  values  are  inserted  in  Equations  (A-26)  and  (A-27)  and  expanded 
by  Taylor's  theorem  for  a  function  of  two  variables.  Neglecting  all  non¬ 
linear  terms,  the  resulting  simultaneous  equations  can  be  solved  by  de¬ 
terminants  to  yield  the  corrections  m  and  n. 
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Successive  applications  of  Equations  (A-30)  and  (A-31)  result  in 
values  of  x  and  y  which  satisfy  Equations  (A-26)  and  (A-27). 

Applying  Equation  (A-24)  to  Equation  (A-15a)  and  Equations  (A-30) 
and  (A-31)  to  Equations  (A-i9a)  and  (A-20a)  and  solving  Equations 
(A-16a),  (A-17a),  (A-18a),  (A-21a)  and  (A-22a)  for  temperature  at  the 
end  of  the  time  increment,  the  following  system  of  equations  results: 
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1)  Fabric  outer  surface 
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2)  First  element  of  fabric 
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3)  General  fabric  element 
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4)  Last  fabric  element 
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5)  Interface  boundary  conditions 


a)  Fabric  inner  skin 
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b)  Skin  surface 
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The  subscript  notations  on  the  above  functions  represent  partial  deriva¬ 
tives  with  respect  to  the  subscript,  or 
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6)  First  element  of  skin 
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7)  General  skin  element 
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The  validity  of  the  proposed  model  and  system  of  equations  pre¬ 
sented  herein  can  be  determined  only  if  experimental  results  obtained 
under  equivalent  conditions  are  available  for  comparison.  The  un¬ 
certainties  regarding  values  of  the  heat  transfer  coefficients  and 
temperature  dependence  of  the  thermal  and  optical  properties  remain 
yet  to  be  resolved.  For  fabrics  which  have  a  low  value  of  thermal 
conductivity,  front  surface  degradation  effects  due  to  the  high  local 
temperatures  can  greatly  change  the  optical  properties.  In  our  flash 
tube  experiments  on  plastics,  it  was  found  that  these  effects  did  occur. 
These  findings  were  substantiated  by  the  theoretical  predictions  of 
Reference  9,  which  indicated  a  maximum  front  surface  temperature  of 
around  600  C. 

A  computer  program  utilizing  the  Daisy  interpretive  routines  for  the 
Control  Data  G-15  computer  was  written  to  obtain  solutions  for  the  system 
of  finite  difference  equations  represented  by  Equations  (A-33)  through 
(A-42).  The  program  is  entirely  general  and  allows  variation  of  all 
the  nondimensional  parameters  presented  in  Table  A.  I  for  any  one 
calculation.  In  addition,  the  initial  temperature  conditions  and  en¬ 
vironmental  temperature  can  vary.  However,  for  diathermanous  cal¬ 
culations,  the  particular  way  in  which  the  energy  is  assumed  to  be 
absorbed  in  depth  restricts  the  calculation  to  opaque  or  only  slightly 
diathermanous  materials.  The  modification  required  to  handle  highly  dia- 
thermanous  materials  is  minor,  but  for  the  present  this  limitation  will  be 
observed. 

A  schematic  of  the  program  is  shown  as  Figure  A-a. 
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III.  COMPARISON  OF  PRESENT  CALCULATION  PROCEDURE  WITH 
OTHER  RESULTS 

Although  experimental  data  are  limited,  it  was  possible  to  check  out 
various  portions  of  our  total  computer  program  by  comparison  with  other 
theoretical  results  available.  In  Reference  7,  computer  solutions  were 
presented  both  for  opaque  fabric  and  diathermanous  fabric -airspace -skin 
systems.  In  this  study,  the  radiation  and  convection  heat  exchange  co¬ 
efficients  were  linearized  into  a  single  term  and  a  square -wave  radiation 

heat  input  term  was  assumed.  The  fabric  was  divided  into  six  slices 

2  . 

and  the  modulus  M,  defined  as  AX  /  A  0,  was  the  same  in  both  the  cloth 
and  skin.  As  the  actual  time  increments  were  also  selected  as  being 
equal,  the  AX  ratio  of  fabric  to  skin  was  found  to  be  inversely  propor¬ 
tional  to  the  square  root  of  the  k/cp  ratio  of  skin  to  fabric.  The  same 
conditions  were  assv  led  to  perform  the  following  calculations. 

A.  Case  I.  Opaque  Fabric -Airspace -Skin  System 

1.  Initial  Conditions 

T1  =  T2  "  T3 . etc*  =  0 

T  =  0 

e 

T,  >  0 
fo 

2,  Values  of  Parameters 


F.  *  0 

J.O 


e 


0 


Hfo  «  0.04 
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0.50 


I  =  1 

AXf  =  1/6, 3  AX  =  1/30 

A  S 

AXf2 

Mf  =  Ms  =  2,  Aef  =  — *r~  -  1/72 

These  calculations  and  equivalent  points  from  the  MIT  calculations^ 
are  shown  in  Figure  A -5  for  comparison. 

B.  Case  II,  Diathermanous  Fabric -Airspace -Skin  System 

1.  Initial  Conditions 

T }  =  T2  *  Tj . etc.  =  0  (only  internal  points) 


Hfo  =  0.04 

H.  =  0.50 
xo 
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A  T 


1 


1 


A  X{  =  1/6,  A  Xg  =  1/30 

AXf2 

M £  s  M ^  s  2,  A0 ^  s  — ^ —  —  1  /  <  2 

These  calculations  and  equivalent  points  from  the  MIT  calculations 
are  shown  in  Figure  A-6  for  comparison. 

The  good  agreementbetween  the  present  results  and  the  MIT  data 
was  sufficient  to  show  that  the  calculation  procedures  are  equivalent. 
Because  of  close  correlation  between  both  sets  of  data  during  the 
earlier  phases  of  the  temperature -time  history,  it  was  not  felt  neces¬ 
sary  to  continue  the  calculation. 

A  good  insight  into  the  effects  of  the  nonlinear  radiation  term  was 
provided  by  comparison  with  the  recent  data  of  Reference  10.  In  this 
study,  a  finite  difference  technique  adapted  for  computer  solution  was 
utilized  to  present  front  and  back  surface  temperature  variation  as  a 
function  of  time  for  a  flat  plate  radiating  into  deep  space  with  an  in¬ 
sulated  back  surface.  For  this  comparison  the  following  calculation 
was  performed. 

C.  Case  III.  Flat  Plate  Radiating  to  a  Deep  Space  Environment 
1.  Initial  Conditions 


T  =  0 
e 


T 


2.  Values  of  Parameters 


0 
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Figure  A- 6.  Temperature  HisU  ry  within  a  Diathermanous - 
Fabric -Skin  System 


e 


=  0 


I  =  0 

AXf  =  1/10 

Mf  *  2,  3 

Rfo  3  0.1,  0.2,  0.4,  1.0,  2.0,  4.0,  10.0,  and  20. 0. 

This  calculation  is  presented  in  Figure  A-7,  where  is  t^/t^. 

Neither  the  modulus  nor  number  of  plate  subdivisions  AXj  were 
soecified  in  Reference  10,  so  values  of  2  and  1/10,  respectively,  were 
utially  chosen.  The  resulting  curves  indicate  that  as  fhe  radiation 
coefficient  increases,  an  instability  appears  in  the  solution  which  gets 
progressively  worse  for  each  increase  in  R^.  However,  the  instability 
is  damped  out  after  a  certain  number  of  nondimensional  time  increments. 
For  R^o  =  20,  the  modulus  was  increased  to  3,  which  improved  the 
accuracy  of  the  calculation.  Division  of  the  plate  into  a  larger  number 
of  increments  also  increases  the  stability  of  the  calculation.  This 
calculation  was  performed  for  values  of  R^q  much  larger  than  those  to 
be  encountered  in  a  fabric  heat  exchange  situation,  and  the  instability 
would  not  normally  be  a  problem.  However,  the  results  obtained  bring 
to  light  a  basic  feature  of  f:nite  difference  calculations  as  regards  in¬ 
stabilities  and  boundary  subdivisions. 

D.  Case  IV.  Flat  Plate  Irradiated  by  a  Time  Dependent  Heat  Pulse 
1.  Initial  Conditions 

T  »  T,  •  T,  *  T, . etc.  ■  0 

o  J  2  3 

T  =  0 
e 
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rature  Change  or  Outer  Surface  of  Plate  Radiating 


2.  Values  of  Parameters 


R,  -  R.  =  0 

fo  10 


■TL 


=  0 


Hf0  =  0.10 


H.  =  0 
10 


I  =  I(AQ) 
A  =  1 


AXf  =  1/6 


AX, 


Mf  =  2,  A  0f  s  — =  1/72 


A  calculation  was  performed  to  investigate  the  behavior  of  the 
explicit  finite  difference  calculation  procedure  for  the  case  of  a  time 
dependent  heat  input.  The  dimensionless  presentation  of  radiation 
intensity  variation  with  time  for  a  nuclear  energy  release*  * 
is  given  by  the  following  equations: 


Region 
0  S  AO  S  0.3 
0.  3  S  AG  S  0.65 
0.65  S  A0  S  1,4 
1.4  f  A0  i  10 


Equation 
I  =  0.  565  A0 
I  =  1.72  A0  -  0.347 
1=1-  1.42(AQ  -  1.  05)£ 
I  =  1.434(AQf 1,64 


where  A  is  the  reciprocal  of  the  Fourier  number  evaluated  at  t  ,  or, 

max 
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and 


A 


T 

max 


f 


T 

max 

Batter**  presented  plotted  data  for  a  flat  plate  with  an  insulated  back 

surface  irradiated  by  a  heat  pulse  with  a  shape  described  by  the  above 

equations.  Comparison  of  the  present  method  with  these  results  was 

made  for  a  selected  A  (proporational  to  explosion  yield)  of  1.  0  and  a 

front  surface  heat  loss  parameter,  H,  of  0.  10.  This  comparison  is 

10  12 

presented  in  Figure  A-8  along  with  data  of  Loria  and  Schmit  who  used 
a  method  very  similar  to  Batter's.  **  Presented  are  the  front  surface, 
back  surface  and  center  temperatures  of  the  plate. 

An  initial  parameter  study  was  performed  to  obtain  a  better 
understanding  of  the  relative  significances  of  the  various  parameters 
isolated  earlier  in  this  appendix. 

To  reduce  the  scope  of  the  study,  only  a  constant-time  heat  input 
and  an  opaque  fabric  were  considered.  To  isolate  and  study  the  radia¬ 
tion  and  convection  heat -transfer  processes  separately,  we  retained 
the  nonlinear  boundary  conditions.  This  complicates  the  calculation 
because  temperature  levels,  as  well  as  temperature  differences,  must 
then  be  considered.  Recall  that  when  linear  heat-transfer  terms  are 
used,  differences  only  need  be  considered,  and  the  initial  or  environ¬ 
mental  temperature  can  be  chosen  as  zero.  This  restriction  reduces 
the  generality  of  the  solution.  Consideration  of  Equation  (A-33)  and 
Table  A.  1,  shows  that  the  environmental  temperature  t  ,  maximum 

values  of  radiation  input  i  ,  thermal  conductivity  k,,  and  thickness 

r  max  1  f 

L  of  the  fabric  must  be  considered  to  properly  select  the  non-dimen¬ 
sional  environmental  (and  initial)  temperatures. 
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Dimensionless  Temperature, 


Figure  A-8a  .  Temperature  Response  of  a  Semi-Infinite  Slab  (with  Front 

Surface  Heat  Loss)  Irradiated  by  a  Time -Dependent  Heat  Pulse 
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0.  1 


i.  r 


’0.0  20  40 


Dimensionless  Time, 


Figure  A-Sb.  Temperature  Res?.:rse  of  a  Semi-Infinite  Slab  (with  Front 
Surface  Heat  Loss/  Irradiated  by  a  Time -Dependent 
Heat  I'uise 
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Dimensionless  Temperature, 


10,  0 


0.01 


Dimensionless  Time,  0, 


Figure  A-$c 


Temperature  Response  of  a  Semi-Infinite  Slab  {with 
Front  Surface  Heat  Loss)  Irradiated  by  a  Time- 
Dependent  Heat  Pulse 
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In  addition,  the  non-dimensional  front-  and  rear -surface  radiation 
and  convection  heat-transfer  coefficients  must  be  specified  separately. 

Because  of  the  appreciable  time  involved  for  each  calculation, 
owing  to  the  relatively  limited  capacity  of  our  computer,  the  para¬ 
meter  variations  selected  for  the  computations  were  chosen  to  satisfy 
only  our  immediate  requirements.  Since  the  fabric-  and  skin-surface 
radiation  heat-transfer  coefficients  were  not  linearized,  and  also  were 
separated  from  the  convection  terms,  the  solutions  are  more  specific 
than  the  generalized  non-dimensional  solution.  Further,  temperature 
levels,  as  well,  as  temperature  differences,  must  be  considered. 

For  this  calculation,  the  following  conditions  were  chosen. 


E.  Case  V.  Opaoue  Fabric -Air space -Skin  System  (Variation  of 
Heat  Transfer  Coefficients 


1.  Initial  Conditions 


T]  =  T2  r  T3 . etc.  =  Te 

T  =  0.  108 


>  T 


e 


2.  Values  of  Constant  Parameters 


AX{  *  1/6 

AX  *  1/30 
s 
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2 


M,  =  M  = 
f  s 

AQf  =  AXf2/2  r  1/72 

3.  Values  of  Variable  Parameters  (Heat-Transfer  Coefficients) 


S 

1, 

10 

R. 

a 

1, 

10 

10 

Hr 

fo 

= 

0. 

01, 

0.  10 

H. 

lO 

s 

0, 

10, 

1 

Because  of  the  nonlinear  radiation  term,  a  non-zero  value  of  the  non- 
dimensional  environmental  temperature  must  be  chosen.  Referring  to 
Table  A.  1,  the  non-dimensional  representation  of  Tg  is 


e  i  L, 
max  f 

This  was  evaluated  by  selecting  a  constant  radiation  input  of 
2 

i  -  10  cal/cm  sec,  an  environmental  temperature  of  t  »  15*0, 

max  e 

and  assuming  a  9  oz  "standard"  cqtton-sateen  cloth  having  k  /L,  * 

°  2  *  ^ 

3.  76  x  10  °  cal/cm  sec  °C.  This  combination  resulted  in  a  T  of 

e 

0.  108. 

The  non-dimensional  heat-transfer  coefficients  (R,  ,  R.  ,  H,  , 

to  xo  to 

and  H;  )  were  selected  in  a  similar  manner,  using  the  best  available 
10 

information  pertaining  to  the  dimensional  convection  coefficients. 
Utilizing  the  above  selected  values,  the  radiation  coefficients  are  of 
the  order  of  5,  the  fabric  front-  and  rear-surface  convection  coef- 
ficents  are  of  the  order  of  0.  05  and  0,50,  respectively. 
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For  our  calculations,  only  envelopes  of  fabric  front-  and  rear- 
surface  and  skin  front-surface  response  are  presented  (Refer  to 
Figure  A-9).  They  show  the  effects  of  order -of -magnitude  changes 
of  the  transfer  coefficients,  selecting  as  a  narameter  the  extreme 
values  of  each  coefficient.  Variation  of  the  other  parameters  in  any 
manner  within  the  prescribed  limits  results  in  curves  lying  within 
the  envelopes,  A  presentation  such  as  this  allows  a  quicker  under¬ 
standing  of  the  significance  of  each  heat-transfer  parameter,  without 
the  necessity  of  considering  a  vast  number  of  curves  resulting  from 
incremental  changes  in  the  parameters. 

An  appreciation  and  an  understanding  of  the  curves  shown  in 
Figure  A-9  can  be  obtained  by  referring  to  the  Nomenclature  on  page 
A-47  and  the  non-dimensional  parameters  presented  in  Table  A.  1. 

The  dimensionless  representation  is  necessary  to  allow  implicit 
consideration  of  many  possible  temperature  responses  for  various 
heat  inputs  and  combinations  of  thermal  and  optical  properties,  and 
surface  heat  loss  rates.  The  exact  response  fer  a  particular  con¬ 
figuration  can  be  obtained  by  inserting  the  proper  thermal  and  optical 
properties  in  the  various  parameters. 

One  obvious  feature  of  these  curves  is  that,  when  the  fabric  front- 
surface  radiation  coefficient  (R^)  1  s  fixed,  the  envelope  for  AT^o  is 
very  small,  which  indicates  that  the  front-surface  temperature  depends 
largely  on  the  radiation  coefficient.  Other  conclusions  can  be  drawn 
as  well.  In  most  cases,  the  percentage  change  in  the  skin  response  is 
larger  than  the  fabric  responses,  except  for  -  1,  where  they  are 
approximately  equal.  The  size  of  the  envelope  increases  progressively 
from  the  fabric  outer  surface  to  the  inner  surface  and  skin  surface, 
except  for  H.  =  i,  where  they  are  of  the  same  order. 
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Figure  A-9a.  Envelopes  of  Temperature  Response  of  an  Opaque  Fabric 
Skin  System 
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Figure  A-9c. 
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Figure  A-9e. 
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Time  Increments 


Figure  A-9g. 
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SYMBOLS  AND  ABBREVIATIONS 


Symbol 

Bi 

c 

D 

F, 

f-8 

Fo 

h 

i 

J 

k 

L 

M. 

m 

n 

t 

x 

a 

r 

A 

€ 

P 

a 

T 


Biot  number 
specific  heat 
diathermancy  factor 

radiation  interchange  factor  between  fabric  and 
skin 

Fourier  number 

heat  transfer  coefficient 

radiant  heat  flux 

Jacobian 

heat  conductivity 

fabric  or  air  gap  thickness  or  characteristic  length 

dimensionless  modulus 

iteration  correction  factor 

iteration  correction  factor  or  series  term 

temperature 

distance 

absorptivity  or  diffusivity 
extinction  coefficient 
time  or  distance  increment 
emissivity 

density  or  reflectivity 
Stefan -Bolt*,  mann  constant 
time  or  transmissivity 
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Subscript 


a 

e 

eff 

r 

1 

j 

k 

max 

o 

s 


air  gap  or  air 
environment 
effective 
fabric 

inner  surface  of  fabric  or  initial  value 
general  space  increment 
last  fabric  element 
maximum 

outer  surface  or  impinging  radiation 
skin 


Superscript 

'  conditions  at  end  of  time  increment 

k  first  approximation  in  iteration 

n  general  time  increment 
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I.  INTERFEROMETER  THEORY 


In  the  Mach-Zehnder  instrument,  shown  in  Figure  B-l,  the  wave 
fronts  from  two  separated  beams  of  light,  both  emanating  from  a  single 
source,  are  made  to  interfere  with  one  another.  The  interference 
phenomena  produces  bright  or  dark  bands  on  the  image  plane  as  re¬ 
flected  by  phase  cancellation  or  amplification  of  the  beams  wherever 
the  wave  fronts  are  in  or  out  of  phase.  If  the  light  beam  passing 
through  the  test  region  of  interest  is  retarded  (or  advanced)  because  of 
changes  in  the  refractive  index  in  that  region,  the  interference  fringes 
can  be  seen  to  shift.  Localized  fringe  shifting  occurs  near  a  heated 
body  since  the  temperature  or  density  gradient  in  a  thermal  boundary 
layer  surrounding  the  heated  body  produces  a  direct  change  in  the 
refractive  index.  In  gases  the  simple  relationship  rj  -  1  =  C  ,  as  shown 
by  Lorenz- Lorentz  describes  the  dependency  of  density  p  on  the  index 
of  refraction  rj.  The  Gladestone-Dale  constant  C  maintains  its  con¬ 
stancy  for  gas  over  a  wide  range  of  pressure  and  temperature. 

Since  the  degree  of  fringe  displacement  is  dependent  on  total  path 
length  through  which  the  refractive  effects  are  occurring,  among  other 
things,  an  equation  or  equations  can  be  derived  expressing  the  varia¬ 
bles  to  whatever  determinant  or  quantity  desired.  For  example,  to 
calculate  temperature  at  any  location  where  refractive  variances  are 
restricted  to  a  two-dimensional  plane  (no  variation  occurs  along  the 
light  path)  the  following  equation  may  be  used: 


t 


2 


__1  _ 

5N  VTK*  tj 

—£ZP~ 


where: 

tj  =  ambient  temperature 
6N  -  fringe  shift 

X  =  wavelength  of  light 
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c 


Gladstone-Dale  constant 


Z  =  interference  path  length 
P  =  pressure. 

The  derivation  of  the  basic  interferometer  equations  are  presented  in 
Reference  1.  Included  also  is  a  discussion  of  measurement  errors.  A 
complete  description  of  the  interferometer,  its  light  source,  alignment 
procedures,  related  equipment  and  system  of  measurements  is  also 
included  in  the  above  reference. 


II .  DEVELOPMENT  OF  EQUATIONS 


A  geometrical  understanding  of  wave-front  interference  phenomena 
can  be  obtained  for  the  condition  of  ’’finite"  fringe  adjustment  of  the 
interferometer.  Consider  that  one  of  the  beam-splitters  has  been 
rotated  by  a  small  angle  in  a  given  plane.  The  emergent  coherent 
beams  are  now  not  parallel  but  intersect  such  that  phase  differences  in 
the  two  beams  appear  as  parallel  dark  or  light  fringes  everywhere 
along  the  intersection  plane  (see  Figure  B-2).  The  width  between 
fringes  will  depend  on  the  angle  of  rotation. 

For  analysis  purposes,  the  wave  fronts  can  be  made  to  intersect 
so  that  horizontal  fringes  will  be  projected  across  the  field  of  view; 
this  creates  a  two-dimensional  representation  of  fringe  shifts  in  the 
x-y  plane  permitting  easier  study.  With  appropriate  adjustment  of  the 
optical  elements,  the  intersection  of  wave  fronts  can  also  be  made  to 
occur  within  the  test  region  where  refraction  effects  take  place  (see 
Reference  1  for  alignment  procedures).  Analysis  is  made  with  the 
assumption  that  variations  in  the  z  direction  do  not  exist. 

For  a  heated  body  the  fringes  are  observed  to  bend  as  shown  in 
Figure  B-3.  The  direction  of  fringe  bending  (up  or  down)  depends  on 
the  rotation  direction  of  the  beam-  splitter.  The  optical  path  difference 
6Q  for  two  light  paths  (projected  at  points  1  and  2)  traveling  in  the 
z  plane  a  distance  Z  is  given  as 

5Q  =  (^  -  n2)  Z  , 

where  Z  is  the  length  of  heated  body  (z  direction),  and  and  are 
the  indices  of  refraction  of  gas  along  the  undisturbed  path  and  along  the 
shifted  or  deviant  path,  respectively.  The  difference  ir.  wavelengths 
&N  from  point  1  to  point  2  can  be  shown  equal  to’ 


Figure  B-2,  Wave  Front  Interference  and  Formation  of  Fringes 


5Q 
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where  X  is  the  wavelength  of  light  in  a  vacuum. 

Considering  the  fringe  bending,  point  2  represents  one  whole  fringe 
shift  from  the  uniform  undisturbed  field  at  point  1,  which  indicates  that 
the  optical  path  at  point  2  is  1  wavelength  different  from  that  at  point  1. 
The  number  of  wavelengths  reaching  each  point,  for  light  traveling  dis¬ 
tance  Z,  is  given  as: 


where  X^  and  X^  represent  the  wavelengths  corresponding  with  light 
traveling  the  undisturbed  and  the  deviant  paths,  respectively.  The 
difference  in  the  number  of  waves  is  expressed  as: 


5N 


Since  the  index  of  refraction  is 


1 


1 


The  wavelength  difference  can  then  be  written  as 

8N  \  (1,  -  v  =  5-£  . 


Independently,  Lorenz  (using  the  elastic-solid  theory)  and  Lorentz 
(using  the  electromagnetic  theory)  arrived  at  the  same  theoretical  law 
related  density  of  ga  &  with  the  refractive  index. 
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where  the  constant  depends  upon  the  molecular  nature  of  the  gas. 

Since  most  gases  have  a  refractive  index  very  nearly  1,  the  factor 
2 

(t|  +  1 )/( t|  +  2)  is  nearly  constant  ann  the  above  equation  can  be  simpli¬ 
fied  to: 
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n-  i 
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where  C  (Gladstone-Dale  constant)  maintains  its  constancy  for  a  gas 
over  a  wide  range  of  pressure  and  temperatures. 

Between  points  1  and  2  in  our  interference  field 


or 


(Tll  '  n2'  "  C  (pl  ~  P2)  • 

The  optical  path  difference  6Q  may  be  written  as: 

5Q  =  Z  ( -  iJ2>  CZ(Pi  -  P2) 

and  the  wavelength  difference  or  fringe  shift  number  5N  ie  then  related 
to  densi.y  as 


CZ 

nr 
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Fringe  shift  number  can  be  measured  at  any  point  in  the  field 
directly  from  the  interferogram.  For  density  measurement  at  the  face 
of  a  heated  body  (point  A  in  Figure  B-3),  fringe  crossings  can  be 
counted  from  the  undisturbed  field  along  a  line  parallel  to  the  undis¬ 
turbed  fringes  with  interpolation  made  between  adjacent  fringes  inter¬ 
cepting  the  surface. 

Fringe  count  can  also  be  made  from  A  to  B,  which  gives  the  same 
density  value  at  A  as  when  counting  directly  across  from  the  undis¬ 
turbed  field  at  1.  The  density  at  B,  however,  may  be  different  from 
that  at  A  and  a  separate  determination  is  necessary  at  that  point. 

Since  the  difference  in  one  fringe  shift  is  equivalent  to  a  difference 
of  one  wavelength,  another  technique  is  to  take  the  ratio  of  physical 
distance  Y,  e.g.  ,  from  A  to  B,  to  actual  fringe  width  b,  which  gives 
the  equivalent: 
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The  measurement  in  this  instance  is  for  point  B.  The  advantage  of 
this  technique  is  that  '  r  adjustment  of  widely  spaced  fringes,  the 
physical  measurements  of  Y  and  b  can  be  made  more  accurately  than 
with  narrow  fringe  adjustment.  (.  lote  that  fringe  count  would  be  the 
same  regardless  of  fringe  spacing.  )  The  choice  of  fringe  spacing 
depends  on  the  degree  of  refraction  occurring  in  the  test  region  and  on 
the  field  of  view  available  for  measurements. 

Since,  in  our  measurements,  we  are  actually  interested  in  the 
temperature  boundary  field  near  a  heated  object  rather  than  density, 
the  gas  law  equation  can  be  applied.  Any  convective  flow  of  gas  near 
the  heated  body  is  of  sufficiently  low  velocity  that  pressure  can  be 
assumed  constant  throughout,  and  //e  arrive  at  the  expression: 
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or,  in  terms  of  the  unknown  temperature  t£ 
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where  R  is  the  gas  constant  and  P  the  ambient  pressure.  For  a  given 
gas  and  for  monochromatic  light  of  given  wavelength,  XR/C  is  constant. 
Ambient  pressure,  temperature,  and  body  length  are  measured  to  solve 
the  expression  for  unknown  temperature  t,.  When  the  quantity 
(6N  XR  t j  )/(CZ  P)  is  small  compared  to  1,  it  is  convenient  for  compu¬ 
tational  purposes  to  convert  the  bracketed  expression  in  the  above 
equation  into  a  series  expansion 
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The  above  relationships  permit  a  localized  density  or  temperature 
determination  in  the  two-dimensional  plane  from  a  single  interferogram 
taken  while  the  refractive  variances  ore  occurring.  The  undisturbed 
interference  field  to  the  right  of  the  heated  body  is  the  reference 
region  of  uniform  density.  However,  if  cne  test  region  is  enclosed  and 
pressure  or  temperature  changed  uniformly  throughout,  the  region  of 
uniform  density  cannot  then  be  taken  as  an  absolute  reference.  The 
entire  fringe  field  will  shift  and  reference  must  necessarily  be  made  to 
the  optical  path  length  in  the  external  compensating  beam. 

If  the  fringe  shifc  in  the  uniform  density  field  can  be  observed  and 
measured  (assume  that  a  temperature  change  occurs  in  the  test  cham¬ 
ber),  then  temperature  can  again  be  determined  by  the  equation: 
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where  the  reference  t  ,  is  itself  determined  from  the  same  equation; 

r 

thus, 
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where  (&N)c  is  the  fringe  shift  in  the  uniform  field  with  respect  to  the 
compensating  optical  path,  t  is  the  external  temperature  in  the  com¬ 
pensated  path,  and  is  the  path  length  through  the  test  chamber. 

With  monochromatic  light  the  fringe  shift  in  a  uniform  field  is 
difficult  to  measure  since  all  fringes  are  alike,  indistinguishable  one 
from  another.  The  use  of  white -light  makes  this  measurement  possible 
since  each  fringe  can  be  identified  and  a  shift  in  the  uniform  density 
field  can  easily  be  determined. 

The  interferometer  equation  for  predicting  temperature  can  be 
plotted  for  different  plate  lengths  Z.  From  these  plots  (and,  as  shown 
by  equations)  the  sensitivity  and  accuracy  of  the  interferometer  to 
predict  high  temperature  decreases.  From  the  above  developed  equations 
when 
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the  predicted  temperature  would  become  infinite.  The  fringe  shift  number 
6N  has  a  discrete  maximum  value  for  a  given  plate  length  Z. 
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III.  DISCUSSION  OF  ERRORS 


Deviation  from  100  percent  correlation  is  shown  for  the  fringe 
shift  temperature  prediction  versus  the  thermocouple -measured 
temperature  in  Figures  B-4  and  B-5.  In  practice  the  calibration 
curve  would  suffice  to  determine  interferometric  fringe -shift  temper¬ 
ature.  However,  since  errors  such  as  end  effects  anci  refraction  effects 
have  been  treated  extensively  in  the  li'erature,  these  and  other  errors 
must  be  considered. 

A.  Fringe  Shift  Measurement 

The  scatter  of  data  points  from  a  mean  datum  line  is  attributed 

to  fringe  count  accuracy.  Our  measurements  were  made  by  projecting 

the  original  interferogram  negative  onto  a  screen.  From  the  enlarged 

interference  field  (25-50X  scale)  extrapolations  between  fringe  lines  were 

made  with  a  ruled  scale.  Variation  between  separate  readings  or  between 

readings  made  by  different  individuals  was  ±0.  1  fringe  with  a  mean 

variation  of  approximately  0.  05  fringe.  An  error  of  ±0.  1  fringe  for  a 

total  shift  of  5  fringes  would  give  ±2  percent  deviation,  and  for  a  shift 

of  a  single  fringe  the  deviation  would  be  ±10  percent.  Better  absolute 

accuracy  could  be  made  by  using  a  microscope  and  mechanical  stage 

2 

arrangement  to  correct  for  projection  distortion.  Kendall  describes 
an  interferometer  fringe  pattern  analyzerused  for  routine  measure- 

3 

ments  which  has  an  accuracy  of  at  least  0.  1  fringe.  Werner  devised 
a  means  to  measure  fringe  shift  with  indicated  accuracy  of  better 
than  1/ 500th  of  a  fringe  separation. 

B.  Ambient  Condition  Errors 

As  shown  in  the  preceding  equations,  relation  of  fringe  shift  data 
to  temperature  measurement  requires  knowledge  of  the  ambient  pressure 
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(Measurements  made  on  heated  silver 
disk  1. 9  cm  x  3. 2  cm  x  0.  3  cm! 


Measured  Fringe 
No.  Versus 
T  emperature 
Variation 


Theoretical  Fringe  No.  Versus 
Temperature  Variation 
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Figure  B-5.  Variation  of  Fringe  Shift  Number  with  Temperature, 
Measured  and  Theoretical  (3.  2  cm  Plate  Length) 
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and  temperature.  The  effects  of  variations  of  each  for  a  given  plate 
width  {3.  18  cm),  are  shown  in  Figures  B-6  and  B-7. 

C.  Refraction  Effect  Errors 

Light  passing  through  the  thermal  boundary  field  near  a  heated 
object  is  bent  away  from  the  surface  toward  the  undisturbed  region  of 
greater  density.  Since  the  path  deviates  through  layers  of  increasing 
density,  the  average  density  through  which  the  beam  passes  is  higher 
than  it  would  be  if  refraction  did  not  occur.  Fringe  shift  count  would  be 
less  Since  the  difference  in  density  between  the  undisturbed  region  and 
the  averaging  density  is  less. 

In  aerodynamic  studies  refraction  errors  were  considered  by 

4  5  6 

Gooderum,  et.  al. ,  Blue.  Ladenburg  and  Bershader,  and  Howes 

7 

and  Buchele,  among  others.  Gooderum  defines  the  problem  quite 
well  and  for  his  study  of  turbulent  mixing  of  a  free  supersonic  jet 
found,  fT-om  developed  expressions,  that  the  effects  of  refraction  were 
negligible.  Blue  suggests  that  a  wind-tunnel  width  dimension  for  a 
particular  boundary  layer  study  be  such  as  to  make  the  sum  of  re¬ 
fraction  error  directly  proportional  to  the  square  of  tunnel  span. 
Ladenburg  and  Howes  and  their  associates  present  very  elegant  ex¬ 
pressions  to  correct  for  refraction  effects  applicable  to  wind  tunnel 
conditions. 

g 

In  heat  transfer  studies  Kennard  did  not  mention  refraction  errors 

9 

or  may  not  have  been  aware  of  them.  Eckert  and  Carlson  in  early 
laminar -free  convection  studies  did  not  consider  the  refraction  errors, 
but  in  a  later  report  Eckert  and  Soehngen^  derived  a  correction  for 
light  ray  refraction  in  determining  heat  transfer  coefficients  around 
cylinders  in  cross  flow.  The  temperature  correction  was  given  as 


B-16 


100  130  200 

Temperature  tt^l  (  0 


Figure  B-6.  Variation  o f  Fringe  Shift  Number  with  Temperature 
for  Changes  in  Ambient  Pressure 
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where  C  is  the  Gladstone -Dale  constant,  Z  the  object  length,  P  refer¬ 
ence  pressure,  R  gas  constant,  and  t,  rtference  temperature.  The 

(dt  I  ^ 

^•1  is  determined  from  the  temperature  boundary 
profile  as  shown  in  Figure  B-8. 


Figure  B-8.  Temperature  Gradient  at  Wall 

Goldstein  and  Eckert*1  used  the  derivation  developed  by  Eckert 
and  Soehngen10  and  found  that  the  refractive  error  was  less  than  1/2 
percent  in  his  free  convective  boundary  layer  measurements  on  a 
uniformly  heated  vertical  plate.  The  length  dimension  of  his  plate 
was  approximately  10  cm.  Carlson1^  applied  the  same  correction 
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equation,  and  for  a  temperature  difference  of  71  C  across  an  air  gap 
between  two  30  cm-wide  plates,  calculated  a  negative  2  percent  error. 
Our  calculations  for  refraction  errors,  using  the  above  relationship, 
show  that  for  even  large  temperature  gradients,  refraction  error  is 
very  small.  The  refraction  errors  may  be  significant  if  Z  is  very  large, 
5  t  is  very  large,  and  a  very  thin  boundary  layer  exists  (forced  con¬ 
vection  flow). 

D.  End  Effect  Errors 

The  most  significant  error  in  funge  shift  prediction  of  temperature 
for  our  investigations  is  the  effect  of  increased  path  length  due  to 
heated  air  extending  beyond  the  ends  of  the  specimen.  The  effect  of 
length  on  fringe  shift  count  for  a  given  temperature  is  shown  in 
Figure  B-9.  Compared  with  variation  in  ambient  temperature  or  am¬ 
bient  pressure,  Figures  B-7  and  B-6,  the  length  characteristic  is 
considerably  more  significant. 

Before  an  end  correction  can  b~  appropriately  applied,  some 
consideration  of  the  boundary  layer  at  the  ends  of  a  heated  specimen 
must  be  made.  The  interferograms  in  Figures  B-10  and  B-ll,  show  iso¬ 
therms  about  a  heated  disc.  It  can  be  imagined  that  the  boundary 
layer  at  the  ends  would  be  somewhat  similar  to  that  showing  at  the 
bottom  of  the  disc.  The  boundary  layer  would  increase  in  thickness 
with  vertical  height  as  the  convective  flew  rises  and  the  correction 
would  be  different  from  the  bottom  to  the  top  of  the  disc  (see  Figure 
B-12). 

In  our  measurements  the  indicated  fringe  shift  temperature  was 
indeed  higher  along  the  vertical  axes  with  fringe  count  increasing  0.  I 
to  0.  i  fringes  from  the  bottom  quarter  to  the  middle  of  the  disc.  From 
th  middle  to  the  top  of  the  disc  the  nc rease  was  less  (0,  1  fringe).  An 
effort  was  made  to  determine  whether  any  temperature  variations 
existed  along  the  plate  by  providing  a  differential  ;hcrmoeouple  circuit 
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Figure  B-9.  Variation  of  Fringe  Shift  Number  with  Temperature 
for  Various  Specimen  Lengths 
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Constant  Temperature  Contours  for  Vertically  Heated  Plate 


to  measure  temperature  difference  from  center  to  edge  of  the  disc. 
Although  the  tests  we  re  inconclusive,  indications  were  that  gradients, 
if  existing,  were  extremely  3mall  because  of  the  high  thermal  dif- 
fusivity  of  silver.  More  significant  gradients  might  be  encountered 
with  fabrics  or  materials  of  low  conductivity. 

With  increased  temperature, or  with  time, the  boundary  layer  thick¬ 
ness  and  end  corrections  would  not  remain  the  same  at  any  path  position. 
This  time -temperature  effect  would  certainly  be  true  for  transient  heating 
experiments,  and  probably  not  as  time -dependent  when  measurements  are 
made  during  free  convective  cooling.  Under  forced  convection  conditions 
the  entire  boundary  layer  could  be  radically  changed. 

g 

Kennard  made  several  assumptions  regarding  the  end  corrections. 
Without  allowance  for  end  effects  he  found  that  the  interference  tempera¬ 
ture  indication  was  approximately  5  percent  higher  than  the  temperature 

in  the  plate  itself.  His  correction  as  one  assumption  amounted  to  simply 

9 

ascribing  a  5  percent  increase  in  length.  Eckert  and  Carlson  made  a 

reasonable  approximation  of  the  end  effects  by  the  method  of  conformal 

11  12 

mapping.  Goldstein  and  Eckert  and  Carlson  all  considered  end  ef¬ 
fects  by  the  method  proposed  by  Eckert  and  found,  for  their  conditions 
of  long  interference  path  (10  cm  to  30  cm)  and  thin  plate,  that  end 
corrections  were  small  and  tended  to  cancel  refraction  errors  of 
opposite  magnitude.  Bevans^  failed  to  arrive  at  any  quantitative 
correction  for  end  effects  in  his  experiments  of  vertical  free  convec¬ 
tion  from  an  isothermal  surface.  All  of  the  above  corrections  were 
made  for  steady -state  or  near  steady -state  conditions. 

Experimental  and  theoretical  investigations  of  the  free  convection 
boundary  layer  growth  with  time  on  semi -infinite  vertical  flat  platet. 
under  different  conditions  of  heat  input  are  described  in  the  review  by 
Goldstein.  We  have  yet  to  correlate  our  experimentally  produced 
boundary  profiles  with  prediction  by  others.  The  variations  with 
time,  temperature,  and  vertical  location  for  a  given  heat  input 
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function  describe  the  difficulty  in  assuming  an  end  correction  based  on 
the  boundary  layer  or  thermal  boundary  condition  at  the  end  of  the  plate. 

E.  Other  Errors 

Other  errors  are  inherent  in  the  system.  For  example,  alignment 
of  the  test  specimen  must  oe  perfect  with  reppect  to  the  collimated  light 
beam,  or  surface  fringe  shift  measurements  camot  be  made  (alignment 
can  be  made  perfect  only  if  the  source  is  infinitely  small  and  the  colli¬ 
mating  lens  perfectly  corrected  for  the  wavelength  of  light).  Also,  refer¬ 
ence  temperature  in  the  undisturbed  region  must  be  correctly  known  or 
readily  determined.  Changes  of  the  refractive  index  due  to  foreign  gases 
have  been  mentioned  as  error  possibilities.  Heat  flux  uniformity  through 
the  flux  redistributor  is  not  optimum  and  non-adherence  to  two- 
dimensionality  causes  averaging  errors.  Heated  air  within  the  dis¬ 
tributor  may  generate  convection  flow  irregularities.  Aberrations  in  the 
interferometer  optics  or  lens  can  distort  fringe  patterns  significantly. 
Enlargement  of  interference  pictures  may  give  distortion  errors. 

There  are  other  possible  errors  but,  with  due  precaution,  these  and 
the  above  sources  of  error  can  be  reduced  or  made  negligible. 


IV.  INTERFEROMETFR  LIGHT  SOURCE  AND  PHOTOGRAPHIC 

REQUIREMENTS 

In  optimizing  performance  characteristics  of  the  interferometer 
various  undisturbed  fields  or  no  flow  tests  were  carried  out  using  dif¬ 
ferent  light  sources  and  optical  filters.  Unfortunately  the  requirements 
of  monochromatic  light  for  optimum  fringe  patterns  and  high  intensity 
light  for  short  duration  exposures  tend  to  be  conflicting  thus  simul¬ 
taneous  optimization  is  extremely  difficult. 

A  Pek  Model  200-watt  high  pressure  mercury  light  source  met  the 
necessary  light  intensity  and  spectral  requirements  and  subsequent 
interference  studies  were  made  with  this  source.  This  lamp  is  capable 
only  of  continuous  operation  and  has  an  effective  source  size  of  0.  25  cm 
x  0.  18  cm. 

In  conjunction  with  the  mercury  source,  nearly  monochromatic 
light  was  obtained  bv  using  a  Baird  Atomic  Corporation  5641  A  typebl 
interference  filter.  With  the  interference  filter  90%  of  the  light  in-4 
tensity  was  lost  either  in  the  basic  filtering  process  or  in  filter  ab¬ 
sorption  or  reflection  in  its  bandpass  range.  Other  factors  reduce 
the  lightintensity  further,  such  as  losses  in  the  interferometer  optics 
(one  beam  splitter  discards  50%  of  the  light  in  each  beam),  the  small 
effective  source  size  required  and  the  magnification  of  image  desired. 
Tho  effective  source  size  (pinhole  diameter)  was  varied  and  fringe 
quality  found  tv  improve  by  decreasing  the  hole  diameter.  Smaller 
pinhole  sizes  than  0. 1  centimeters  did  not  appreciably  improve  the 
fringe  rendition. 

The  light  intensity  at  the  image  plane  was  measured  with  a  light 
dependent  resistor  having  a  negative  coefficient  of  resistance  linear 
over  an  extremely  wide  range  of  light  intensities.  The  light  intensity  - 
resistance  values  for  different  films  were  determined  and  relation¬ 
ships  for  image  magnification,  pinhole  sise  and  shutter  speed  ex¬ 
pressed. 
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V.  TEMPERATURE  CALIBRATION  OF  INTERFEROMETER 

Tests  were  conducted  to  calibrate  and  verify  experimentally  the 
ability  of  the  Mach-Zehnder  interferometer  to  indicate  temperature 
directly.  The  investigation  of  a  free  convection  boundary  layer  on  a 
heated  plate  is  a  convenient  way  to  check  the  validity  of  data.  A 
calorimeter  with  front  and/or  rear  surface  mounted  thermocouple 
was  positioned  at  the  flux  redistributor  exit.  These  discs  were  made 
of  silver  and  had  widths  corresponding  to  the  actual  fabric  test  pieces 
(3.  2  centimeters,  3.  8  centimeter  wide).  Interferograms  were  taken 
of  the  arc  heated  disc  at  different  temperatures  and  fringe  shift  cal¬ 
culations  of  surface  temperatures  were  compared  with  thermocouple 
measurements. 

For  interference  fringe  shift  measurement  of  temperature  the 
undisturbed  field  behind  the  calorimeter  was  taken  as  an  ambient 
reference  temperature.  The  field  between  the  top  and  the  bottom  re- 
distributor  plates  was  found  to  seek  a  higher  temperature  level  after 
each  e  posure  and  was  not  used  as  ambient  reference.  Representa¬ 
tive  interferograms,  Figure  B-13,  show  silver  discs  at  different 
temperatures.  Figure  B-14  shows  a  comparison  of  fringe  shift  tempera¬ 
ture  prediction  of  3.8  centimeter  disc  with  the  thermocouple  measure¬ 
ments  for  temperatures  to  150'  C.  Figure  B-15  shows  a  correlation 
comparison  for  a  3.  2  centimeter  disc  up  to  3^0CC  in  terms  of  fringe 
shift  number. 


Figure  B~13a.  Interferogram  of 
Thin  Heated  Platt 
(275  C) 


Figure  B-13b.  Interferogram  of 

Thick  Heated  Plate 
(350  a 


Figure  B-13c.  toterferogram  of  Two 
Thin  Heated  Plate# 
(Front  Plate  375  Q 


Figure  B-13d.  Interferogram  of  Thick 
and  Thin  Heated  Plate# 
Infinite  Fringe 
(Front  Plate  325  C) 
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NOMENCLATURE 


Symbol 


b 

C 

N 

BN 

(BN)C 

n 

p 

BQ 

R 

t 


x,y,z 


P 

X 


fringe  width 

Gladstone -Dale  constant 

number  of  wavelengths 

wavelength  difference  or  fringe  shift  count 

fringe  shift  referenced  to  compensated  optical  path 

refractive  index 

ambient  pressure 

optical  path  difference 

gas  constant 

temperature 

rectangular  coordinates 

fringe  shift 

length  of  heated  body 

length  of  test  chamber 

density 

wavelength  o£  light  in  vacuum 


Subscript 

1,2  location 

r  reference 
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GENERAL 


Quantitative  evaluation  of  fringe  shift  data  yields  an  accurate  indi¬ 
cation  of  temperature  provided  additional  information  is  available  during 
release  of  degradation  gases.  Among  various  constants  and  physical 
parameters  that  enter  into  the  equations  the  gas  constant  and  variation  of 
index  of  refraction  with  density  must  be  known. 

As  fabric  absorbs  high  intensity  thermal-radiation  a  condition  is  soon 
reached  during  which  degradation  processes  start.  Changes  are  evidenced 
by  release  of  various  degradation  products  accompanied  by  both  exothermic 
and  endothermic  heat  exchanges.  These  reactions  are  sensitive  not  only 
to  temperature  level  but  to  temperature  change  rates  as  well.  Knowledge 
of  the  various  gas  concentrations  in  the  free  convective  boundary  layers 
adjacent  to  the  heated  fabric  is  required  to  allow  the  assignment  of  ef¬ 
fective  gas  properties  to  be  used  in  the  temperature  evaluation  equations. 

From  preliminary  motion  picture  sequences  it  was  observed  that 
gaseous  products  were  emitted  from  the  fabric  at  various  instances, 
particularly  during  the  final  heating  phase.  These  gases  showed  as 
light  or  dark  areas  in  the  interference  pattern  which  would  develop 
initially  at  the  fabric  surface  and  subsequently  move  out  into  the  boundary 
layer  surrounding  the  fabric.  Although  no  sharp  discontinuities  were 
displayed  in  the  interference  fringe  patterns,  it  was  recognized  that 
the  actual  fabric  temperature  or  temperatures  i  lhe  thermal  boundary 
layer  differed  from  the  indicated  values  due  to  the  presence  of  gas 
having  an  index  of  refraction  different  from  that  of  air.  It  was  also 
recognized  that  if  the  reactions  occurring  have  any  resemblance  to 
chemical  temperature  effects  of  known  constituents,  the  interferometer 
technique  would  be  ideal  for  indicating  composition. 

Two  interferometer  techniques  were  considered  to  determine  the 
extent  to  which  interference  fringe  shifts  temperature  data  varied  from 
actual  values.  A  two  or  more  wavelength  technique  and  a  multi- 


environmental  experimental  technique  were  considered,  the  latter  ex 
tensively  utilized  in  this  present  program. 
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MULTI-WAVELENGTH  INTERFEROMETER  TECHNIQUE 


II. 

A,  General 

Olsen*  rtpi.aca  this  interferometric  method  to  measure 
chemical  composition  and  temperature  of  a  multi -component  gas 
mixture.  This  can  be  accomplished  even  when  chemical  changes  occur 
in  the  gas  mixture.  As  Olsen  points  out,  the  application  of  this 
method  to  combustion  processes  is  complex  because  of  the  many  gas 
components  present  and  because  of  the  high  temperatures  and  temper¬ 
ature  gradients  encountered.  No  information  was  found  as  to  how  the 
refractive  index  for  different  gases  varies  with  temperatures.  In  his 
study  of  steH&y -state  combustion  occurring  in  Bunsen  flames..  Olsen  first 
worked  out  a  derivation  to  apply  to  a  binary  gas  mixture.  His  analysis 
was  limited  to  a  two-dimensional  air -oxygen  jet,  but  indications  show 
that  a  more  complex  mixture  can  be  studied. 

2,  3 

Ross,  El-Wakil  and  C.  Jaeck  used  the  two  wavelength  interfer¬ 
ometric  technique  for  the  study  of  vaporization  and  combustion  of  fuels 
and  indicated  the  usefulness  of  this  method  for  heat  and/or  mass  trans¬ 
fer  investigations. 

In  our  studies  such  a  technique  appeared  very  useful.  Initially, 
a  water -saturated  fabric  was  heated  with  an  arc  furnace  in  which 
release  of  water  vapor  associated  with  the  heats  of  vaporization 
could  be  studied  It  was  conceived  that  the  water  v/ould  be  evaporated 
from  the  fabric  surface  at  or  somewhat  above  100°C.  A  period  of 
constant  temperature  exists  during  vaporization  and  until  the  water  is 
fully  released,  as  evidenced  by  the  temperature  plateaus  in  the  temper¬ 
ature-time  recordings.  The  indicated  fabric  surface  temperature 
during  vaporization  was  approximately  150  to  ITS'C. 

To  evaluate  the  two  or  more  wavelength  technique  a  more  rigorous 
investigation  of  refractive  variances  associated  with  the  release  of  water 
vapor  or  other  foreign  gases  was  conducted. 
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B.  General  Temperature  Relationship 


The  relationship  expressing  temperature  changes  in  air  surrounding 

a  heated  fabric  in  terms  of  fringe  distortion  has  been  worked  out  for  our 

4 

particular  studies.  The  expression  is  given  as 
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The  reference  density  P  ^  is  determined  from  the  perfect  gas  relationship 


pi  *  slf 


(C-2) 


where  pressure  p  and  temperature  tj  are  easily  measured  in  the  un¬ 
disturbed  field. 

C.  Constant  of  Proportionality 

The  basic  expression^relating  the  index  of  refraction  to  gas  density 
tor  a  given  wavelength  of  light  includes  a  constant  of  proportionality  C 
that  would  be  different  for  the  different  wavelengths  of  light  used  in  the 
multi -wave length  technique. 

Since  this  number  is  important  it  can  be  considered  in  some  detail. 

The  proportionality  constant  is  arrived  at  through  the  general  dis¬ 
persion  expression  developed  by  Lorentz  (1936): 
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where  N  is  the  number  of  molecules  par  unit  volume  of  the  gas  and, 
therefore,  is  proportional  to  density.  The  remainder  of  terms  in  the 
summation  of  the  right-hand  side  of  the  expression  is  equal  to  a  con¬ 
stant  K,  relating  density  and  the  index  of  refraction,  and  includes  the 
terms  e  and  m,  electron  charge  and  mass;  V  ,  velocity  of  light  in 
vacuum;  X  ,  wavelength  of  light  in  vacuum;  and  X,  the  variable  wave¬ 
length  in  the  summation  of  dispersion  values.  For  one  wavelength  of 
light  and  for  one  medium  (air),  a  linear  approximation  of  the  Lorentz- 
Lorenz  equation,  based  on  the  fact  that  the  index  of  air  is  nearly  equal 
to  1,  yields  the  proportionality  constant  C: 


C  = 


(C-4) 


where  C  is  simply  3/2  K. 

This  equation  (Gladstone -Dale  expression  1858,  1863)  was  experi¬ 
mentally  verified  and  actually  preceded  the  Lorentz-Lorenz  relationship, 
which  corrected  deviations  resulting  when  a  liquid  passes  to  a  gaseous 
state.  The  constancy  of  C  or  K  for  the  Lorentz-Lorenz  expression 
is  markedly  superior  over  that  of  the  Gladstone -Dale  approximation  when 
the  pressure  or  density  of  gas  increases  greatly.  However,  in  the  present 
study  the  linear  expression  for  density  and  refractive  index  provides  a 
constant  of  proportionality  unchanged  for  the  range  of  conditions  encountered. 

The  exactness  in  the  measure  of  C  is  dependent  on  how  accurately 
can  be  measured  and  how  precisely  is  known  for  the  particular  wave¬ 
length  of  light.  For  our  tests5  C  was  calculated  from  tabulated  index 

5,6 

values  predicted  and  measured  by  Edlen  for  standard  conditions. 

Edlen’s  index  of  refraction  expression  is  given  as 


where  tj  is  the  ref- 
8 

measured  in  mic 


1)  108  .  6432.  8  ♦  M15L810  +  J«40 

146  -  -ij  41  -  -!j- 

tive  index  for  standard  air,  and  wavelength  X  is 


(C-5) 
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D.  Presence  of  CO2  in  Air 

Standard  air,  as  considered  in  the  above  expression,  is  dry  air  con¬ 
taining  0.  03%  CC>2  by  volume  at  normal  pressure  of  760  mm  Hg  and  having 
an  air  temperature  of  15*C.  Earlier  measurements  generally  neglected 
the  small  quantity  of  in  air. 

Laboratory  air  may  contain  considerably  greater  concentrations  of 
CO2  than  normal  air  designated  as  standard.  The  index  of  refraction  of 
CC>2  at  20”C  and  for  X*  0.  5461  |i  is  1,  00045060  and  for  air  at  the  same 
wavelength  and  temperature  is  1.00027310.  A  concentration  change  of 
CO2  in  air  from  0.  03  to  0.  3%  would  increase  the  air  index  value  to 
1.  00027364.  Thus,  for  exact  measurements  the  amount  of  CO2  present 
and  its  increase  in  confined  areas  where  refractive  studies  are  being 
made  should  be  considered. 

E.  Moist  Air  Conditions 

The  above  discussion  pertains  to  a  dry -air  condition  which  is  not 
generally  manifest  in  laboratory  experiments,  Lorenz  determined  a 
correction,  general  for  all  wavelengths  of  visible  light,  and  applied  it 
to  the  refractive  index  of  dry  air  n  to  obtain  the  index  for  moist  air  n  : 

\r>  *  ng  -  0.556  X  10’6  (pwv)  .  (C-6) 

Vhe  vapor  pressure  p^  is  given  in  mm  of  Hg.  From  the  table 
below  it  can  be  seen  that  the  influence  of  water  vapor  even  at  saturated 
conditions  is  only  felt  at  higher  temperatures  when  the  vapor  pressure 
becomes  a  substantial  portion  of  the  total  pressure.  Presumably  if  the 
relationship  holds  throughout  the  temperature  range  imposed,  at  100*C, 
when  the  water  vapor  constitutes  the  entire  gaseou»  phase,  the  differ¬ 
ence  in  refractive  index  is  the  difference  between  the  indices  of  water 
vapor  and  air. 
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Table  C-l.  Change  in  Index  of  Air  for  Saturated  Water  Vapor  Conditions 


Temperature 

•c 

T)  «  1.  00027773 

S 

Saturated 

Vapor  Pressure 
mm  Hg* 

Change  in  Index  at  Air 
due  to  Water  Vapor 
Af)(at  given  temperature) 

0 

4.  579 

0. 000000255 

15 

12.788 

0.000000711 

50 

92.51 

0.000005140 

100 

760.  00** 

0.  000052250 

♦values  obtained  from  steam  tables 
**100  percent  water  vapor 


E.  Effects  of  water  Vapor 

For  a  wetted  fabric  the  surface  temperature  history  after  initiation 
of  the  radiant  pulse  shows  a  temperature  plateau  near  100°C  where  the 
release  of  water  vapor  and  the  vaporization  heat  exchange  occurs.  If 
we  consider  that  saturated  water  vapor  exists  near  the  surface  of  the 
fabric,  the  change  in  indices  is  as  follows: 


%  "  ^  i  00 

T) 

100 


x  100 


^100  "  ^  wv 


wv 


x  100 


tj  (15' C  Dry  Air)-l.  0002777  \  „ 

t)  ( 100#C  Dry  Air)  -  1.  0002145/ 

T)  (100*C  Water  Vapor)  -  1. 00017222^  19.  6% 


Calculating  an  interferometric -indicated  temperature  under  these 
conditions,  the  actual  temperature  should  be  lower  by  as  much  as  70*C. 


G.  Dispersion  Values  of  Water  Vapor 


The  Lorenz  correction  earlier  stated  applies  uniformly  over  the 
visible  spectrum  and  is  based  on  the  consideration  that  the  dispersion 
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of  water  vapor  is  the  same  as  air.  In  reality  the  indices  of  refraction  vary 

7 

from  one  wavelength  to  another  for  different  gases.  Cuthbertson  and 

8 

Barrell  and  Sears  developed  expressions  governing  the  dispersion  of 
water  vapor  in  air -vapor  mixtures.  Barrell  and  Sears1  expression  (not 
included  here)  was  derived  for  moist  air  conditions  at  approximately  30* C 
and  relative  humidity  of  80%;  it  is  intended  for  use  at  normal  atmospheric 
conditions  where  vapor  pressures  of  10  mm  is  average  and  rarely  exceeds 
20  mm.  As  with  the  Lorenz  correction,  the  presence  of  water  vapor  at 
standard  conditions  changes  the  index  value  of  air  only  slightly.  Whether 
these  expressions  hold  for  the  possible  higher  concentrations  of  water 
vapor  associated  with  higher  temperatures  is  not  indicated. 

r 

H.  Measurement  of  Vapor  Concentrations  for  Dispersion  Values 

The  possibility  of  measuring  the  dispersion  of  water  vapor  and 
air -vapor  mixtures  to  arrive  at  concentration  profiles  of  the  water 
vapor  at  higher  temperatures  can  be  postulated.  The  change  of  re¬ 
fractive  index  from  pure  water  vapor  at  100*C  to  vapor -free  air  at  the 
same  temperature  and  wavelength  is  quite  substantial  as  indicated 
earlier.  However,  the  dispersions  of  gases  at  different  wavelengths 
are  small.  If  the  interference  events  are  recorded  at  different  wave¬ 
lengths  simultaneously,  the  indicated  temperatures  for  each  inter  - 
ferogram  would  be  different,  since  the  dispersions  from  X  ^  to  X  ^ 

•.re  not  equal: 

1 

The  subscripts  A  and  wv  denote  air  and  water  vapor,  respectively.  Using 
Barrell  and  Sears*  dispersion  formulation  for  air -water -vapor  mixtures, 
and  assuming  that  it  applies  at  100*C  and  for  higher  water  vapor  concen¬ 
trations,  one  can  ascribe  different  concentration  percentages  of  water 
vapor  in  the  air -vapor  mixture  until  the  simultaneously  recorded 
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interferograms  at  the  different  wavelengths  yield  the  same  temperature. 

At  that  point,  the  ascribed  concentrations  percentage  would  apply.  Since 
the  concentrations  and  temperatures  would  be  different  at  every  location  in 
the  thermal  or  diffusion  boundary  layers,  a  series  of  nomographs  re- 
lating  vapor -air  concentrations  (for  the  two  wavelengths)  with  temperature 
and  fringe  shift  number  v/ould  eliminate  the  necessity  of  solving  simul¬ 
taneous  expressions  for  every  location. 


I.  Accuracy  in  Predicted  Concentrations 


In  Table  C-2  the  dispersion  correction  coefficients  by  Barrell  and 
Sears  are  listed  (k  values)  for  various  wavelengths.  The  indices  of  dry 
air  at  15°C  and  100°C,  and  the  computed  vapor  indices  - 

k(Pwv)  ]  for  pure  water  vapor  (100%  saturation  of  water  vapor  in  air 
at  100°C)  are  listed.  As  an  example,  if  we  record  simultaneous  inter¬ 
ference  fringe  fields  using  the  Hg  lines  at  0.546lp  and  0.  4358  p  ,  we 
can  note  that  the  dispersion  for  air  at  100UC  is 

1.0002146  -  1.  0002170  =  0.0000034 
and  for  water  vapor  at  the  same  wavelengths 


1.  0001754  -  1.0001720  =  0.0000024  . 

The  difference  between  these  dispersion  values  in  terms  of  an  equivalent 
change  in  air  temperature  is  approximately  1.  5UC.  This  means  that  if 
we  had  100%  water  vapor  at  the  fabric  surface  but  assumed  a  dry  air 
index  relationship,  the  indicated  temperatures  from  both  interferograms 
would  differ  by  only  1.  5°C.  The  indic  ated  temperature,  however,  would 
be  higher  by  70°C  since  indices  of  air  rather  than  of  water  vapor  were 
used.  For  example,  a  measure  of  temperature  to  0,  1"C  from  each 
interferogram  would  still  result  in  an  approximate  6.  7%  error  in  the 
determination  of  concentration  or  in  temperature  difference  between 


Table  C-2.  Refractive  Index  Values  for  Air  at  15°C  and  100°C,  and  for 
Water  Vapor  at  100°C  at  Various  Wavelengths 
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indicated  and  actual  values.  Therefore,  the  accuracy  of  measuring 
fringe  shift  differences  between  interferograms  reflects  considerably 
not  only  on  the  accuracy  of  predicted  concentration  but  on  the  actual 
temperature  as  well. 

J.  Expressions  Governing  the  Two  or  More  Wavelength  Inter¬ 
ferometric  Techniques 

From  the  elastic  solid  and  electromagnetic  theories  relationships 
of  the  index  of  refraction  for  a  nonabsorbing  medium  to  density  and  gas 
constituents  has  evolved  (see  earlier  discussion). 


-  1 
+  2 


(C-8) 


3 

where  P  ./M.  are  the  concentrations  (moles /cm  )  assumed  fcr  the  con- 

ii 

stituents  and  is  the  molecular  refractivity  of  each  constituent  given. 


*  ~j~~  No.  (C-9) 

where  N  is  the  Avogadro  number  and  cr  (cro^)  is  a  frequency -dependent 
"polarizability"  factor  (electromagnetic  theory  of  light)  which  gives  sig¬ 
nificance  to  the  individual  atoms  and  molecules  in  relating  induced  dipole 
moment  and  the  effective  electric  field.  The  molar  refractivity  P  is  un¬ 
effected  by  a  change  of  state  and  is  independent  of  density. 

Utilizing  the  Gladstone -Dale  appr ox;mat:on  and  expressing  Equation 
(C-8)  in  terms  of  the  mol  fraction  (i.’M. )  cf  each  constituent 

(9-1)  *  j  E  Xj  Pj  (C-10) 

1 

where  x^  is  the  mol  fraction  (moles/gm)  for  the  i  component  and  P.  is  the 
refractivity  of  that  component. 
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The  relationship  between  fringe  count  and  change  in  indices  (from 
the  undisturbed  interference  field  to  the  surface  of  heated  specimen)  is 
given  for  the  one  wavelength  interference  event  as 


SN  «  -f-  (,,  -i)  2) 


(C-ll) 


where  Z  is  the  path  length  through  which  refraction  of  light  occurs  (width 
of  fabric)  and  x  is  the  wavelength  of  light.  The  index  value  for  the  un¬ 
disturbed  condition  is  given  as 

<n  -  *>  ■  I  ife  ?  hpi),  <c‘iz> 

1  i  1  *  l 

and  for  the  changed  or  disturbed  condition  as 


D  « 


I  Rt, 


(C- 13) 


If  these  expressions  for  and  tj are  substituted  into  Equation  (C-ll) 
the  fringe  shift  equation  for  the  one  wavelength  then  becomes 


5N 


3  Z  p 


7  x  IT 


(C-14) 


where  fringe  shift  count  5N  is  now  a  function  of  gas  constituents  and 
temperature.  This  equation  would  predict  gas  composition  or  concen¬ 
tration  in  an  isothermal  environment,  or  alternatively  temperature, 
if  the  gas  concentrations  did  not  change. 


Since  the  molar  refractivity  of  a  constituent  is  wavelength  - 
dependent  and  has  dispersion  values  that  differ  from  other  constituents, 
it  can  be  seen  that  to  measure  temperature  which  varies  and  gas 
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concentrations  which  also  vary  it  is  necessary  to  take  as  many  simul¬ 
taneous  interferograms  at  different  wavelengths  as  there  are  gas  com¬ 
ponents.  Expressions  similar  to  Equation  (C-  14)  for  each  fringe  count 
determination  are  all  solved  simultaneously  to  arrive  at  the  concen¬ 
trations  and  temperature  information. 

For  gas  mixtures  of  more  than  two  constituents  the  expressions 
become  quite  imposing.  As  earlier  indicated,  using  two  simultaneous 
wavelength  interference  recordings  to  accurately  measure  concentra¬ 
tions  or  temperature  for  air -water -vapor  mixture  requires  extremely 
precise  fringe  shift  measurements.  The  method  becomes  tenuous  where 
dispersion  values  for  the  gases  are  unknown,  *  where  the  dispersions  them¬ 
selves  are  very  similar,  or  where  many  gases  are  present  at  one  instant. 
Gas  analysis  can  be  more  assuredly  accomplished  by  improving  the 
measure  of  refractive  variances  or  using  variations  in  the  analysis  for 
whatever  conditions  exist. 

If  gases  arc  individually  released  as  a  function  cf  heat  loading  and 
time,  the  technique  reduces  to  the  two -wavelength  interferometric 
approach  where  the  concentration  of  gas  in  air  (or  in  any  background 
gaseous  field)  is  more  readily  measurable. 

The  interferometric  arrangement  used  to  record  simultaneous 
interference  fringe  fields  at  two  wavelengths  is  shown  in  Figure  C-l. 

White  light  division  and  recombination  of  wave  amplitude  is  accom¬ 
plished  in  the  normal  manner  using  the  Mach-Zehnder  interferometric 
model.  The  light,  however,  is  again  divided  by  an  equal;  path  -  splitting 
prism  in  which  each  beam  is  filtered  for  the  appropriate  wavelength 
desired.  For  an  increased  number  of  gas  components,  the  divided 
beams  could  conceivably  be  divided  again  in  similar  fashion  to  pro¬ 
vide  additional  wavelength  interference  fields. 

*The  International  Critical  Tables  ( i 9 3  3 )  has  dispersion  formulae  for 
various  gases.  Seemingly  little  information  has  been  published  since 
1933,  especially  for  refractivities  of  gases  and  their  dispersions  at 
higher  temperatures. 
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Figure  C-l.  Schematic  of  Mach-Zchnder  Two -Wavelength  Interferometer 


K.  Means  of  Improving  Fringe  Count  Accuracy 

Gas  analysis  using  this  method  requires  stringent  measurement 
accuracies  as  indicated  in  the  foregoing  discussion.  Fringe-count 
accuracy  can  be  improved  greatly  by  using  refined  techniques  in 
measuring  small  fractions  of  fringe  separation  itself  and  by  increasing 
the  length  through  which  the  light  is  refracted.  The  gas  environment 
in  which  tests  are  conducted  can  be  varied  to  provide  the  greatest  dis¬ 
persive  powers  for  identifying  and  quantifying  emitted  gases.  The 
wavelengths  of  light  used  in  a  two-wavelength  interferometric  study 
can  be  separated  widely  in  the  visible  spectrum  also  to  permit  as  great 
a  dispersion  difference  as  possible. 

Our  present  fringe  count  accuracy  is  approximately  0.  1  of  a  fringe 

separation  which  can  readily  be  improved  to  0.  01  fringe  by  methods 

a 

indicated  in  the  literature.  Werner  and  Leadon  suggested  a  method  to 
obtain  fringe  measurement  accuracy  to  l  '500  fringe  spacing.  Since 
variations  in  fringe  spacing  occur  from  the  undisturbed  to  the  dis¬ 
placed  fringe  fields,  this  accuracy  for  their  particular  measurements 
has  been  questioned  by  Winkler. 

■j 

Peters  and  Stroke  proposed  a  mea^s  «..f  electronic  location  and 
measurement  of  fringes  with  a  capability  of  discerning  a  displacement 
of  1/400  of  a  fringe.  With  irregularities  in  peripheral  servo  equipment 
to  restore  zerc  position,  an  accuracy  c  f  l  i  00  fringe  was  possible. 

12 

More  recently  a  photoelectric  scanning  device  described  by  Dew 
evaluated  inte rferograms  to  l.’iOO  fnrge  or  better.  The  same  accuracy 
was  claimed  by  Dipon^  for  a  device  used  to  measure  fringe  displace¬ 
ment  from  streak  photographs. 

In  Figure  C-2  it  can  be  nvt^d  that  fringe  count  accuracy  is  improved 
considerably  by  increasing  the  path  length  through  which  refraction  occurs. 
The  curves  can  be  read  in  terms  of  absolute  fringe  count  b  N  or  as  a 
measure  of  changes  in  fringe  count  A6N  since  a  linear  relationship  exists. 
The  corresponding  temperature  difference  associated  with  change  in 
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index  of  refraction  and  fringe  shift  number  (constant  pressure)  is  also 
shown.  From  this  plot  it  is  seen  that  for  our  presently  obtained  fringe 
count  accuracy  (0.  i  fringe)  and  for  our  path  length  (fabric  width)  of 
approximately  3  cm,  our  temperature  measurement  accuracy  is  some¬ 
what  less  than  2°C.  If  fringe  count  accuracy  can  be  improved  to  0.  01 
fringe  and  an  accuracy  of  0.  1°C  or  better  is  required,  a  path  length  of 
approximately  10  cm  must  be  used. 


Although  the  multi -wave length  technique  does  net  appear  to  lend 
itself  to  precise  measurements  of  water  vapor  in  an  air  environment 
(since  the  dispersion  values  of  water  vapor  are  not  considerably  dif¬ 
ferent  from  that  of  air),  the  technique  would  be  sensitive  in  detecting 
gases  having  high  dispersion  characteristics.  If  the  temperature  was 
measured  by  other  means  the  interferometer  would  predict  concen¬ 
tration  levels  accurately.  This,  of  course,  would  involve  only  one 
wavelength  of  light,  and  it  would  not  be  necessary  to  consider  sepa¬ 
ration  of  light  into  two  or  more  wavelengths  except  for  correlation 
measurements. 


It  can  be  recognized  that  if  the  wavelengths  of  light  used  in  the 
multi -wavelength  technique  are  at  extremes  in  the  visible  spectrum 
or,  if  possible  extended  well  into  the  infrared  at  one  extreme  and 
into  the  ultra  violet  at  the  other,  maximum  consideration  of  the  gas 
dispersion  characteristics  is  possible.  Howevt  r,  operational  limi¬ 
tation  of  the  Mach-Zehnder  interferometer  prevent  wide  separation  of 
wavelengths  in  utilizing  the  two  or  more  wavelength  technique. 

Usually,  optimum  interference  fringes  are  obtained  from  one  wave¬ 
length  of  light,  and  for  the  second  wavelength  an  effective  "out  of 
focus"  condition  exists.  This  results  in  part  because  of  the  finite 
size  of  the  interferometer  apferture  and  because  of  the  non-monochro- 
macity  of  the  light  source.  Thfe  optical  path  adjustment  elements  are 
generally  set  in  the  compensated  leg  of  the  interferometer  and  any 
scheme  to  compensate  optical  path  length  for  more  than  one  wavelength 
is  difficult. 
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Several  methods  were  devised  to  permit  compensation  of  two  or 
more  wavelengths.  Since  modification  of  our  Mach-Zehnder  inter¬ 
ferometer  would  be  necessary,  the  multi -wavelength  technique  was  not 
used  in  this  phase  our  program. 
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III,  MULTI-ENVIRONMENTAL  INTERFEROMETER  TECHNIQUE 

A.  General 

Separate  but  identical  fabric  heating  tests  in  different  environments 
provides  the  basis  for  this  technique.  The  unknown  variables  are  sepa¬ 
rated  by  solving  simultaneously  the  fringe  shift  expressions  for  the 
tests  in  the  various  environmental  gases. 


B.  Expressions  Governing  the  Two  or  More  Environmental 
Interferometer  Technique 


For  a  single  environment  the  basic  equation  of  the  Mach-Zehnder 

2 

interferometer  has  been  previously  worked  out  .  The  expression  is 
given  as 

C7 

5N  *  ^  (Pl  -p2)  (C- 15) 


where  8N  is  the  observed  fringe  shift,  C  the  Gladstone -Dale  constant  for 
the  particular  environmental  gas,  p  ^  the  gas  density  in  the  undisturbed 
portion  of  the  interference  fie'd,  p  ^  the  gas  density  at  the  point  under 
investigation,  Z  the  length  of  the  heated  fabric,  and  X  the  wavelength 
in  vacuum  of  the  light  used.  From  the  ideal  gas  law  we  can  express 
the  density  as 

p  *  pM/R't  (C- 16) 

where  p  is  the  pressure,  t  the  temperature,  M  the  molecular  weight  of 
the  gas,  and  R'the  universal  gas  constant.  Equation  (C-j5)  then  becomes 
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In  a  single  gas  environment 
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and  for  a  constant  pressure  condition 


When  a  foreign  gas  is  introduced  a  universal  Gladstone -Dale  constant 
cannot  be  used  and  Equation  (C-17)  thus  becomes 


5N 


(C  - 18) 


where  the  quantity  CM  is  defined  as  the  molecular  refractivity  of  a  gas. 

If  we  assume  that  a  gas  at  any  point  consists  of  a  mixture  of  the 
environmental  gas  and  of  the  foreign  gases  then 


(CM)2  *  £  ^  (CM).  (C-19) 

i 

where  X.  is  the  percentage  concentration  of  gas  component  i.  With  only 
one  species  of  foreign  gas  present,  for  example,  the  expression  may  be 
written 


6  N 


(CM)j 


XlCM)1  +  (1-X)  <CM)2 


(C-20) 


The  three  unknowns  in  the  above  expression,  molecular  refractivity 
(CM^i  temperature  ( )*  and  gas  concentration  can  be  determined 

by  solving  simultaneously  the  fringe  shift  values  (as  a  function  of  these 
variables)  obtained  from  tests  in  three  different  environmental  gases. 

Nitrogen,  having  essentially  the  same  molecular  weight  and  re> 
fractivity  as  air,  was  used  in  our  tests  principally  to  determine  fabric 
degradation  characteristics  that  may  be  different  from  oxidation  or 
combustion  in  air.  Argon  and  helium,  readily  obtainable  inert  gases, 
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wer*  chosen  to  represent  molecular  weights  above  and  below  that  of 
air  or  nitrogen.  As  shown  below,  air  and  argon  have  quite  similar 
molecular  refractivities  despite  the  fact  that  each  has  relatively  large 
differences  in  molecular  weight. 

Table  C-3.  Molecular  Weight  and  Refractivity 
Values  for  Various  Gases 

^If^vironmeni  Index  of  Gladstone-Dale  Molecular  Molecular 
Gas  Refraction  Constant  Weight  Refractivity 

n*  C  M  CM 


Air 

1. 000*936 

0. 2262 

29 

6.  56 

Nitrogen 

1. 0002998 

0.2397 

28 

6.  71 

Argon 

1. 0C0281 

0. 1583 

39.9 

6.32 

Helium 

1. 000036 

0. 1958 

4 

0.  783 

Water  Vapor 

1. 00025 

0.  3143 

18 

5.  66 

Carbon  Dioxide 

1.  00045 

0.2280 

44 

10.  03 

'’'Standard  Conditions:  0°C,  760  mm  Hg,  X  =  5461  A. 


C.  Theoretical  Fringe  Shift  Variations 

Theoretical  predictions  of  fringe  shift  values  versus  temperature 
for  various  gases  are  shown  in  Figure  C-3.  These  curves  display  the 
marked  effect  different  gases  have  on  fringe  count.  If  temperature 
measurement  sensitivity  is  required,  a  gas  environment  with  a  high 
molecuiar  refractivity  may  be  chosen  with  light  interference  occurring 
over  a  long  optical  path  (Figure  C-3  b). 

If  a  foreign  gas  is  introduced  across  the  interference  light  path, 
either  a  positive  or  ?.  negative  fringe  count  can  occur,  depending  on  the 
molecular  refractivity  of  both  foreign  and  environmental  gas. 
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The  sensitivity  of  this  multi-gas  interferometric  technique  for 
detecting  water  vapor,  its  concentration,  and  the  temperature  of  a 
heated  specimen,  can  be  considered  from  the  fringe  number  associated 
with  the  several  environmental  gases.  In  a  background  of  carbon 
dioxide,  for  example,  the  introduction  of  water  vapor  would  increase 
fringe  count,  whereas  in  an  environment  of  helium  the  fringe  count 
becomes  negative.  If  fringe  count  measurements  are  made  to  a  small 
fraction  of  a  fringe  separation,  then  the  two  equations  expressing  the 
fringe  values  for  tests  in  both  gases  would  normally  permit  considera¬ 
ble  accuracy  in  temperature  and  in  vapor  concentration  determinations. 
If  the  identity  of  the  foreign  gas  is  not  known,  three  or  more  tests  and 
at  least  three  equations  are  required. 

A  computerized  program  with  readout  similar  to  that  shown  in 
Figure  C-4  may  be  utilized  to  rapidly  determine  molecular  weight, 
concentration  and  temperature  from  measured  fringe  shift  values 
in  the  different  environments.  This  may  be  plotted  as  a  series  of 
graphs  for  more  rapid  determination  of  degradation  gases  in  the 
multi -environmental  test  technique.  This  method  has  several  pro¬ 
nounced  limitations.  It  would  be  inconceivable  that  exactly  identical 
tests  could  be  conducted  since  variances  exist  in  the  test  procedures; 
inhomogeneties  exist  in  the  arc  furnace  heating  rates,  in  the  fabric 
specimens  and  in  the  environmental  conditions.  Thus,  for  rapidly 
changing  events  interference  fringe  shift  correspondence  for 
separate  tests  in  different  environments  is  not  possible  and  analysis 
therefore  limited  to  quasi- steady  state  conditions.  A  further  limitation 
of  this  test  technique  occurs  since  the  fabric  thermal  boundary  layer 
profiles  are  different  for  tests  in  each  environment,  and  thus  the  free 
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convective  heat  transfer  characteristics  are  also  different.  It  is  ex¬ 
pected  that  the  mass  diffusion  processes  differ  for  the  separate  en¬ 
vironments  such  that  measurements  are  only  possible  at  the  gas- 
fabric  interfaces  where  zero  velocity  is  postulated.  This  condition  is 
not  altogether  true,  and  since  the  convective  heat  transfer  coefficients 
would  vary  for  each  test  as  a  function  of  time  and  heat  loading  the  actual 
temperature  of  the  fabric  and  possible  the  water  vapor  concentration 
would  be  somewhat  different  one  test  to  another. 

Another  question  arises  as  to  what  degradation  gases  are  present 
in  the  different  heating  tests.  The  surrounding  environment  would 
influence  to  some  extent  the  type  and  degree  of  fabric  degradation  that 
would  occur.  It  is  anticipated,  however,  that  the  behavior  of  a  heated 
fabric  would  be  somewhat  similar  for  tests  in  inert  environments.  Our 
subsequent  tests  were  made  in  helium,  nitrogen,  argon  and  CO^. 

Although  the  above  restrictions  are  quite  imposing,  several  im¬ 
portant  findings  are  possible.  For  a  water  saturated  fabric  a  temper¬ 
ature  plateau  exists  at  or  near  100*0  until  vaporization  is  complete. 
Thereafter  fabric  surface  temperatures  increase  with  time  until  more 
complete  degradation  or  combustion  occurs.  During  the  vaporization 
period  a  quasi-steady  state  condition  exists  permitting  ready  analysis 
of  concentrations.  The  fringe  shift  values  in  the  different  environ¬ 
ments  are  quite  different  and  an  extremely  sensitive  measure  of 
temperature  and  vapor  concentration  from  separate  environmental 
tests  is  possible. 

Although  the  vapor  concentrations  and  temperature  indications  would 
be  somewhat  different  in  the  combinational  evaluation  of  separate  environ¬ 
ment  tests,  a  consistent  variance  would  be  recognized  from  repeated 
experiments  and  comparisons  between  these  different  environmental 
combinations  (i.e.,  comparison  ox  nitrogen-helium  with  nitrogen-CO^) 
would  be  possible.  An  averaging  procedure  based  on  the  contributing 
effects  of  the  molecular  weights  and  refractive  index  constants  of 


environment  and  foreign  (water  vapor)  gases  could  be  conceived. 
Comparisons  were  made  from  our  tests  on  water -satm ated  fabrics; 
however,  no  averaging  procedure  was  devised. 

Extensive  heating  tests  were  made  in  different  environments  on  a 

single  wool  type  and  on  two  different  cotton  fabrics  under  both  dry  and 

wet  fabric  conditions.  For  the  water -saturated  tests,  concentrations 

of  water  vapor  at  the  fabric  surfaces  were  determined  for  a  single  heat 

2 

ing  load  (approximately  6  cal/cm  -sec).  To  facilitate  determinations 
of  molecular  weight,  concentration,  and  temperature  for  tests  of  the 
various  gases,  the  equations  expressing  interference  fringe  shift 
values  as  a  function  of  above  variables  were  programmed  for  solution 
on  a  Control  Data  1604  digital  computer. 
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BN 

(BN)c 

*1 

P 

Q 

R 

t 

x,  y,  z 

Y 

Z 

Z _ 

TC 

P 

X 

M 

MC 

Xx 

R' 

a 

v 

o 

e 

m 

P 

K,  k 


fringe  width 

Gladstone -Da  le  constant 

number  of  wavelengths  or  number  of  molecules  per 
unit  volume 

wavelength  difference  or  fringe  shift  count 
fringe  shift  referenced  to  compensated  optical  path 
refractive  index 

ambient  pressure  or  molar  refractivity 

optical  path  difference 

gas  constant 

temperature 

rectangular  coordinates 

fringe  shift 

length  of  heated  body 

length  of  test  chamber 

density 

wavelength  of  light  in  vacuum 

molecular  weight 

molecular  refractivity 

concentration  percentage,  mol  fraction 

universal  gas  constant 

polarizability  factor 

velocity  of  light  in  vacuum 

electron  charge 

electron  mass 

partial  pressure 

constants 


C-ilH 


Subscript 

1,2 


>$  f 


r 

s 

m 

wv 

c 
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location  or  condition 
reference 

standard  condition  (dry  air) 
moist  air 
water  vapor 
referenced  to  vacuum 
component 
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